Changeset 741 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jul 25, 2012, 3:20:21 PM (12 years ago)

Author:

jleconte

Message:

Model now works without aerosols at all. You still need to compile it with -s 1 option. LK12

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• aeropacity.F90 (modified) (2 diffs)
• aerosol_mod.F90 (modified) (1 diff)
• iniaerosol.F (modified) (2 diffs)

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

@@ -116 +116 @@
         enddo
 
-        if (naerkind.ne.0) then
+        if (noaero) then
+          print*, "No active aerosols found in aeropacity"
+        else
           print*, "If you would like to use aerosols, make sure any old"
           print*, "start files are updated in newstart using the option"
           print*, "q=0"
-          write(*,*) "Aerosols found in aeropacity:"
-        endif
-
-        if (iaero_co2.ne.0) then
+          write(*,*) "Active aerosols found in aeropacity:"
+        endif
+
+        if ((iaero_co2.ne.0).and.(.not.noaero)) then
           print*, 'iaero_co2=  ',iaero_co2
         endif
@@ -150 +152 @@
             aerosol(1:ngridmx,1:nlayermx,iaer)=0.0
 !       2. Opacity calculation
-            if (aerofixco2.or.(i_co2ice.eq.0)) then !  CO2 ice cloud prescribed
+            if (noaero) then ! aerosol set to a very low value
+             aerosol(1:ngridmx,1:nlayermx,iaer)=1.e-9
+            elseif (aerofixco2.or.(i_co2ice.eq.0)) then !  CO2 ice cloud prescribed
                aerosol(1:ngridmx,1:nlayermx,iaer)=1.e-9
                !aerosol(1:ngridmx,12,iaer)=4.0 ! single cloud layer option

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

@@ -9 +9 @@
       integer, save :: iaero_dust
       integer, save :: iaero_h2so4
-
+      logical, save :: noaero
+      
 !==================================================================
 end module aerosol_mod

trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F

@@ -28 +28 @@
       iaero_dust  =0
       iaero_h2so4 =0
-      
+      noaero = .false.
+ 
       ia=0
       if (aeroco2) then
-      ia=ia+1
-      iaero_co2=ia
-      write(*,*) 'CO2 aerosol = ', iaero_co2
+         ia=ia+1
+         iaero_co2=ia
+         write(*,*) 'CO2 aerosol = ', iaero_co2
       endif
  
       if (aeroh2o) then
-      ia=ia+1
-      iaero_h2o=ia
-      write(*,*) 'H2O aerosol = ', iaero_h2o
+         ia=ia+1
+         iaero_h2o=ia
+         write(*,*) 'H2O aerosol = ', iaero_h2o
       endif
 
       if (dusttau.gt.0) then
-      ia=ia+1
-      iaero_dust=ia
-      write(*,*) 'Dust aerosol = ', iaero_dust
+         ia=ia+1
+         iaero_dust=ia
+         write(*,*) 'Dust aerosol = ', iaero_dust
       endif
 
       if (aeroh2so4) then
-      ia=ia+1
-      iaero_h2so4=ia
-      write(*,*) 'H2SO4 aerosol = ', iaero_h2so4
+         ia=ia+1
+         iaero_h2so4=ia
+         write(*,*) 'H2SO4 aerosol = ', iaero_h2so4
       endif
       
       write(*,*) 'Number of aerosols= ', ia
       
-      if (ia.eq.0) then  !For the zero aerosol case
-      ia = 1
+! For the zero aerosol case, we currently make a dummy co2 aerosol which is zero everywhere.
+! (See aeropacity.F90 for how this works). A better solution would be to turn off the 
+! aerosol machinery in the no aerosol case, but this would be complicated. LK
+
+      if (ia.eq.0) then  !For the zero aerosol case. 
+         ia = 1
+         noaero = .true.
+         iaero_co2=ia
       endif
 
@@ -65 +72 @@
           print*, 'or change options in callphys.def'
           print*, 'Abort in iniaerosol.F'
-      call abort
+          call abort
       endif