AWIESM 2.5 FESOM branch

README.md

@@ -28,6 +28,8 @@ References
 
 * **[Version with ICEPACK sea ice thermodynamics]** Zampieri, Lorenzo, Frank Kauker, Jörg Fröhle, Hiroshi Sumata, Elizabeth C. Hunke, and Helge Goessling. Impact of Sea-Ice Model Complexity on the Performance of an Unstructured-Mesh Sea-ice/ocean Model Under Different Atmospheric Forcings. Washington: American Geophysical Union, 2020. https://dx.doi.org/10.1002/essoar.10505308.1.
 
+* **[Version coupled with OpenIFS atmosphere]** Streffing, J., Sidorenko, D., Semmler, T., Zampieri, L., Scholz, P., Andrés-Martínez, M., et al ( 2022). AWI-CM3 coupled climate model: description and evaluation experiments for a prototype post-CMIP6 model. Geoscientific Model Development, 15. https://doi.org/10.5194/gmd-15-6399-2022
+
 Documentation for FESOM2
 ************************
 

config/namelist.dyn

@@ -19,5 +19,6 @@ use_freeslip = .false. ! Switch on free slip
 use_wsplit   = .false. ! Switch for implicite/explicte splitting of vert. velocity
 wsplit_maxcfl= 1.0     ! maximum allowed CFL criteria in vertical (0.5 < w_max_cfl < 1.) 
                        ! in older FESOM it used to be w_exp_max=1.e-3
+ldiag_KE=.false.       ! activates energy diagnostics
 /
 

config/namelist.io

@@ -2,7 +2,8 @@
 ldiag_solver     =.false.
 lcurt_stress_surf=.false.
 ldiag_curl_vel3  =.false.
-ldiag_energy     =.false.
+ldiag_Ri         =.false.
+ldiag_turbflux   =.false.
 ldiag_salt3D     =.false.
 ldiag_dMOC       =.false.
 ldiag_DVD        =.false.

config/namelist.io.nextGEMS

@@ -0,0 +1,53 @@
+&diag_list
+ldiag_solver     =.false.
+lcurt_stress_surf=.false.
+ldiag_curl_vel3  =.false.
+ldiag_Ri         =.false.
+ldiag_turbflux   =.false.
+ldiag_salt3D     =.false.
+ldiag_dMOC       =.false.
+ldiag_DVD        =.false.
+ldiag_forc       =.false.
+ldiag_extflds    =.false.
+/
+
+&nml_general
+io_listsize    =100 !number of streams to allocate. shallbe large or equal to the number of streams in &nml_list
+vec_autorotate =.false.
+lnextGEMS=.true.
+nlev_upper=10
+/
+
+! for sea ice related variables use_ice should be true, otherewise there will be no output
+! for 'curl_surf' to work lcurt_stress_surf must be .true. otherwise no output
+! for 'fer_C', 'bolus_u', 'bolus_v', 'bolus_w', 'fer_K' to work Fer_GM must be .true. otherwise no output
+! 'otracers' - all other tracers if applicable
+! for 'dMOC' to work ldiag_dMOC must be .true. otherwise no output
+&nml_list
+io_list =  'sst       ',1, 'm', 4,
+           'sss       ',1, 'm', 4,
+    	   'ssh       ',1, 'm', 4,
+           'uice      ',1, 'm', 4,
+           'vice      ',1, 'm', 4,
+           'a_ice     ',1, 'm', 4,
+           'm_ice     ',1, 'm', 4,
+           'm_snow    ',1, 'm', 4,
+           'MLD1      ',1, 'm', 4,
+           'MLD2      ',1, 'm', 4,
+           'MLD3      ',1, 'm', 4,
+           'tx_sur    ',1, 'm', 4,
+           'ty_sur    ',1, 'm', 4,
+           'temp      ',3, 'h', 4,
+           'salt      ',3, 'h', 8,
+           'N2        ',1, 'y', 4,
+           'Kv        ',1, 'y', 4,
+           'u         ',3, 'h', 4,
+           'v         ',3, 'h', 4,
+           'unod      ',1, 'y', 4,
+           'vnod      ',1, 'y', 4,
+           'w         ',3, 'h', 4,
+           'Av        ',1, 'y', 4,
+           'bolus_u   ',1, 'y', 4,
+           'bolus_v   ',1, 'y', 4,
+           'bolus_w   ',1, 'y', 4,
+/

setups/test_pi/setup.yml

@@ -29,8 +29,6 @@ namelist.ice:
         evp_rheol_steps: 120
 
 namelist.io:
-    diag_list:
-        ldiag_energy: False
     nml_list:
         io_list:
             "sst       ":

setups/test_souf/setup.yml

@@ -58,8 +58,6 @@ namelist.ice:
         evp_rheol_steps: 150
 
 namelist.io:
-    diag_list:
-        ldiag_energy: False
     nml_list:
         io_list:
             "sst       ":

src/CMakeLists.txt

@@ -32,9 +32,19 @@ if(ALEPH_CRAYMPICH_WORKAROUNDS)
 endif()
 if(ALBEDO_INTELMPI_WORKAROUNDS)
    add_compile_options(-DENABLE_ALBEDO_INTELMPI_WORKAROUNDS)
+   add_compile_options(-DDISABLE_PARALLEL_RESTART_READ)
+endif()
+
+if(${FESOM_PLATFORM_STRATEGY} STREQUAL  levante.dkrz.de )
+   message(STATUS "multithreading disabled for Levante") # multithreading suddenly produces an error, disable it until a fix is found. issue #413
+   option(DISABLE_MULTITHREADING "disable asynchronous operations" ON)
+elseif(${FESOM_PLATFORM_STRATEGY} STREQUAL  albedo )
+   message(STATUS "multithreading disabled for Albedo") # multithreading suddenly produces an error, disable it until a fix is found. issue #413
+   option(DISABLE_MULTITHREADING "disable asynchronous operations" ON)
+else()
+   option(DISABLE_MULTITHREADING "disable asynchronous operations" OFF)
 endif()
 
-option(DISABLE_MULTITHREADING "disable asynchronous operations" ON)
 option(ENABLE_OPENACC "compile with OpenACC support" OFF)
 set(NV_GPU_ARCH "cc80" CACHE STRING "GPU arch for nvfortran compiler (cc35,cc50,cc60,cc70,cc80,...)")
 

src/MOD_DYN.F90

@@ -67,7 +67,7 @@ MODULE MOD_DYN
     real(kind=WP), allocatable, dimension(:,:,:):: uvnode
 
     ! instant vertical vel arrays
-    real(kind=WP), allocatable, dimension(:,:)  :: w, w_e, w_i, cfl_z, fer_w
+    real(kind=WP), allocatable, dimension(:,:)  :: w, w_e, w_i, w_old, cfl_z, fer_w
 
     ! sea surface height arrays
     real(kind=WP), allocatable, dimension(:)    :: eta_n, d_eta, ssh_rhs, ssh_rhs_old
@@ -110,8 +110,21 @@ MODULE MOD_DYN
     logical                                     :: use_wsplit    = .false.
     ! maximum allowed CFL criteria in vertical (0.5 < w_max_cfl < 1.) 
     ! in older FESOM it used to be w_exp_max=1.e-3
-    real(kind=WP)                               :: wsplit_maxcfl= 1.0     
-
+    real(kind=WP)                               :: wsplit_maxcfl= 1.0
+    ! energy diagnostic part: will be computed inside the model ("hard integration"):
+    logical                                      :: ldiag_ke       = .true.
+    ! different contributions to velocity change. will be computed inside the code.
+    real(kind=WP), allocatable, dimension(:,:,:) :: ke_adv, ke_cor, ke_pre, ke_hvis, ke_vvis, ke_du2, ke_umean, ke_u2mean
+    real(kind=WP), allocatable, dimension(:,:)   :: ke_wind, ke_drag
+    ! same as above but multiplied by velocity. we need both for later computation of turbulent fluxes
+    real(kind=WP), allocatable, dimension(:,:,:) :: ke_adv_xVEL, ke_cor_xVEL, ke_pre_xVEL, ke_hvis_xVEL, ke_vvis_xVEL
+    real(kind=WP), allocatable, dimension(:,:)   :: ke_wind_xVEL, ke_drag_xVEL
+    real(kind=WP), allocatable, dimension(:,:)   :: ke_wrho         !we use pressure to compute (W*dens) as it appeares much easier to compute (P*dW) instead of (dP*w)
+    real(kind=WP), allocatable, dimension(:,:)   :: ke_dW, ke_Pfull !for later computation of turbulent fluxes from the term above
+    real(kind=WP), allocatable, dimension(:,:,:) :: ke_adv_AB, ke_cor_AB
+    real(kind=WP), allocatable, dimension(:,:,:) :: ke_rhs_bak
+    ! surface fields to compute APE generation
+    real(kind=WP), allocatable, dimension(:)     :: ke_J, ke_D, ke_G, ke_D2, ke_n0, ke_JD, ke_GD, ke_swA, ke_swB
     !___________________________________________________________________________
     contains
 #if defined(__PGI)

src/gen_model_setup.F90

@@ -6,7 +6,7 @@ subroutine setup_model(partit)
 !   use i_therm_param
   use g_forcing_param
   use g_config
-  use diagnostics, only: ldiag_solver,lcurt_stress_surf,lcurt_stress_surf, ldiag_energy, &
+  use diagnostics, only: ldiag_solver,lcurt_stress_surf,lcurt_stress_surf, ldiag_Ri, ldiag_TurbFlux, &
                          ldiag_dMOC, ldiag_DVD, diag_list
   use g_clock,     only: timenew, daynew, yearnew
   use g_ic3d 
@@ -28,10 +28,8 @@ subroutine setup_model(partit)
   read (fileunit, NML=geometry)
   read (fileunit, NML=calendar)
   read (fileunit, NML=run_config)
-!!$  read (fileunit, NML=machine)
   close (fileunit)
 
-!PS   partit%flag_debug=flag_debug
 
   ! ==========
   ! compute dt