Changeset 2672

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/inichim_newstart.F90 (modified) (3 diffs)
• libf/phymars/dyn1d/testphys1d.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3701 +3701 @@
 == 25/04/2022 == AB
 Put the CO2 ice scavenging ratio for dust and water ice to 20 instead of 100, to prevent drying of the water cycle.
+
+== 27/04/2022 == EM
+Fixed option to initialize chemistry in testphys1d (inichim_newstart call
+arguments were wrong).
+Turned inichim_newstart.F90 into a module.

trunk/LMDZ.MARS/libf/aeronomars/inichim_newstart.F90

@@ +1 @@
+module inichim_newstart_mod
+
+implicit none
+
+contains
+
       subroutine inichim_newstart(ngrid, nq, pq, qsurf, ps, &
                                   flagh2o, flagthermo)
@@ -49 +55 @@
 ! outputs :
 
-      real,intent(out) :: pq(nbp_lon+1,nbp_lat,nbp_lev,nq)  ! advected fields, ie chemical species
+      real,intent(inout) :: pq(nbp_lon+1,nbp_lat,nbp_lev,nq)  ! advected fields, ie chemical species
       real,intent(out) :: qsurf(ngrid,nq)     ! surface values (kg/m2) of tracers
 
@@ -810 +816 @@
       deallocate(vmrint)
 
-      end
+      end subroutine inichim_newstart
+
+end module inichim_newstart_mod

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