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

Timestamp:

Apr 29, 2019, 2:52:42 PM (6 years ago)

Author:

jvatant

Message:

+ Add diagnostics of optical thickness, in 1D, if 'diagdtau' key is activated, it

outputs dtaui/v(altitude) in diagfi.nc for every narrowband (could be done with one var
but would require to be able to have writediag in 5D

--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90 (modified) (1 diff)

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

@@ -336 +336 @@
      call getin_p("specOLR",specOLR)
      write(*,*)" specOLR = ",specOLR
+
+     write(*,*)"Output diagnostic optical thickness?"
+     diagdtau=.false.
+     call getin_p("diagdtau",diagdtau)
+     write(*,*)" diagdtau = ",diagdtau
+     ! sanity check
+     if (diagdtau.and.(ngrid.ne.1)) then
+       print*,"Diagnostic optical thickness can be output in 1D only !"
+       print*,"Start again with diagdtau=.false."
+       stop
+     endif
+
 
      write(*,*)"Uniform absorption in radiative transfer?"