Changeset 2428 for LMDZ5/trunk/libf/phylmd/cosp/zeff.F90

Timestamp:

Jan 27, 2016, 10:42:32 AM (8 years ago)

Author:

idelkadi

Message:

Mise a jour du simulateur COSP (passage de la version v3.2 a la version v1.4) :

• mise a jour des sources pour ISCCP, CALIPSO et PARASOL
• prise en compte des changements de phases pour les nuages (Calipso)
• rajout de plusieurs diagnostiques (fraction nuageuse en fonction de la temperature, ...)

http://lmdz.lmd.jussieu.fr/Members/aidelkadi/cosp

File:

• LMDZ5/trunk/libf/phylmd/cosp/zeff.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/cosp/zeff.F90

@@ +1 @@
+! $Revision: 88 $, $Date: 2013-11-13 15:08:38 +0100 (mer. 13 nov. 2013) $
+! $URL: http://cfmip-obs-sim.googlecode.com/svn/stable/v1.4.0/quickbeam/zeff.f90 $
   subroutine zeff(freq,D,N,nsizes,k2,tt,ice,xr,z_eff,z_ray,kr,qe,qs,rho_e)
   use math_lib
@@ -15 +17 @@
 !   [nsizes]    number of discrete drop sizes
 !   [k2]        |K|^2, -1=use frequency dependent default 
-!   [tt]        hydrometeor temperature (C)
+!   [tt]        hydrometeor temperature (K)
 !   [ice]       indicates volume consists of ice
 !   [xr]        perform Rayleigh calculations?
@@ -42 +44 @@
 ! ----- INTERNAL -----
   integer :: &
-  correct_for_rho               ! correct for density flag
+  correct_for_rho        ! correct for density flag
   real*8, dimension(nsizes) :: &
-  D0, &                         ! D in (m)
-  N0, &                         ! N in m^-3 m^-1
-  sizep, &                      ! size parameter
-  qext, &                       ! extinction efficiency
-  qbsca, &                      ! backscatter efficiency
-  rho_ice, &                    ! bulk density ice (kg m^-3)
-  f                             ! ice fraction
+  D0, &                  ! D in (m)
+  N0, &                  ! N in m^-3 m^-1
+  sizep, &               ! size parameter
+  qext, &           ! extinction efficiency
+  qbsca, &               ! backscatter efficiency
+  rho_ice, &             ! bulk density ice (kg m^-3)
+  f                 ! ice fraction
+  real*8, dimension(nsizes) :: xtemp
   real*8 :: &
-  wl, &                         ! wavelength (m)
+  wl, &                  ! wavelength (m)
   cr                            ! kr(dB/km) = cr * kr(1/km)
   complex*16 :: &
-  m                             ! complex index of refraction of bulk form
+  m                 ! complex index of refraction of bulk form
   complex*16, dimension(nsizes) :: &
-  m0                            ! complex index of refraction
+  m0                ! complex index of refraction
   
   integer*4 :: i,one
   real*8 :: pi
   real*8 :: eta_sum, eta_mie, const, z0_eff, z0_ray, k_sum, &
-            n_r, n_i, dqv(1), dqxt, dqsc, dbsc, dg, dph(1)
+            n_r, n_i, dqv(1), dqsc, dg, dph(1)
   integer*4 :: err
   complex*16 :: Xs1(1), Xs2(1)
@@ -72 +75 @@
 
 ! // conversions
-  D0 = d*1E-6                   ! m
-  N0 = n*1E12                   ! 1/(m^3 m)
-  wl = 2.99792458/(freq*10)     ! m
+  D0 = d*1E-6            ! m
+  N0 = n*1E12            ! 1/(m^3 m)
+  wl = 2.99792458/(freq*10)   ! m
   
 ! // dielectric constant |k^2| defaults
@@ -127 +130 @@
     eta_sum = qbsca(1)*(n(1)*1E6)*D0(1)**2
   else
-    call avint(qbsca*N0*D0**2,D0,nsizes,D0(1),D0(size(D0,1)),eta_sum)
+    xtemp = qbsca*N0*D0**2
+    call avint(xtemp,D0,nsizes,D0(1),D0(size(D0,1)),eta_sum)
   endif
  
@@ -140 +144 @@
     k_sum = qext(1)*(n(1)*1E6)*D0(1)**2
   else
-    call avint(qext*N0*D0**2,D0,nsizes,D0(1),D0(size(D0,1)),k_sum)
+    xtemp = qext*N0*D0**2
+    call avint(xtemp,D0,nsizes,D0(1),D0(size(D0,1)),k_sum)
   endif
   cr = 10./log(10.)
   kr = k_sum*0.25*pi*(1000.*cr)
-        
+     
 ! // z_ray = sum[D^6*N(D)*deltaD]
   if (xr == 1) then
@@ -151 +156 @@
       z0_ray = (n(1)*1E6)*D0(1)**6
     else
-      call avint(N0*D0**6,D0,nsizes,D0(1),D0(size(D0)),z0_ray)
+      xtemp = N0*D0**6
+      call avint(xtemp,D0,nsizes,D0(1),D0(size(D0)),z0_ray)
     endif
   endif