- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_sno_albedo.f90
r5113 r5116 67 67 use VARySV 68 68 use VARtSV 69 USE surface_data, only: iflag_albcalc,correc_alb69 USE surface_data, ONLY: iflag_albcalc,correc_alb 70 70 71 71 IMPLICIT NONE … … 73 73 74 74 ! + -- INPUT 75 integer:: jjtime75 INTEGER :: jjtime 76 76 77 77 ! +--Internal Variables 78 78 ! + ================== 79 79 80 real:: coalb1(knonv) ! weighted Coalbedo, Vis.81 real:: coalb2(knonv) ! weighted Coalbedo, nIR 182 real:: coalb3(knonv) ! weighted Coalbedo, nIR 283 real:: coalbm ! weighted Coalbedo, mean84 real:: sExt_1(knonv) ! Extinction Coeff., Vis.85 real:: sExt_2(knonv) ! Extinction Coeff., nIR 186 real:: sExt_3(knonv) ! Extinction Coeff., nIR 287 real:: SnOpSV(knonv, nsno) ! Snow Grain optical Size80 REAL :: coalb1(knonv) ! weighted Coalbedo, Vis. 81 REAL :: coalb2(knonv) ! weighted Coalbedo, nIR 1 82 REAL :: coalb3(knonv) ! weighted Coalbedo, nIR 2 83 REAL :: coalbm ! weighted Coalbedo, mean 84 REAL :: sExt_1(knonv) ! Extinction Coeff., Vis. 85 REAL :: sExt_2(knonv) ! Extinction Coeff., nIR 1 86 REAL :: sExt_3(knonv) ! Extinction Coeff., nIR 2 87 REAL :: SnOpSV(knonv, nsno) ! Snow Grain optical Size 88 88 ! #AG real agesno 89 89 90 integer:: isn ,ikl ,isn1, i91 real:: sbeta1,sbeta2,sbeta3,sbeta4,sbeta592 real:: AgeMax,AlbMin,HSnoSV,HIceSV,doptmx,SignG1,Sph_OK93 real:: dalbed,dalbeS,dalbeW94 real:: bsegal,czemax,csegal,csza95 real:: RoFrez,DiffRo,SignRo,SnowOK,OpSqrt96 real:: albSn1,albIc1,a_SnI1,a_SII197 real:: albSn2,albIc2,a_SnI2,a_SII298 real:: albSn3,albIc3,a_SnI3,a_SII399 real:: albSno,albIce,albSnI,albSII,albWIc,albmax100 real:: doptic,Snow_H,SIce_H,SnownH,SIcenH101 real:: exarg1,exarg2,exarg3,sign_0,sExt_0102 real:: albedo_old,albCor103 real:: ro_ave,dz_ave,minalb104 real:: l1min,l1max,l2min,l2max,l3min,l3max105 real:: l6min(6), l6max(6), albSn6(6), a_SII6(6)106 real:: lmintmp,lmaxtmp,albtmp90 INTEGER :: isn ,ikl ,isn1, i 91 REAL :: sbeta1,sbeta2,sbeta3,sbeta4,sbeta5 92 REAL :: AgeMax,AlbMin,HSnoSV,HIceSV,doptmx,SignG1,Sph_OK 93 REAL :: dalbed,dalbeS,dalbeW 94 REAL :: bsegal,czemax,csegal,csza 95 REAL :: RoFrez,DiffRo,SignRo,SnowOK,OpSqrt 96 REAL :: albSn1,albIc1,a_SnI1,a_SII1 97 REAL :: albSn2,albIc2,a_SnI2,a_SII2 98 REAL :: albSn3,albIc3,a_SnI3,a_SII3 99 REAL :: albSno,albIce,albSnI,albSII,albWIc,albmax 100 REAL :: doptic,Snow_H,SIce_H,SnownH,SIcenH 101 REAL :: exarg1,exarg2,exarg3,sign_0,sExt_0 102 REAL :: albedo_old,albCor 103 REAL :: ro_ave,dz_ave,minalb 104 REAL :: l1min,l1max,l2min,l2max,l3min,l3max 105 REAL :: l6min(6), l6max(6), albSn6(6), a_SII6(6) 106 REAL :: lmintmp,lmaxtmp,albtmp 107 107 108 108 ! +--Local DATA … … 464 464 ! prescription for each time step with NEMO values 465 465 466 ! #AO if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. .TRUE.) then467 ! #AO if (AOmask(ikl) .eq. 0) then466 ! #AO if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. .TRUE.) THEN 467 ! #AO if (AOmask(ikl) .eq. 0) THEN 468 468 ! #AO albisv(ikl) = (1.-AOmask(ikl))* albAOsisv(ikl) 469 469 ! #AO. +(AOmask(ikl)*albisv(ikl)) … … 569 569 570 570 571 end subroutinesnoptp571 END SUBROUTINE snoptp 572 572 573 573
Note: See TracChangeset
for help on using the changeset viewer.