Changeset 368 for trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

Timestamp:

Nov 9, 2011, 4:40:11 PM (13 years ago)

Author:

rwordsworth

Message:

.nc saves in physiq.F90 adapted to new diagfi.def format.
More info added to gases.def warnings.

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -2005 +2005 @@
         call writediagfi(ngrid,"tsurf","Surface temperature","K",2,tsurf)
         call writediagfi(ngrid,"ps","Surface pressure","Pa",2,ps)
-
-        if(1.eq.2)then
         call writediagfi(ngrid,"temp","temperature","K",3,zt)
         call writediagfi(ngrid,"u","Zonal wind","m.s-1",3,zu)
@@ -2012 +2010 @@
         call writediagfi(ngrid,"w","Vertical wind","m.s-1",3,pw)
         call writediagfi(ngrid,'p','Pressure','Pa',3,pplay)
-        endif
 
 !     Total energy balance diagnostics
@@ -2036 +2033 @@
         !call writediagfi(ngrid,"vdifNC","H2O loss vdifc","kg m-2 s-1",2,vdifcncons)
         !call writediagfi(ngrid,"cadjNC","H2O loss convadj","kg m-2 s-1",2,cadjncons)
-        call writediagfi(ngrid,'rnat','Terrain type',' ',2,real(rnat))
+        !call writediagfi(ngrid,'rnat','Terrain type',' ',2,real(rnat))
         !call writediagfi(ngrid,'pphi','Geopotential',' ',3,pphi)
-        !print*,vdifcncons
-        !call abort
 
 !     Output aerosols
-!        call writediagfi(ngridmx,'CO2ice_reff','CO2ice_reff','m',3,reffrad(1,1,1))
-!        call writediagfi(ngridmx,'H2Oice_reff','H2Oice_reff','m',3,reffrad(1,1,2))
-!        call writediagfi(ngridmx,'CO2ice_reffcol','CO2ice_reffcol','um kg m^-2',2,reffcol(1,1))
-!        call writediagfi(ngridmx,'H2Oice_reffcol','H2Oice_reffcol','um kg m^-2',2,reffcol(1,2))
+        call writediagfi(ngridmx,'CO2ice_reff','CO2ice_reff','m',3,reffrad(1,1,1))
+        call writediagfi(ngridmx,'H2Oice_reff','H2Oice_reff','m',3,reffrad(1,1,2))
+        call writediagfi(ngridmx,'CO2ice_reffcol','CO2ice_reffcol','um kg m^-2',2,reffcol(1,1))
+        call writediagfi(ngridmx,'H2Oice_reffcol','H2Oice_reffcol','um kg m^-2',2,reffcol(1,2))
 
 !     Output tracers
        if (tracer) then
         do iq=1,nq
-!          call writediagfi(ngrid,noms(iq),noms(iq),'kg/kg',3,zq(1,1,iq))
+          call writediagfi(ngrid,noms(iq),noms(iq),'kg/kg',3,zq(1,1,iq))
 !          call writediagfi(ngridmx,trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  & 
 !                          'kg m^-2',2,qsurf(1,iq) )