Changeset 3067 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

Timestamp:

Oct 3, 2023, 11:21:28 AM (15 months ago)

Author:

jbclement

Message:

Mars PCM:
In 1D, 'q' has been converted from dimension (:,:) to (1,:,:) and 'q2' is now got through the module 'turb_mod'. It allows more generalization and to match dimension in the subroutines.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -5 +5 @@
 contains
 
-SUBROUTINE init_testphys1d(pem1d,ngrid,nlayer,odpref,nq,q,time,psurf,u,v,temp,startfiles_1D,therestart1D,therestartfi, &
-                           ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,q2,play,plev,   &
-                           latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
+SUBROUTINE init_testphys1d(pem1d,ngrid,nlayer,odpref,nq,q,time,psurf,u,v,temp,startfiles_1D,therestart1D, &
+                           therestartfi,ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,        &
+                           play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
 
 use ioipsl_getincom,          only: getin ! To use 'getin'
@@ -13 +13 @@
 use time_phylmdz_mod,         only: daysec, day_step, ecritphy, iphysiq
 use planete_h,                only: year_day, periheli, aphelie, peri_day, obliquit, emin_turb, lmixmin
-use surfdat_h,                only: albedodat, z0_default, z0, emissiv, emisice, albedice, iceradius, dtemisice, &
+use surfdat_h,                only: albedodat, z0_default, z0, emissiv, emisice, albedice, iceradius, dtemisice,      &
                                     zmea, zstd, zsig, zgam, zthe, hmons, summit, base, phisfi, watercaptag, watercap, &
                                     tsurf, emis, qsurf
@@ -37 +37 @@
 use mod_grid_phy_lmdz,        only: regular_lonlat
 use phys_state_var_init_mod,  only: phys_state_var_init
+use turb_mod,                 only: q2
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -48 +49 @@
 ! Arguments
 !=======================================================================
-integer,            intent(in) :: ngrid, nlayer
-real,               intent(in) :: odpref        ! DOD reference pressure (Pa)
-logical,            intent(in) :: pem1d         ! If initialization for the 1D PEM
+integer, intent(in) :: ngrid, nlayer
+real,    intent(in) :: odpref        ! DOD reference pressure (Pa)
+logical, intent(in) :: pem1d         ! If initialization for the 1D PEM
 
 integer, intent(inout) :: nq
 
-real, dimension(:,:), allocatable, intent(out) :: q     ! tracer mixing ratio (e.g. kg/kg)
-real,                              intent(out) :: time  ! time (0<time<1; time=0.5 at noon)
-real,                              intent(out) :: psurf ! Surface pressure
-real, dimension(nlayer),           intent(out) :: u, v  ! zonal, meridional wind
-real, dimension(nlayer),           intent(out) :: temp  ! temperature at the middle of the layers
-logical,                           intent(out) :: startfiles_1D, therestart1D, therestartfi ! Use of starting files for 1D
-integer,                           intent(out) :: ndt
-real,                              intent(out) :: ptif, pks
-real,                              intent(out) :: dttestphys                   ! testphys1d timestep
-real, dimension(:),   allocatable, intent(out) :: zqsat                        ! useful for (atm_wat_profile=2)
-real, dimension(:,:), allocatable, intent(out) :: dq, dqdyn                    ! Physical and dynamical tandencies
-integer,                           intent(out) :: day0                         ! initial (sol ; =0 at Ls=0) and final date
-real,                              intent(out) :: day                          ! date during the run
-real,                              intent(out) :: gru, grv                     ! prescribed "geostrophic" background wind
-real, dimension(nlayer),           intent(out) :: w                            ! "Dummy wind" in 1D
-real, dimension(nlayer + 1),       intent(out) :: q2                           ! Turbulent Kinetic Energy
-real, dimension(nlayer),           intent(out) :: play                         ! Pressure at the middle of the layers (Pa)
-real, dimension(nlayer + 1),       intent(out) :: plev                         ! intermediate pressure levels (pa)
-real, dimension(1),                intent(out) :: latitude, longitude, cell_area
-real,                              intent(out) :: atm_wat_profile, atm_wat_tau ! Force atmospheric water profiles
+real, dimension(:,:,:), allocatable, intent(out) :: q     ! tracer mixing ratio (e.g. kg/kg)
+real,                                intent(out) :: time  ! time (0<time<1; time=0.5 at noon)
+real,                                intent(out) :: psurf ! Surface pressure
+real, dimension(nlayer),             intent(out) :: u, v  ! zonal, meridional wind
+real, dimension(nlayer),             intent(out) :: temp  ! temperature at the middle of the layers
+logical,                             intent(out) :: startfiles_1D, therestart1D, therestartfi ! Use of starting files for 1D
+integer,                             intent(out) :: ndt
+real,                                intent(out) :: ptif, pks
+real,                                intent(out) :: dttestphys                   ! testphys1d timestep
+real, dimension(:),     allocatable, intent(out) :: zqsat                        ! useful for (atm_wat_profile=2)
+real, dimension(:,:,:), allocatable, intent(out) :: dq, dqdyn                    ! Physical and dynamical tandencies
+integer,                             intent(out) :: day0                         ! initial (sol ; =0 at Ls=0) and final date
+real,                                intent(out) :: day                          ! date during the run
+real,                                intent(out) :: gru, grv                     ! prescribed "geostrophic" background wind
+real, dimension(nlayer),             intent(out) :: w                            ! "Dummy wind" in 1D
+real, dimension(nlayer),             intent(out) :: play                         ! Pressure at the middle of the layers (Pa)
+real, dimension(nlayer + 1),         intent(out) :: plev                         ! intermediate pressure levels (pa)
+real, dimension(1),                  intent(out) :: latitude, longitude, cell_area
+real,                                intent(out) :: atm_wat_profile, atm_wat_tau ! Force atmospheric water profiles
 
 !=======================================================================
@@ -237 +237 @@
 
 ! allocate arrays:
-allocate(tname(nq),q(nlayer,nq),zqsat(nlayer))
-allocate(dq(nlayer,nq),dqdyn(nlayer,nq),tnom_transp(nq))
+allocate(tname(nq),q(1,nlayer,nq),zqsat(nlayer))
+allocate(dq(1,nlayer,nq),dqdyn(1,nlayer,nq),tnom_transp(nq))
 
 ! read tracer names from file traceur.def
@@ -289 +289 @@
 write(*,*) 'nqfils=',nqfils
 
-
-
 #ifdef CPP_XIOS
     call init_physics_distribution(regular_lonlat,4,1,1,1,nlayer,COMM_LMDZ)
@@ -455 +453 @@
     do iq = 1,nq
         open(3,file = start1Dname,status = "old",action = "read")
-        read(3,*) header, qsurf(1,iq,1),(q(ilayer,iq), ilayer = 1,nlayer)
+        read(3,*) header, qsurf(1,iq,1),(q(1,ilayer,iq), ilayer = 1,nlayer)
         if (trim(tname(iq)) /= trim(header)) then
             write(*,*) 'Tracer names not compatible for initialization with "'//trim(start1Dname)//'"!'
@@ -554 +552 @@
 w = 0. ! default: no vertical wind
 
-! Initialize turbulente kinetic energy
+! Initialize turbulent kinetic energy
 q2 = 0.
 
@@ -690 +688 @@
     ! q & psurf arrays are on the dynamics scalar grid
     allocate(qdyn(2,1,llm,nq),psdyn(2,1))
-    qdyn(1,1,1:llm,1:nq) = q(1:llm,1:nq)
+    qdyn(1,1,1:llm,1:nq) = q(1,1:llm,1:nq)
     psdyn(1:2,1) = psurf
     call inichim_newstart(ngrid,nq,qdyn,qsurf(1,:,1),psdyn,flagh2o,flagthermo)
-    q(1:llm,1:nq) = qdyn(1,1,1:llm,1:nq)
+    q(1,1:llm,1:nq) = qdyn(1,1,1:llm,1:nq)
 endif