Changeset 1529 for trunk/LMDZ.GENERIC/libf/phystd/radii_mod.F90

Timestamp:

Apr 5, 2016, 10:51:51 AM (9 years ago)

Author:

emillour

Message:

Generic GCM: Towards a cleaner separation between physics and dynamics

• Got rid of references to "dimensions.h" from physics packages: use nbp_lon (=iim), nbp_lat (=jjp1) and nbp_lev from module mod_grid_phy_lmdz (in phy_common) instead.
• Removed module "comhdiff_mod.F90", as it is only used by module surf_heat_transp_mod.F90, moved module variables there.
• Addedin "surf_heat_transp_mod" local versions of some arrays and routines (from dyn3d) required to compute gradient, divergence, etc. on the global dynamics grid. As before, the slab ocean only works in serial.

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/radii_mod.F90 (modified) (8 diffs)

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

@@ -8 +8 @@
 !     water cloud optical properties
 
-      use callkeys_mod, only: radfixed,Nmix_co2,                        &
-                pres_bottom_tropo,pres_top_tropo,size_tropo,    &
-                pres_bottom_strato,size_strato
+      use callkeys_mod, only: radfixed,Nmix_co2,                    &
+                pres_bottom_tropo,pres_top_tropo,size_tropo,        &
+                pres_bottom_strato,size_strato
       
       real, save ::  rad_h2o
@@ -21 +21 @@
 
 
-      contains
+contains
 
 
@@ -38 +38 @@
       use ioipsl_getin_p_mod, only: getin_p
       use radinc_h, only: naerkind
-      use aerosol_mod
-!      USE tracer_h
-      Implicit none
-
-!      include "dimensions.h"
-!      include "dimphys.h"
+      use aerosol_mod, only: iaero_back2lay, iaero_co2, iaero_dust, &
+                             iaero_h2o, iaero_h2so4
+      Implicit none
 
       integer,intent(in) :: ngrid
@@ -81 +78 @@
                nueffrad(1:ngrid,1:nlayer,iaer) = 0.1 
             endif
-            
-            if(iaer.eq.iaero_back2lay)then ! Two-layer aerosols
+            
+            if(iaer.eq.iaero_back2lay)then ! Two-layer aerosols
                reffrad(1:ngrid,1:nlayer,iaer) = 2.e-6
                nueffrad(1:ngrid,1:nlayer,iaer) = 0.1 
@@ -101 +98 @@
          if (radfixed) then
 
-            write(*,*)"radius of H2O water particles:"
+            write(*,*)"radius of H2O water particles:"
             rad_h2o=13. ! default value
             call getin_p("rad_h2o",rad_h2o)
             write(*,*)" rad_h2o = ",rad_h2o
 
-            write(*,*)"radius of H2O ice particles:"
+            write(*,*)"radius of H2O ice particles:"
             rad_h2o_ice=35. ! default value
             call getin_p("rad_h2o_ice",rad_h2o_ice)
             write(*,*)" rad_h2o_ice = ",rad_h2o_ice
 
-         else
+         else
 
             write(*,*)"Number mixing ratio of H2O water particles:"
@@ -122 +119 @@
             call getin_p("Nmix_h2o_ice",Nmix_h2o_ice)
             write(*,*)" Nmix_h2o_ice = ",Nmix_h2o_ice
-         endif
+         endif
 
       print*,'exit su_aer_radii'
@@ -173 +170 @@
                zfice = 1.0 - (pt(ig,l)-T_h2O_ice_clouds) / (T_h2O_ice_liq-T_h2O_ice_clouds)
                zfice = MIN(MAX(zfice,0.0),1.0)
-               zrad_liq  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o*pi*rhowater) )
-               zrad_ice  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o_ice*pi*rhowaterice) )
+               zrad_liq  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o*pi*rhowater) )
+               zrad_ice  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o_ice*pi*rhowaterice) )
                nueffrad(ig,l) = coef_chaud * (1.-zfice) + coef_froid * zfice
                zrad = zrad_liq * (1.-zfice) + zrad_ice * zfice
 
-               reffrad(ig,l) = min(max(zrad,1.e-6),1000.e-6)
+               reffrad(ig,l) = min(max(zrad,1.e-6),1000.e-6)
                enddo
             enddo      
@@ -213 +210 @@
 
       if (radfixed) then
-         reffliq(1:ngrid,1:nlayer)= rad_h2o
-         reffice(1:ngrid,1:nlayer)= rad_h2o_ice
+         reffliq(1:ngrid,1:nlayer)= rad_h2o
+         reffice(1:ngrid,1:nlayer)= rad_h2o_ice
       else
-         do k=1,nlayer
-           do i=1,ngrid
-             reffliq(i,k) = CBRT(3*pql(i,k)/(4*Nmix_h2o*pi*rhowater))
-             reffliq(i,k) = min(max(reffliq(i,k),1.e-6),1000.e-6)
-           
-             reffice(i,k) = CBRT(3*pql(i,k)/(4*Nmix_h2o_ice*pi*rhowaterice))
-             reffice(i,k) = min(max(reffice(i,k),1.e-6),1000.e-6)
-           enddo
-         enddo
+         do k=1,nlayer
+           do i=1,ngrid
+             reffliq(i,k) = CBRT(3*pql(i,k)/(4*Nmix_h2o*pi*rhowater))
+             reffliq(i,k) = min(max(reffliq(i,k),1.e-6),1000.e-6)
+           
+             reffice(i,k) = CBRT(3*pql(i,k)/(4*Nmix_h2o_ice*pi*rhowaterice))
+             reffice(i,k) = min(max(reffice(i,k),1.e-6),1000.e-6)
+           enddo
+         enddo
       endif