Changeset 900 for LMDZ4/trunk


Ignore:
Timestamp:
Feb 7, 2008, 3:26:01 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Inclusion de la constante pmagic dans le calcul de l'albedo MAFO/JLD
LF

Location:
LMDZ4/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/albedo.F

    r830 r900  
    7575      ENDDO
    7676      IF (srmu .NE. 0.0) THEN
    77          albedo(i) = salb / srmu * fmagic
     77         albedo(i) = salb / srmu * fmagic+pmagic
    7878      ELSE ! nuit polaire (on peut prendre une valeur quelconque)
    7979         albedo(i) = fmagic
     
    120120      ENDDO
    121121      IF (srmu .NE. 0.0) THEN
    122          albedo(i) = salb / srmu * fmagic
     122         albedo(i) = salb / srmu * fmagic+pmagic
    123123      ELSE ! nuit polaire (on peut prendre une valeur quelconque)
    124124         albedo(i) = fmagic
     
    170170c
    171171         fauxo = ( 1.47 - ACOS( rmu0(i) ) )/0.15
    172          albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))
     172         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))+pmagic
    173173         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
    174174      ENDDO
     
    181181         rmu0(i) = MAX(rmu0(i),0.0)
    182182cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)
    183          albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)
     183         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)+pmagic
    184184         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
    185185      ENDDO
  • LMDZ4/trunk/libf/phylmd/clesphys.h

    r889 r900  
    3131!IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH
    3232       LOGICAL ok_kzmin
    33 !IM fmagic : parametre pour regler l albedo sur ocean
    34        REAL fmagic
     33!IM, MAFo fmagic, pmagic : parametres - additionnel et multiplicatif -
     34!                          pour regler l albedo sur ocean
     35       REAL fmagic, pmagic
    3536! Hauteur (imposee) du contenu en eau du sol
    3637           REAL qsol0
     
    5354     &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &
    5455     &     , top_height, overlap, cdmmax, cdhmax, ksta, ksta_ter        &
    55      &     , ok_kzmin, fmagic, lev_histhf, lev_histday, lev_histmth     &
     56     &     , ok_kzmin, fmagic, pmagic                                   &
     57     &     , lev_histhf, lev_histday, lev_histmth                       &
    5658     &     , type_run, ok_isccp, ok_regdyn                              &
    5759     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
  • LMDZ4/trunk/libf/phylmd/conf_phys.F90

    r889 r900  
    101101  REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
    102102  LOGICAL,SAVE :: ok_kzmin_omp
    103   REAL, SAVE :: fmagic_omp
     103  REAL, SAVE ::  fmagic_omp, pmagic_omp
    104104  INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
    105105  CHARACTER*4, SAVE :: type_run_omp
     
    685685!
    686686!Config Key  = fmagic
    687 !Config Desc =
     687!Config Desc = additionnal multiplicator factor used for albedo
    688688!Config Def  = 1.
    689 !Config Help = Used in albedo.F
     689!Config Help = additionnal multiplicator factor used in albedo.F
    690690!
    691691  fmagic_omp = 1.
    692692  call getin('fmagic',fmagic_omp)
     693
     694!
     695!Config Key  = pmagic
     696!Config Desc = additional factor used for albedo
     697!Config Def  = 0.
     698!Config Help = additional factor used in albedo.F
     699!
     700  pmagic_omp = 0.
     701  call getin('pmagic',pmagic_omp)
     702
    693703
    694704!Config Key = ok_lic_melt
     
    941951    ok_kzmin = ok_kzmin_omp
    942952    fmagic = fmagic_omp
     953    pmagic = pmagic_omp
    943954    iflag_pbl = iflag_pbl_omp
    944955    lev_histhf = lev_histhf_omp
     
    10431054  write(numout,*)' ok_kzmin = ',ok_kzmin
    10441055  write(numout,*)' fmagic = ',fmagic
     1056  write(numout,*)' pmagic = ',pmagic
    10451057  write(numout,*)' ok_ade = ',ok_ade
    10461058  write(numout,*)' ok_aie = ',ok_aie
Note: See TracChangeset for help on using the changeset viewer.