Changeset 2789 for trunk/LMDZ.MARS/libf

Timestamp:

Aug 31, 2022, 6:30:47 PM (2 years ago)

Author:

jnaar

Message:

1d flag "atm_wat_profile" can be used to force atmospheric water vapor profile :

• atm_wat_profile = 0 (default case) : free vapor & ice profiles
• atm_wat_profile = 1 : dry atmosphere for both vapor and ice
• atm_wat_profile = 2 : wet water vapor profile (0.7*qsat for now) and no ice

Also corrected mistake in README file where EM was responsable for 01/07/2022 commit (but it was me).
JN

File:

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

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

@@ -15 +15 @@
       use phyredem, only: physdem0,physdem1
       use geometry_mod, only: init_geometry
+      use watersat_mod, only: watersat
+      use tracer_mod, only: igcm_h2o_vap,igcm_h2o_ice
       use planete_h, only: year_day, periheli, aphelie, peri_day,
      &                     obliquit, emin_turb, lmixmin
@@ -128 +130 @@
       REAL halfaxe, excentric, Lsperi
       Logical paleomars
+c   JN : Force atmospheric water profiles
+      INTEGER atm_wat_profile
+      REAL,ALLOCATABLE :: zqsat(:) ! useful for (atm_wat_profile=2)
+
 
 c   MVals: isotopes as in the dynamics (CRisi)
@@ -246 +252 @@
         allocate(tname(nq))
         allocate(q(nlayer,nq))
+        allocate(zqsat(nlayer))
         allocate(qsurf(nq))
         allocate(dq(nlayer,nq))
@@ -678 +685 @@
       write(*,*) " watercaptag = ",watercaptag
       
+c Check if the atmospheric water profile is specified
+c ---------------------------------------------------
+      ! Adding an option to force atmospheric water values JN
+      atm_wat_profile=0 ! Default: no water ice reservoir
+      print *,'Force atmospheric water vapor profile ?'
+      call getin("atm_wat_profile",atm_wat_profile)
+      write(*,*) "atm_wat_profile = ", atm_wat_profile
+      if (atm_wat_profile.EQ.0) then
+        write(*,*) "Free atmospheric water vapor profile"
+      else if (atm_wat_profile.EQ.1) then
+        write(*,*) "Dry atmospheric water vapor profile"
+      else if (atm_wat_profile.EQ.2) then
+        write(*,*) "Wet atmospheric water vapor profile"
+      endif 
+
 
 c    Initialization for GRADS outputs in "g1d.dat" and "g1d.ctl"
@@ -742 +764 @@
       ENDDO
 
+c       Force atmospheric water if asked
+c       Added "atm_wat_profile" flag (JN)
+c       -------------------------------------
+        call watersat(nlayer,temp,play,zqsat)
+        DO iq = 1, nq
+          IF ((iq.eq.igcm_h2o_vap).and.(atm_wat_profile.eq.1)) THEN
+            DO ilayer=1,nlayer
+             q(ilayer,igcm_h2o_vap)=0.
+c             write(*,*) "atm_wat_profile dry"
+c             write(*,*) "iq (igcm_h2o_vap)", iq
+            ENDDO! ilayer=1,nlayer
+          ELSE IF((iq.eq.igcm_h2o_vap).and.(atm_wat_profile.eq.2)) THEN
+            DO ilayer=1,nlayer
+             q(ilayer,igcm_h2o_vap)=zqsat(ilayer)*0.7
+c             write(*,*) "atm_wat_profile wet"
+c             write(*,*) "iq (igcm_h2o_vap)", iq
+c             q(ilayer,igcm_h2o_vap)=zqsat(ilayer)-0.9*zqsat(ilayer)
+            ENDDO! ilayer=1,nlayer
+          ELSE IF ((iq.eq.igcm_h2o_ice).and.(atm_wat_profile.ne.0)) THEN
+            DO ilayer=1,nlayer
+             q(ilayer,igcm_h2o_ice)=0.
+c             write(*,*) "atm_wat_profile : reset ice"
+c             write(*,*) "iq (igcm_h2o_ice)", iq
+            ENDDO! ilayer=1,nlayer
+          ENDIF  !((iq.eq.igcm_h2o_vap).and.(atm_wat_profile.eq.2)) THEN
+        ENDDO
+
+
+!      write(*,*) "testphys1d avant q", q(1,:)
 c       call physics
 c       --------------------
-!      write(*,*) "testphys1d avant q", q(1,:)
       CALL physiq (1,llm,nq,
      ,     firstcall,lastcall,