Ignore:
Timestamp:
Dec 18, 2015, 8:27:39 PM (9 years ago)
Author:
oboucher
Message:

limiting iflag_rrtm=0 to NSW=2 case (NSW=4 and 6 were possible with slightly different averaging procedures)
limiting iflag_rrtm choices to 0 and 1
limiting iflag_albedo choices to 0 and 1
cleaning up surface ocean albedo parametrisations
fmagic and pmagic corrections are now done in surf_ocean_albedo.F90
fmagic and pmagic corrections also possible for iflag_rrtm=1 now

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/albedo.F90

    r2346 r2413  
    2222    include "clesphys.h"
    2323
    24     ! fmagic -> clesphys.h/.inc
    25     ! REAL fmagic ! un facteur magique pour regler l'albedo
    26     ! cc      PARAMETER (fmagic=0.7)
    27     ! ccIM => a remplacer
    28     ! PARAMETER (fmagic=1.32)
    29     ! PARAMETER (fmagic=1.0)
    30     ! PARAMETER (fmagic=0.7)
    3124    INTEGER npts ! il controle la precision de l'integration
    3225    PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes
     
    7467          END DO
    7568          IF (srmu/=0.0) THEN
    76              albedo(i) = salb/srmu*fmagic + pmagic
     69             albedo(i) = salb/srmu
    7770          ELSE ! nuit polaire (on peut prendre une valeur quelconque)
    78              albedo(i) = fmagic
     71             albedo(i) = 1.0
    7972          END IF
    8073       END DO
     
    119112          END DO
    120113          IF (srmu/=0.0) THEN
    121              albedo(i) = salb/srmu*fmagic + pmagic
     114             albedo(i) = salb/srmu
    122115          ELSE ! nuit polaire (on peut prendre une valeur quelconque)
    123              albedo(i) = fmagic
     116             albedo(i) = 1.0
    124117          END IF
    125118       END DO
     
    146139    real, intent(out):: albedo(klon)
    147140
    148     ! REAL fmagic ! un facteur magique pour regler l'albedo
    149     ! cc      PARAMETER (fmagic=0.7)
    150     ! ccIM => a remplacer
    151     ! PARAMETER (fmagic=1.32)
    152     ! PARAMETER (fmagic=1.0)
    153     ! PARAMETER (fmagic=0.7)
    154 
    155141    REAL fauxo
    156142    INTEGER i
     
    161147       DO i = 1, klon
    162148          fauxo = (1.47-acos(max(rmu0(i), 0.0)))/0.15
    163           albedo(i) = fmagic*(.03+.630/(1.+fauxo*fauxo)) + pmagic
     149          albedo(i) = 0.03+.630/(1.+fauxo*fauxo)
    164150          albedo(i) = max(min(albedo(i),0.60), 0.04)
    165151       END DO
    166152    ELSE
    167153       DO i = 1, klon
    168           albedo(i) = fmagic*0.058/(max(rmu0(i), 0.0)+0.30) + pmagic
     154          albedo(i) = 0.058/(max(rmu0(i), 0.0)+0.30)
    169155          albedo(i) = max(min(albedo(i),0.60), 0.04)
    170156       END DO
Note: See TracChangeset for help on using the changeset viewer.