Changeset 1264

Timestamp:

May 16, 2014, 11:49:52 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS. Two modifications related to outputs. 1. corrected a bug for thermospheric outputs (writediagfi under a callstats flag) 2. (only for mesoscale) made adaptations to allow for freedust runs with MESOSCALE.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• conf_phys.F (modified) (2 diffs)
• physiq.F (modified) (6 diffs)

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

@@ -370 +370 @@
            stop
          endif
+#ifndef MESOSCALE
+         ! this test is valid in GCM case
+         ! ... not in mesoscale case, for we want to activate mesoscale lifting
          if (freedust.and.lifting) then
            print*,'if freedust is used, then lifting should not be used'
@@ -375 +378 @@
            lifting=.false.
          endif
+#endif
 
 ! callddevil

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

@@ -1857 +1857 @@
 !                   call writediagfi(ngrid,"rho_"//trim(noms(iq)),
 !     $                     "Number density","cm-3",3,rhopart)
-                   if ((noms(iq).eq."o") .or. (noms(iq).eq."co2").or.
-     $                 (noms(iq).eq."o3")) then                     
-                      call writediagfi(ngrid,"vmr_"//trim(noms(iq)), 
-     $                         "Volume mixing ratio","mol/mol",3,vmr)
-                   end if
                    do ig = 1,ngrid
                       colden(ig,iq) = 0.                           
@@ -1874 +1869 @@
                    call wstats(ngrid,"c_"//trim(noms(iq)),           
      $                         "column","mol cm-2",2,colden(1,iq))  
-                   call writediagfi(ngrid,"c_"//trim(noms(iq)),  
-     $                             "column","mol cm-2",2,colden(1,iq))
                   end if ! of if (noms(iq) .ne. "dust_mass" ...)
                 end do ! of do iq=1,nq
@@ -2134 +2127 @@
 
          if (tracer.and.(dustbin.ne.0)) then
+
           call WRITEDIAGFI(ngrid,'tau','taudust','SI',2,tau(1,1))
+
+#ifndef MESOINI
            if (doubleq) then
 c            call WRITEDIAGFI(ngrid,'qsurf','qsurf',
@@ -2164 +2160 @@
 c             call WRITEDIAGFI(ngrid,"tauscaling",
 c     &                    "dust conversion factor"," ",2,tauscaling)
-#ifdef MESOINI
+           else
+             do iq=1,dustbin
+               write(str2(1:2),'(i2.2)') iq
+               call WRITEDIAGFI(ngrid,'q'//str2,'mix. ratio',
+     &                         'kg/kg',3,zq(1,1,iq))
+               call WRITEDIAGFI(ngrid,'qsurf'//str2,'qsurf',
+     &                         'kg.m-2',2,qsurf(1,iq))
+             end do
+           endif ! (doubleq)
+
+           if (scavenging) then
+             call WRITEDIAGFI(ngrid,'ccnq','CCN mass mr',
+     &                        'kg/kg',3,qccn)
+             call WRITEDIAGFI(ngrid,'ccnN','CCN number',
+     &                        'part/kg',3,nccn)
+           endif ! (scavenging)
+
+c          if (submicron) then
+c            call WRITEDIAGFI(ngrid,'dustsubm','subm mass mr',
+c    &                        'kg/kg',3,pq(1,1,igcm_dust_submicron))
+c          endif ! (submicron)
+
+#else
 !     !!! to initialize mesoscale we need scaled variables
 !     !!! because this must correspond to starting point for tracers
@@ -2175 +2193 @@
 !             call WRITEDIAGFI(ngrid,'ccnN','Nuclei number',
 !     &           'part/kg',3,pq(1:ngrid,1:nlayer,igcm_ccn_number))
-              call WRITEDIAGFI(ngrid,'dustq','Dust mass mr',
-     &                        'kg/kg',3,pq(1,1,igcm_dust_mass))
-              call WRITEDIAGFI(ngrid,'dustN','Dust number',
-     &                        'part/kg',3,pq(1,1,igcm_dust_number))
-              call WRITEDIAGFI(ngrid,'ccn','Nuclei mass mr',
-     &                        'kg/kg',3,pq(1,1,igcm_ccn_mass))
-              call WRITEDIAGFI(ngrid,'ccnN','Nuclei number',
-     &                        'part/kg',3,pq(1,1,igcm_ccn_number))
-#endif
-           else
-             do iq=1,dustbin
-               write(str2(1:2),'(i2.2)') iq
-               call WRITEDIAGFI(ngrid,'q'//str2,'mix. ratio',
-     &                         'kg/kg',3,zq(1,1,iq))
-               call WRITEDIAGFI(ngrid,'qsurf'//str2,'qsurf',
-     &                         'kg.m-2',2,qsurf(1,iq))
-             end do
-           endif ! (doubleq)
-
-           if (scavenging) then
-             call WRITEDIAGFI(ngrid,'ccnq','CCN mass mr',
+      if (freedust) then
+             call WRITEDIAGFI(ngrid,'dustq','Dust mass mr',
+     &                        'kg/kg',3,qdust)
+             call WRITEDIAGFI(ngrid,'dustN','Dust number',
+     &                        'part/kg',3,ndust)
+             call WRITEDIAGFI(ngrid,'ccn','CCN mass mr',
      &                        'kg/kg',3,qccn)
              call WRITEDIAGFI(ngrid,'ccnN','CCN number',
      &                        'part/kg',3,nccn)
-           endif ! (scavenging)
-
-c          if (submicron) then
-c            call WRITEDIAGFI(ngrid,'dustsubm','subm mass mr',
-c    &                        'kg/kg',3,pq(1,1,igcm_dust_submicron))
-c          endif ! (submicron)
+      else
+             call WRITEDIAGFI(ngrid,'dustq','Dust mass mr',
+     &                        'kg/kg',3,pq(1,1,igcm_dust_mass))
+             call WRITEDIAGFI(ngrid,'dustN','Dust number',
+     &                        'part/kg',3,pq(1,1,igcm_dust_number))
+             call WRITEDIAGFI(ngrid,'ccn','Nuclei mass mr',
+     &                        'kg/kg',3,pq(1,1,igcm_ccn_mass))
+             call WRITEDIAGFI(ngrid,'ccnN','Nuclei number',
+     &                        'part/kg',3,pq(1,1,igcm_ccn_number))
+      endif
+#endif
+
          end if  ! (tracer.and.(dustbin.ne.0))
 
@@ -2222 +2231 @@
             call WRITEDIAGFI(ngrid,"qnir","NIR heating","K/s",
      $           3,zdtnirco2)
+
+
+             if (thermochem.or.photochem) then
+                do iq=1,nq
+                  if (noms(iq) .ne. "dust_mass" .and.
+     $                 noms(iq) .ne. "dust_number" .and.
+     $                 noms(iq) .ne. "ccn_mass" .and.
+     $                 noms(iq) .ne. "ccn_number" .and.
+     $                 noms(iq) .ne. "h2o_vap" .and.
+     $                 noms(iq) .ne. "h2o_ice") then
+                    vmr(1:ngrid,1:nlayer)=zq(1:ngrid,1:nlayer,iq)
+     &                          *mmean(1:ngrid,1:nlayer)/mmol(iq)
+                    rhopart(1:ngrid,1:nlayer)=zq(1:ngrid,1:nlayer,iq)
+     &                          *rho(1:ngrid,1:nlayer)*n_avog/
+     &                           (1000*mmol(iq))
+                   call writediagfi(ngrid,"rho_"//trim(noms(iq)),
+     $                     "Number density","cm-3",3,rhopart)
+                   if ((noms(iq).eq."o") .or. (noms(iq).eq."co2").or.
+     $                 (noms(iq).eq."o3")) then
+                      call writediagfi(ngrid,"vmr_"//trim(noms(iq)),
+     $                         "Volume mixing ratio","mol/mol",3,vmr)
+                   end if
+                   do ig = 1,ngrid
+                      colden(ig,iq) = 0.
+                   end do
+                   do l=1,nlayer
+                      do ig=1,ngrid
+                         colden(ig,iq) = colden(ig,iq) + zq(ig,l,iq)
+     $                                  *(zplev(ig,l)-zplev(ig,l+1))
+     $                                  *6.022e22/(mmol(iq)*g)
+                      end do
+                   end do
+                   call writediagfi(ngrid,"c_"//trim(noms(iq)),
+     $                             "column","mol cm-2",2,colden(1,iq))
+                  end if ! of if (noms(iq) .ne. "dust_mass" ...)
+                end do ! of do iq=1,nq
+             end if ! of if (thermochem.or.photochem)
+
 
          endif  !(callthermos)