Changeset 3067 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.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/testphys1d.F90 (modified) (6 diffs)

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

@@ -1 +1 @@
 PROGRAM testphys1d
 
-use comsoil_h,               only: inertiedat, inertiesoil, nsoilmx, tsoil
-use surfdat_h,               only: albedodat, perenial_co2ice, watercap, tsurf, emis, qsurf
-use comslope_mod,            only: def_slope, subslope_dist
-use phyredem,                only: physdem0, physdem1
-use watersat_mod,            only: watersat
-use tracer_mod,              only: igcm_h2o_vap, igcm_h2o_ice, noms
-use comcstfi_h,              only: pi, rad, omeg, g, mugaz, rcp, r, cpp
-use time_phylmdz_mod,        only: daysec, day_step
-use dimradmars_mod,          only: tauvis, totcloudfrac, albedo
-use dust_param_mod,          only: tauscaling
-use comvert_mod,             only: ap, bp, aps, bps, pa, preff, sig
-use physiq_mod,              only: physiq
-use phyetat0_mod,            only: phyetat0
-use write_output_mod,        only: write_output
-use init_testphys1d_mod,     only: init_testphys1d
+use comsoil_h,           only: inertiedat, inertiesoil, nsoilmx, tsoil
+use surfdat_h,           only: albedodat, perenial_co2ice, watercap, tsurf, emis, qsurf
+use comslope_mod,        only: def_slope, subslope_dist
+use phyredem,            only: physdem0, physdem1
+use watersat_mod,        only: watersat
+use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, noms
+use comcstfi_h,          only: pi, rad, omeg, g, mugaz, rcp, r, cpp
+use time_phylmdz_mod,    only: daysec, day_step
+use dimradmars_mod,      only: tauvis, totcloudfrac, albedo
+use dust_param_mod,      only: tauscaling
+use comvert_mod,         only: ap, bp, aps, bps, pa, preff, sig
+use physiq_mod,          only: physiq
+use turb_mod,            only: q2
+use write_output_mod,    only: write_output
+use init_testphys1d_mod, only: init_testphys1d
 ! Mostly for XIOS outputs:
-use mod_const_mpi,           only: init_const_mpi
-use parallel_lmdz,           only: init_parallel
+use mod_const_mpi,       only: init_const_mpi
+use parallel_lmdz,       only: init_parallel
 
 implicit none
@@ -58 +58 @@
 ! Declarations
 !--------------------------------------------------------------
-integer, parameter                :: ngrid = 1     ! (2+(jjm-1)*iim - 1/jjm)
-integer, parameter                :: nlayer = llm
-real, parameter                   :: odpref = 610. ! DOD reference pressure (Pa)
-integer                           :: unitstart     ! unite d'ecriture de "startfi"
-integer                           :: ndt, ilayer, ilevel, isoil, idt, iq
-logical                           :: firstcall, lastcall
-integer                           :: day0          ! initial (sol ; =0 at Ls=0)
-real                              :: day           ! date during the run
-real                              :: time          ! time (0<time<1 ; time=0.5 a midi)
-real                              :: dttestphys    ! testphys1d timestep
-real, dimension(nlayer)           :: play          ! Pressure at the middle of the layers (Pa)
-real, dimension(nlayer + 1)       :: plev          ! intermediate pressure levels (pa)
-real                              :: psurf         ! Surface pressure
-real, dimension(nlayer)           :: u, v          ! zonal, meridional wind
-real                              :: gru, grv      ! prescribed "geostrophic" background wind
-real, dimension(nlayer)           :: temp          ! temperature at the middle of the layers
-real, dimension(:,:), allocatable :: q             ! tracer mixing ratio (e.g. kg/kg)
-real, dimension(1)                :: wstar = 0.    ! Thermals vertical velocity
-real, dimension(nlayer + 1)       :: q2            ! Turbulent Kinetic Energy
+integer, parameter                  :: ngrid = 1     ! (2+(jjm-1)*iim - 1/jjm)
+integer, parameter                  :: nlayer = llm
+real, parameter                     :: odpref = 610. ! DOD reference pressure (Pa)
+integer                             :: unitstart     ! unite d'ecriture de "startfi"
+integer                             :: ndt, ilayer, ilevel, isoil, idt, iq
+logical                             :: firstcall, lastcall
+integer                             :: day0          ! initial (sol ; =0 at Ls=0)
+real                                :: day           ! date during the run
+real                                :: time          ! time (0<time<1 ; time=0.5 a midi)
+real                                :: dttestphys    ! testphys1d timestep
+real, dimension(nlayer)             :: play          ! Pressure at the middle of the layers (Pa)
+real, dimension(nlayer + 1)         :: plev          ! intermediate pressure levels (pa)
+real                                :: psurf         ! Surface pressure
+real, dimension(nlayer)             :: u, v          ! zonal, meridional wind
+real                                :: gru, grv      ! prescribed "geostrophic" background wind
+real, dimension(nlayer)             :: temp          ! temperature at the middle of the layers
+real, dimension(:,:,:), allocatable :: q             ! tracer mixing ratio (e.g. kg/kg)
+real, dimension(1)                  :: wstar = 0.    ! Thermals vertical velocity
 
 ! Physical and dynamical tandencies (e.g. m.s-2, K/s, Pa/s)
-real, dimension(nlayer)           :: du, dv, dtemp, dudyn, dvdyn, dtempdyn
-real, dimension(1)                :: dpsurf
-real, dimension(:,:), allocatable :: dq, dqdyn
+real, dimension(nlayer)             :: du, dv, dtemp, dudyn, dvdyn, dtempdyn
+real, dimension(1)                  :: dpsurf
+real, dimension(:,:,:), allocatable :: dq, dqdyn
 
 ! Various intermediate variables
@@ -115 +114 @@
 !call initcomgeomphy
 
-call init_testphys1d(.false.,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)
+call init_testphys1d(.false.,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)
 
 ! Write a "startfi" file
@@ -171 +170 @@
         ! If atmospheric water is monitored
             if (atm_wat_tau < 0.) then ! Prescribed atm_wat_profile: wet if >0, dry if =0
-                q(:,igcm_h2o_vap) = min(zqsat(:),atm_wat_profile*g/psurf)
-                q(:,igcm_h2o_ice) = 0. ! reset h2o ice
+                q(1,:,igcm_h2o_vap) = min(zqsat(:),atm_wat_profile*g/psurf)
+                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
             else ! Relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
-                q(:,igcm_h2o_vap) = atm_wat_profile*g/psurf + (q(:,igcm_h2o_vap) - atm_wat_profile*g/psurf)*dexp(-dttestphys/atm_wat_tau)
-                q(:,igcm_h2o_vap) = min(zqsat(:),q(:,igcm_h2o_vap))
-                q(:,igcm_h2o_ice) = 0. ! reset h2o ice
+                q(1,:,igcm_h2o_vap) = atm_wat_profile*g/psurf + (q(1,:,igcm_h2o_vap) - atm_wat_profile*g/psurf)*dexp(-dttestphys/atm_wat_tau)
+                q(1,:,igcm_h2o_vap) = min(zqsat(:),q(1,:,igcm_h2o_vap))
+                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
             endif
         endif
@@ -195 +194 @@
     !enddo
     ! For some tests: No coriolis force at equator
-    !if(latitude(1) == 0.) then
+    !if (latitude(1) == 0.) then
     du(:) = du(:) + (gru - u(:))/1.e4
     dv(:) = dv(:) + (grv - v(:))/1.e4
@@ -222 +221 @@
 
     ! Increment tracers
-    q(:,:) = q(:,:) + dttestphys*dq(:,:)
+    q(1,:,:) = q(1,:,:) + dttestphys*dq(1,:,:)
 enddo ! End of time stepping loop (idt=1,ndt)