Changeset 299

Timestamp:

Sep 21, 2011, 12:12:23 PM (13 years ago)

Author:

aslmd

Message:

LMDZ.MARS [associated to previous commit]

21/09/11 == AS

--> Added MESOINI precompiler flag so that all fields needed to initialize the mesoscale model

are being output in the diagfi.nc file. Previously it was made through a separate physiq.F
which needed to be updated every time...

--> This is completely transparent to the casual GCM users and only appears in the WRITEDIAGFI section of physiq.F

Location:

trunk

Files:

• LMDZ.MARS/README (modified) (1 diff)
• LMDZ.MARS/libf/phymars/physiq.F (modified) (3 diffs)
• MESOSCALE/LMDZ.MARS.new/libf (added)

trunk/LMDZ.MARS/README

@@ -1006 +1006 @@
  ^-----------------> HFX and USTM computations now uses the different betas for heat and momentum.
 
+== 21/09/11 == AS
+--> Added MESOINI precompiler flag so that all fields needed to initialize the mesoscale model
+     are being output in the diagfi.nc file. Previously it was made through a separate physiq.F
+     which needed to be updated every time...
+--> This is completely transparent to the casual GCM users and only appears in the WRITEDIAGFI section of physiq.F

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -1480 +1480 @@
 c     &                         '',2,tauTES)
 
+#ifdef MESOINI
+        call WRITEDIAGFI(ngrid,"emis","Surface emissivity","w.m-1",2,
+     &                  emis)
+        call WRITEDIAGFI(ngrid,"rho","density","none",3,rho)
+        call WRITEDIAGFI(ngrid,"tsoil","Soil temperature",
+     &                       "K",3,tsoil)
+        call WRITEDIAGFI(ngrid,"inertiedat","Soil inertia",
+     &                       "K",3,inertiedat)
+#endif
+
 
 c        ----------------------------------------------------------
@@ -1508 +1518 @@
          if (tracer) then
            if (water) then
+
+#ifdef MESOINI
+            !!!! waterice = q01, voir readmeteo.F90
+            call WRITEDIAGFI(ngridmx,'q01',noms(igcm_h2o_ice),
+     &                      'kg/kg',3,
+     &                       zq(1:ngridmx,1:nlayermx,igcm_h2o_ice))
+            !!!! watervapor = q02, voir readmeteo.F90
+            call WRITEDIAGFI(ngridmx,'q02',noms(igcm_h2o_vap),
+     &                      'kg/kg',3,
+     &                       zq(1:ngridmx,1:nlayermx,igcm_h2o_vap))
+            !!!! surface waterice qsurf02 (voir readmeteo)
+            call WRITEDIAGFI(ngridmx,'qsurf02','surface tracer',
+     &                      'kg.m-2',2,
+     &                       qsurf(1:ngridmx,igcm_h2o_ice))
+#endif
 
              CALL WRITEDIAGFI(ngridmx,'mtot',
@@ -1588 +1613 @@
 c            call WRITEDIAGFI(ngridmx,'dustN','Dust number',
 c    &                        'part/kg',3,pq(1,1,igcm_dust_number))
+#ifdef MESOINI
+             call WRITEDIAGFI(ngridmx,'dustN','Dust number',
+     &                        'part/kg',3,pq(1,1,igcm_dust_number))
+#endif
            else
              do iq=1,dustbin