Changeset 899 for LMDZ4/branches
- Timestamp:
- Feb 7, 2008, 1:12:46 PM (17 years ago)
- 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 75 75 ENDDO 76 76 IF (srmu .NE. 0.0) THEN 77 albedo(i) = salb / srmu * fmagic 77 albedo(i) = salb / srmu * fmagic+pmagic 78 78 ELSE ! nuit polaire (on peut prendre une valeur quelconque) 79 79 albedo(i) = fmagic … … 120 120 ENDDO 121 121 IF (srmu .NE. 0.0) THEN 122 albedo(i) = salb / srmu * fmagic 122 albedo(i) = salb / srmu * fmagic+pmagic 123 123 ELSE ! nuit polaire (on peut prendre une valeur quelconque) 124 124 albedo(i) = fmagic … … 170 170 c 171 171 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 173 173 albedo(i) = MAX(MIN(albedo(i),0.60),0.04) 174 174 ENDDO … … 181 181 rmu0(i) = MAX(rmu0(i),0.0) 182 182 cIM: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 184 184 albedo(i) = MAX(MIN(albedo(i),0.60),0.04) 185 185 ENDDO -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/clesphys.h
r857 r899 31 31 !IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH 32 32 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 35 36 !IM lev_histhf : niveau sorties 6h 36 37 !IM lev_histday : niveau sorties journalieres … … 53 54 & , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt & 54 55 & , 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 & 56 58 & , type_run, ok_isccp, ok_regdyn & 57 59 & , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins & -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/conf_phys.F90
r857 r899 83 83 REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp 84 84 LOGICAL,SAVE :: ok_kzmin_omp 85 REAL, SAVE :: fmagic_omp 85 REAL, SAVE :: fmagic_omp, pmagic_omp 86 86 INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp 87 87 CHARACTER*4, SAVE :: type_run_omp … … 603 603 ! 604 604 !Config Key = fmagic 605 !Config Desc = 605 !Config Desc = additionnal multiplicator factor used for albedo 606 606 !Config Def = 1. 607 !Config Help = Used in albedo.F607 !Config Help = additionnal multiplicator factor used in albedo.F 608 608 ! 609 609 fmagic_omp = 1. 610 610 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 611 621 612 622 !Config Key = ok_lic_melt … … 822 832 ok_kzmin = ok_kzmin_omp 823 833 fmagic = fmagic_omp 834 pmagic = pmagic_omp 824 835 iflag_pbl = iflag_pbl_omp 825 836 lev_histhf = lev_histhf_omp … … 915 926 write(numout,*)' ok_kzmin = ',ok_kzmin 916 927 write(numout,*)' fmagic = ',fmagic 928 write(numout,*)' pmagic = ',pmagic 917 929 write(numout,*)' ok_ade = ',ok_ade 918 930 write(numout,*)' ok_aie = ',ok_aie
Note: See TracChangeset
for help on using the changeset viewer.