Changeset 3100 for trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd

Timestamp:

Oct 25, 2023, 10:34:54 AM (15 months ago)

Author:

bhatnags

Message:

Generic-PCM

This commit updates the slab ocean module to a parallelisable dynamic slab ocean module. This is particularly relevant if you want to be able to use oceanic heat transport in parallel mode.

It has the following features:

(a) Computes sea ice creation and evolution.
(b) Snow has thermodynamic properties.
(c) Computes oceanic horizontal transport (diffusion & surface-wind driven Ekman transport).
(d) Can be used in parallel mode.

Required callphys.def flags:
The slab ocean and its dependencies can be activated with the following flags (already added to deftank):
## Ocean options
## ~
# Model slab-ocean (Main flag for slab ocean)
ok_slab_ocean = .true.
# The following flags can only be set to true if ok_slab_ocean is true
# Ekman transport
slab_ekman = .true.
# Ekman zonal advection
slab_ekman_zonadv = .true.
# Horizontal diffusion (default coef_hdiff=25000., can be changed)
slab_hdiff = .true.
# Slab-ocean timestep (in physics timesteps)
cpl_pas = 1
# Gent-McWilliams? Scheme (can only be true if slab_ekman is true)
slab_gm = .true.

Notes:
In the current state, the model crashes if moistadjustment = .true. Unsure whether this is due to the slab or is an inherent issue with moistadj (under investigation).

SB and EM

Location:

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd

Files:

• ini_archive.F (modified) (3 diffs)
• iniphysiq_mod.F90 (modified) (2 diffs)
• lect_start_archive.F (modified) (6 diffs)
• newstart.F (modified) (1 diff)
• start2archive.F (modified) (5 diffs)
• write_archive.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/ini_archive.F

@@ -36 +36 @@
 
       USE comsoil_h
-      USE slab_ice_h, only: noceanmx
+!      USE slab_ice_h, only: noceanmx
+      USE ocean_slab_mod, ONLY: nslay
 !      use control_mod
       USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt
@@ -81 +82 @@
       INTEGER idim_tim
       INTEGER idim_nsoilmx ! "subsurface_layers" dimension ID #
-      INTEGER idim_noceanmx ! "ocean_layers" dimension ID #
+      INTEGER idim_nslay ! "ocean_layers" dimension ID #
       INTEGER nid,nvarid
       real sig_s(llm),s(llm)
@@ -164 +165 @@
       ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
       ierr = NF_DEF_DIM (nid,"subsurface_layers",nsoilmx,idim_nsoilmx)
-      ierr = NF_DEF_DIM (nid,"ocean_layers",noceanmx,idim_noceanmx)
+      ierr = NF_DEF_DIM (nid,"ocean_layers",nslay,idim_nslay)
 
       ierr = NF_DEF_DIM (nid,"index", length, idim_index)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/iniphysiq_mod.F90

@@ -11 +11 @@
 
 use control_mod, only: nday 
-use surf_heat_transp_mod, only: ini_surf_heat_transp
+!use surf_heat_transp_mod, only: ini_surf_heat_transp
+USE slab_heat_transp_mod, ONLY: ini_slab_transp_geom
 use infotrac, only : nqtot ! number of advected tracers
 use planete_mod, only: ini_planete_mod
@@ -83 +84 @@
 call getin_p("ok_slab_ocean",ok_slab_ocean)
 if (ok_slab_ocean) then
-  call ini_surf_heat_transp(ip1jm,ip1jmp1,unsairez,fext,unsaire, &
-                            cu,cuvsurcv,cv,cvusurcu,aire,apoln,apols, &
-                            aireu,airev)
+!  call ini_surf_heat_transp(ip1jm,ip1jmp1,unsairez,fext,unsaire, &
+!                            cu,cuvsurcv,cv,cvusurcu,aire,apoln,apols, &
+!                            aireu,airev)
+   CALL ini_slab_transp_geom(ip1jm,ip1jmp1,unsairez,fext,unsaire,&
+                                  cu,cuvsurcv,cv,cvusurcu, &
+                                  aire,apoln,apols, &
+                                  aireu,airev,rlatvdyn)
 endif
 

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F

@@ -9 +9 @@
       USE tracer_h, ONLY: igcm_co2_ice
       USE infotrac, ONLY: tname, nqtot
-      USE slab_ice_h, ONLY: noceanmx
+!      USE slab_ice_h, ONLY: noceanmx
+      USE ocean_slab_mod, ONLY: nslay
       USE callkeys_mod, ONLY: ok_slab_ocean
       USE comvert_mod, ONLY: ap,bp,aps,bps,preff
@@ -68 +69 @@
       REAL,INTENT(OUT) :: emis(ngrid)
       REAL,INTENT(OUT) :: q2(ngrid,llm+1),qsurf(ngrid,nqtot)
-      REAL,INTENT(OUT) :: tslab(ngrid,noceanmx)
+      REAL,INTENT(OUT) :: tslab(ngrid,nslay)
       REAL ,INTENT(OUT) ::rnat(ngrid),pctsrf_sic(ngrid)
       REAL,INTENT(OUT) :: tsea_ice(ngrid),sea_ice(ngrid)
@@ -98 +99 @@
       real emisS(iip1,jjp1)
       REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
-      real tslabS(iip1,jjp1,noceanmx)
+      real tslabS(iip1,jjp1,nslay)
       real pctsrf_sicS(iip1,jjp1),tsea_iceS(iip1,jjp1)
       real rnatS(iip1,jjp1), sea_iceS(iip1,jjp1)
@@ -287 +288 @@
       allocate(mlayerold(nsoilold))
       allocate(qsurfold(imold+1,jmold+1,nqtot))
-      allocate(tslabold(imold+1,jmold+1,noceanmx))
+      allocate(tslabold(imold+1,jmold+1,nslay))
       allocate(rnatold(imold+1,jmold+1))
       allocate(pctsrf_sicold(imold+1,jmold+1))
@@ -684 +685 @@
       if(ok_slab_ocean) then
       start=(/1,1,1,memo/)
-      count=(/imold+1,jmold+1,noceanmx,1/)
+      count=(/imold+1,jmold+1,nslay,1/)
 
        ierr=NF_INQ_VARID(nid,"tslab",nvarid)
@@ -1298 +1299 @@
 c 6.3   Slab Ocean :
 c-----------------------------------------------------------------------
-      call interp_horiz (tslabold,tslabs,imold,jmold,iim,jjm,noceanmx,
-     &                   rlonuold,rlatvold,rlonu,rlatv)
-      call gr_dyn_fi (noceanmx,iim+1,jjm+1,ngrid,tslabs,tslab)
+      call interp_horiz (tslabold,tslabs,imold,jmold,iim,jjm,nslay,
+     &                   rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (nslay,iim+1,jjm+1,ngrid,tslabs,tslab)
 
       call interp_horiz (rnatold,rnats,imold,jmold,iim,jjm,1,

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F

@@ -30 +30 @@
       use phyredem, only: physdem0, physdem1
       use iostart, only: open_startphy
-      use slab_ice_h, only:noceanmx
+!      use slab_ice_h, only:noceanmx
+      USE ocean_slab_mod, ONLY: nslay
       use filtreg_mod, only: inifilr
       USE mod_const_mpi, ONLY: COMM_LMDZ

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/start2archive.F

@@ -22 +22 @@
       USE comsoil_h
       
-      use slab_ice_h, only: noceanmx
+!      use slab_ice_h, only: noceanmx
+      USE ocean_slab_mod, ONLY: nslay
       USE ioipsl_getincom, only: getin
       USE planete_mod, only: year_day
@@ -84 +85 @@
 !     added by BC for slab ocean
       REAL rnat(ngridmx),pctsrf_sic(ngridmx),sea_ice(ngridmx)
-      REAL tslab(ngridmx,noceanmx),tsea_ice(ngridmx)
+      REAL, ALLOCATABLE :: tslab(:,:)
+      REAL tsea_ice(ngridmx)
 
 
@@ -104 +106 @@
 !     added by BC for slab ocean
       REAL rnatS(ip1jmp1),pctsrf_sicS(ip1jmp1),sea_iceS(ip1jmp1)
-      REAL tslabS(ip1jmp1,noceanmx),tsea_iceS(ip1jmp1)
+      REAL, ALLOCATABLE :: tslabS(:,:)
+      REAL tsea_iceS(ip1jmp1)
 
 !     For non-orographic GW
@@ -156 +159 @@
       allocate(qsurf(ngridmx,nqtot))
       allocate(qsurfS(ip1jmp1,nqtot))
+      allocate(tslab(ngridmx,nslay)) !Added by SB for slab ocean
+      allocate(tslabS(ip1jmp1,nslay)) !Added by SB for slab ocean
 ! other array allocations:
 !      call ini_comsoil_h(ngridmx) ! done via iniphysiq
@@ -356 +361 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,tsea_ice,tsea_iceS)
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,sea_ice,sea_iceS)
-      call gr_fi_dyn(noceanmx,ngridmx,iip1,jjp1,tslab,tslabS)
+      call gr_fi_dyn(nslay,ngridmx,iip1,jjp1,tslab,tslabS)
 
       call gr_fi_dyn(llm,ngridmx,iip1,jjp1,du_nonoro_gwd,du_nonoro_gwdS)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/write_archive.F

@@ -33 +33 @@
 
       use comsoil_h, only: nsoilmx
-      use slab_ice_h, only: noceanmx
+!      use slab_ice_h, only: noceanmx
+      USE ocean_slab_mod, ONLY: nslay
 
       implicit none
@@ -195 +196 @@
         edges(1)=iip1
         edges(2)=jjp1
-        edges(3)=noceanmx
+        edges(3)=nslay
         edges(4)=1
 #ifdef NC_DOUBLE