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

Timestamp:

Apr 26, 2023, 3:08:01 PM (20 months ago)

Author:

evos

Message:

adding the options to write profiles at the end of the 1D simulation (write_prof=.true.)

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• testphys1d.F (modified) (3 diffs)
• write_profile.F90 (added)

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

@@ -18 +18 @@
       use geometry_mod, only: init_geometry
       use watersat_mod, only: watersat
-      use tracer_mod, only: igcm_h2o_vap,igcm_h2o_ice,igcm_co2
+      use tracer_mod, only: igcm_h2o_vap,igcm_h2o_ice,igcm_co2,noms
       use planete_h, only: year_day, periheli, aphelie, peri_day,
      &                     obliquit, emin_turb, lmixmin
@@ -87 +87 @@
       INTEGER ilayer,ilevel,isoil,idt,iq
       LOGICAl firstcall,lastcall
+      LOGICAL write_prof
 c
       real,parameter :: odpref=610. ! DOD reference pressure (Pa)
@@ -987 +988 @@
 
       ENDDO   ! of idt=1,ndt ! end of time stepping loop
+      
+      ! update the profiles files at the end of the run
+      write_prof=.false.
+      call getin("write_prof",write_prof)
+      IF (write_prof.eq..true.) then
+              DO iq = 1, nq
+                call writeprofile(nlayer,q(:,iq),noms(iq),iq)
+              ENDDO
+      ENDIF
+
 
 c    ========================================================