Ignore:
Timestamp:
Feb 7, 2008, 1:12:46 PM (17 years ago)
Author:
Laurent Fairhead
Message:

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

Location:
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/albedo.F

    r845 r899  
    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/branches/LMDZ4_V3_patches/libf/phylmd/clesphys.h

    r857 r899  
    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!IM lev_histhf  : niveau sorties 6h
    3637!IM lev_histday : niveau sorties journalieres
     
    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/branches/LMDZ4_V3_patches/libf/phylmd/conf_phys.F90

    r857 r899  
    8383  REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
    8484  LOGICAL,SAVE :: ok_kzmin_omp
    85   REAL, SAVE ::  fmagic_omp
     85  REAL, SAVE ::  fmagic_omp, pmagic_omp
    8686  INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
    8787  CHARACTER*4, SAVE :: type_run_omp
     
    603603!
    604604!Config Key  = fmagic
    605 !Config Desc =
     605!Config Desc = additionnal multiplicator factor used for albedo
    606606!Config Def  = 1.
    607 !Config Help = Used in albedo.F
     607!Config Help = additionnal multiplicator factor used in albedo.F
    608608!
    609609  fmagic_omp = 1.
    610610  call getin('fmagic',fmagic_omp)
     611
     612!
     613!Config Key  = pmagic
     614!Config Desc = additional factor used for albedo
     615!Config Def  = 0.
     616!Config Help = additional factor used in albedo.F
     617!
     618  pmagic_omp = 0.
     619  call getin('pmagic',pmagic_omp)
     620
    611621
    612622!Config Key = ok_lic_melt
     
    822832    ok_kzmin = ok_kzmin_omp
    823833    fmagic = fmagic_omp
     834    pmagic = pmagic_omp
    824835    iflag_pbl = iflag_pbl_omp
    825836    lev_histhf = lev_histhf_omp
     
    915926  write(numout,*)' ok_kzmin = ',ok_kzmin
    916927  write(numout,*)' fmagic = ',fmagic
     928  write(numout,*)' pmagic = ',pmagic
    917929  write(numout,*)' ok_ade = ',ok_ade
    918930  write(numout,*)' ok_aie = ',ok_aie
Note: See TracChangeset for help on using the changeset viewer.