Changeset 2672 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Apr 27, 2022, 2:37:49 PM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Fixed option to initialize chemistry in testphys1d (inichim_newstart call
arguments were wrong).
Turned inichim_newstart.F90 into a module.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (5 diffs)

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

@@ -32 +32 @@
       USE physiq_mod, ONLY: physiq
       USE read_profile_mod, ONLY: read_profile
+      use inichim_newstart_mod, only: inichim_newstart
       IMPLICIT NONE
 
@@ -107 +108 @@
 
 c   Various intermediate variables
-      INTEGER thermo
+      INTEGER flagthermo,flagh2o
       REAL zls
       REAL phi(nlayer),h(nlayer),s(nlayer)
@@ -117 +118 @@
       integer :: nq=1 ! number of tracers
       real :: latitude(1), longitude(1), cell_area(1)
+      ! dummy variables along "dynamics scalar grid"
+      real,allocatable :: qdyn(:,:,:,:),psdyn(:,:)
 
       character*2 str2
@@ -132 +135 @@
       INTEGER ierr0
       LOGICAL :: continu
+
+      logical :: present
 
 c=======================================================================
@@ -644 +649 @@
       if (photochem.or.callthermos) then
          write(*,*) 'Initializing chemical species'
-         ! thermo=0: initialize over all atmospheric layers
-         thermo=0
-         call inichim_newstart(q,psurf,sig,nq,latitude,longitude,
-     $   cell_area,thermo,qsurf)
+         ! flagthermo=0: initialize over all atmospheric layers
+         flagthermo=0
+         ! check if "h2o_vap" has already been initialized
+         ! (it has been if there is a "profile_h2o_vap" file around)
+         inquire(file="profile_h2o_vap",exist=present)
+         if (present) then
+           flagh2o=0 ! 0: do not initialize h2o_vap
+         else
+           flagh2o=1 ! 1: initialize h2o_vap in inichim_newstart
+         endif
+         
+         ! hack to accomodate that inichim_newstart assumes that
+         ! 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)
+         psdyn(1:2,1)=psurf
+         call inichim_newstart(ngrid,nq,qdyn,qsurf,psdyn,
+     $                         flagh2o,flagthermo)
+         q(1:llm,1:nq)=qdyn(1,1,1:llm,1:nq)
       endif