Changeset 1947 for trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

Timestamp:

Jun 13, 2018, 5:22:44 PM (6 years ago)

Author:

jvatant

Message:

Enables altitude-depending gravity field g(z) (glat->gzlat) in physics
+ Can be dangerous ( disagreement with dyn) but important (compulsive !) to have correct altitudes in the chemistry
+ Can be activated with eff_gz flag in callphys.def
-- JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

@@ -181 +181 @@
      call getin_p("Rmean", Rmean)
      write(*,*) "Rmean = ", Rmean
+     
+     write(*,*) "Compute effective altitude-dependent gravity field?"
+     eff_gz = .false.
+     call getin_p("eff_gz", eff_gz)
+     write(*,*) "eff_gz = ", eff_gz
          
 ! Test of incompatibility:
@@ -344 +349 @@
        stop
      endif
+     
+     ! sanity check warning
+     if (callchim.and.(.not.eff_gz)) then
+       print*,"WARNING : You run chemistry without effective altitude-dependent gravity field !!"
+       print*,"You will have wrong vertical photodissociations rates, that's crazy !!"
+       print*,"I let you continue but you should rather set eff_gz =.true. ..."
+       !stop
+     endif
+
 
      write(*,*)"Chemistry is computed every ichim", &