Changeset 1942 for LMDZ5/branches/LMDZ5-DOFOCO

Timestamp:

Jan 21, 2014, 4:06:36 PM (11 years ago)

Author:

jghattas

Message:

Passing the cosine of solar zenith angle from pbl_surface into orchidee.
Matthew McGrath?(LSCE)

Location:

LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd

Files:

• pbl_surface_mod.F90 (modified) (4 diffs)
• surf_land_mod.F90 (modified) (3 diffs)
• surf_land_orchidee_mod.F90 (modified) (5 diffs)

LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/pbl_surface_mod.F90

@@ -443 +443 @@
     REAL, DIMENSION(klon,klev+1)       :: ytke
     REAL, DIMENSION(klon,nsoilmx)      :: ytsoil
+    REAL, DIMENSION(klon)              :: yrmu0
     CHARACTER(len=80)                  :: abort_message
     CHARACTER(len=20)                  :: modname = 'pbl_surface'
@@ -554 +555 @@
     d_t_diss= 0.0 ;d_u = 0.0     ; d_v = 0.0        ; yqsol = 0.0    
     ytherm = 0.0  ; ytke=0.
+    yrmu0 = 0.0
     
     tke(:,:,is_ave)=0.
@@ -710 +712 @@
           ypaprs(j,klev+1) = paprs(i,klev+1)
           ywindsp(j) = SQRT(u10m(i,nsrf)**2 + v10m(i,nsrf)**2 )
+          yrmu0(j)=rmu0(i)
        END DO
 
@@ -856 +859 @@
                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
                yqsurf, ytsurf_new, y_dflux_t, y_dflux_q, &
-               y_flux_u1, y_flux_v1 )
+               y_flux_u1, y_flux_v1, yrmu0 )
                
      

LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_mod.F90

@@ -19 +19 @@
        z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
        qsurf, tsurf_new, dflux_s, dflux_l, &
-       flux_u1, flux_v1 ) 
+       flux_u1, flux_v1, rmu0) 
 
     USE dimphy
@@ -61 +61 @@
                                                          ! corresponds to previous sollwdown
     REAL, DIMENSION(klon), INTENT(IN)       :: q2m, t2m
+    REAL, DIMENSION(klon), INTENT(IN)       :: rmu0      ! cosine of solar zenith angle
 
 ! In/Output variables
@@ -132 +133 @@
             evap, fluxsens, fluxlat, &              
             tsol_rad, tsurf_new, alb1_new, alb2_new, &
-            emis_new, z0_new, qsurf)       
+            emis_new, z0_new, qsurf, rmu0)       
 
 !  

LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_orchidee_mod.F90

@@ -39 +39 @@
        evap, fluxsens, fluxlat, &              
        tsol_rad, tsurf_new, alb1_new, alb2_new, &
-       emis_new, z0_new, qsurf)
+       emis_new, z0_new, qsurf, rmu0)
 
     USE mod_surf_para
@@ -124 +124 @@
     REAL, DIMENSION(klon), INTENT(IN)         :: lwdown, swnet, swdown, ps
     REAL, DIMENSION(klon), INTENT(IN)         :: q2m, t2m
+    REAL, DIMENSION(klon), INTENT(IN)         :: rmu0
 
 ! Parametres de sortie
@@ -198 +199 @@
     INTEGER :: orch_omp_rank
     INTEGER :: orch_omp_size
+
+    REAL, DIMENSION(knon)                     :: sinang
+
 !
 ! Fin definition
@@ -402 +406 @@
                evap, fluxsens, fluxlat, coastalflow, riverflow, &
                tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
-               lon_scat, lat_scat, q2m, t2m)
+               lon_scat, lat_scat, q2m, t2m, sinang=rmu0)
 #endif         
        ENDIF
@@ -426 +430 @@
             evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
             tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), &
-            lon_scat, lat_scat, q2m, t2m)
+            lon_scat, lat_scat, q2m, t2m, sinang=rmu0(1:knon))
 #endif       
     ENDIF