Changeset 888 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Feb 4, 2008, 5:24:28 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Modifications sur l'albedo JG
LF

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (13 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -690 +690 @@
       REAL,save ::  solarlong0
 c
-      REAL,allocatable,save :: falbe(:,:)
-c$OMP THREADPRIVATE(falbe)       ! albedo par type de surface
-c
-      REAL,allocatable,save :: falblw(:,:)
-c$OMP THREADPRIVATE(falblw)       ! albedo par type de surface
+      REAL,allocatable,save :: falb1(:,:)
+c$OMP THREADPRIVATE(falb1)       ! albedo par type de surface pour SW visible 
+c
+      REAL,allocatable,save :: falb2(:,:)
+c$OMP THREADPRIVATE(falb2)       ! albedo par type de surface pour SW proche IR
 
 c
@@ -969 +969 @@
 cIM
 cym      SAVE pctsrf                 ! sous-fraction du sol
-      REAL,allocatable,save :: albsol(:)
-c$OMP THREADPRIVATE(albsol)
-cym      SAVE albsol                 ! albedo du sol total
-      REAL,allocatable,save :: albsollw(:)
-c$OMP THREADPRIVATE(albsollw)      
-cym      SAVE albsollw                 ! albedo du sol total
+
+      REAL,allocatable,save :: albsol1(:) ! albedo du sol total pour SW visible
+c$OMP THREADPRIVATE(albsol1)
+      REAL,allocatable,save :: albsol2(:) ! albedo du sol total pour SW proche IR
+c$OMP THREADPRIVATE(albsol2)      
 
       REAL,allocatable,save :: wo(:,:)
@@ -1527 +1526 @@
       allocate( ftsol(klon,nbsrf))
       allocate( deltat(klon))
-      allocate( falbe(klon,nbsrf))
-      allocate( falblw(klon,nbsrf))
+      allocate( falb1(klon,nbsrf))
+      allocate( falb2(klon,nbsrf))
       allocate( zmea(klon))
       allocate( zstd(klon))
@@ -1573 +1572 @@
       allocate( total_rain(klon), nday_rain(klon))
       allocate( pctsrf(klon,nbsrf))
-      allocate( albsol(klon))
-      allocate( albsollw(klon))
+      allocate( albsol1(klon))
+      allocate( albsol2(klon))
       allocate( wo(klon,klev))
       allocate( clwcon(klon,klev),rnebcon(klon,klev))
@@ -1731 +1730 @@
 c        nhistoW(:,:,:,:) = 0.0
 c        histoW(:,:,:,:) = 0.0
+! fin anne
 
 cIM      
@@ -1773 +1773 @@
      .       rlat,rlon,pctsrf, ftsol,
      .       ocean, ok_veget,
-     .       falbe, falblw, rain_fall,snow_fall,
+     .       falb1, falb2, rain_fall,snow_fall,
      .       solsw, sollw,
      .       radsol,clesphy0,
@@ -2318 +2318 @@
      e     t_seri,    q_seri,    u_seri,  v_seri,   
      e     pplay,     paprs,     pctsrf,            
-     +     ftsol,     falbe,     falblw,  u10m,   v10m,
+     +     ftsol,     falb1,     falb2,   u10m,   v10m,
      s     sollwdown, cdragh,    cdragm,  yu1,    yv1,
-     s     albsol,    albsollw,  sens,    evap,  
+     s     albsol1,   albsol2,   sens,    evap,  
      s     zxtsol,    zxfluxlat, zt2m,    qsat2m, 
      s     d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf,
@@ -3233 +3233 @@
      $                          zxsnow_dummy,
      $                          solsw,
-     $                          albsol,
+     $                          albsol1,
      $                          rain_fall,
      $                          snow_fall,
@@ -3294 +3294 @@
 
       DO i = 1, klon
-         albsol(i) = falbe(i,is_oce) * pctsrf(i,is_oce)
-     .             + falbe(i,is_lic) * pctsrf(i,is_lic)
-     .             + falbe(i,is_ter) * pctsrf(i,is_ter)
-     .             + falbe(i,is_sic) * pctsrf(i,is_sic)
-         albsollw(i) = falblw(i,is_oce) * pctsrf(i,is_oce)
-     .               + falblw(i,is_lic) * pctsrf(i,is_lic)
-     .               + falblw(i,is_ter) * pctsrf(i,is_ter)
-     .               + falblw(i,is_sic) * pctsrf(i,is_sic)
+         albsol1(i) = falb1(i,is_oce) * pctsrf(i,is_oce)
+     .             + falb1(i,is_lic) * pctsrf(i,is_lic)
+     .             + falb1(i,is_ter) * pctsrf(i,is_ter)
+     .             + falb1(i,is_sic) * pctsrf(i,is_sic)
+         albsol2(i) = falb2(i,is_oce) * pctsrf(i,is_oce)
+     .               + falb2(i,is_lic) * pctsrf(i,is_lic)
+     .               + falb2(i,is_ter) * pctsrf(i,is_ter)
+     .               + falb2(i,is_sic) * pctsrf(i,is_sic)
       ENDDO
 
@@ -3335 +3335 @@
       CALL radlwsw ! nouveau rayonnement (compatible Arpege-IFS)
      e            (dist, rmu0, fract, 
-     e             paprs, pplay,zxtsol,albsol, albsollw, t_seri,q_seri,
+     e             paprs, pplay,zxtsol,albsol1, albsol2, t_seri,q_seri,
      e             wo,
      e             cldfra, cldemi, cldtau,
@@ -3574 +3574 @@
      I                   pphis,
      I                   pphi,
-     I                   albsol,
+     I                   albsol1,
      I                   qx(1,1,1), 
      I                   rhcl,
@@ -3827 +3827 @@
 c pour ecriture dans histxxx.nc 
 
-      ! Get some variables from module mod_fonte_neige
+      ! Get some variables from module fonte_neige_mod
       CALL fonte_neige_get_vars(pctsrf, 
      .     zxfqcalving, zxfqfonte, zxffonte)
 
       IF (ocean == 'slab') THEN
-         ! Get some variables from module oceanslab
+         ! Get some variables from module ocean_slab_mod
          CALL ocean_slab_get_vars(tslab, seaice, fluxo, fluxg)
       ELSEIF (ocean == 'couple') THEN
-         ! Get some variables from module oceancpl
+         ! Get some variables from module ocean_cpl_mod
          CALL ocean_cpl_get_vars(fluxo, fluxg)
       ELSE 
-         ! Get some variables from module oceanforced
+         ! Get some variables from module ocean_forced_mod
          CALL ocean_forced_get_vars(fluxo, fluxg)
       ENDIF
@@ -3884 +3884 @@
          CALL phyredem ("restartphy.nc",dtime,radpas,ocean,
      .      rlat, rlon, pctsrf, ftsol,
-     .      falbe,falblw, rain_fall, 
+     .      falb1, falb2, rain_fall, 
      .      snow_fall,
      .      solsw, sollw,