Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Timestamp:
- Jul 23, 2024, 3:29:36 PM (13 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Files:
-
- 14 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARdSV.F90
r5099 r5103 9 9 ! + ======================= 10 10 11 LOGICAL :: INI_SV=. false. ! Initialisation Switch11 LOGICAL :: INI_SV=.FALSE. ! Initialisation Switch 12 12 REAL,PARAMETER :: eps_21=1.e-21 ! Arbitrary very small value 13 13 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/inlandsis.F
r5099 r5103 1 subroutineINLANDSIS(SnoMod,BloMod,jjtime,debut)1 SUBROUTINE INLANDSIS(SnoMod,BloMod,jjtime,debut) 2 2 3 3 USE dimphy … … 415 415 . eta_SV(1,max(1,isnoSV(1)))<epsi) then 416 416 C + ********** 417 callSISVAT_BSn417 CALL SISVAT_BSn 418 418 endif 419 419 else 420 callSISVAT_BSn420 CALL SISVAT_BSn 421 421 C + ********** 422 422 endif … … 659 659 660 660 ! C.Agosta option for snow density, same as for BS i.e. 661 ! is_ok_density_kotlyakov=. false.662 c #BS density_kotlyakov = . false. !C.Amory BS 2018661 ! is_ok_density_kotlyakov=.FALSE. 662 c #BS density_kotlyakov = .FALSE. !C.Amory BS 2018 663 663 C + ... Fallen Snow Density, Adapted for Antarctica 664 664 if (is_ok_density_kotlyakov) then … … 898 898 if(isnoSV(1)>=1.or.NLaysv(1)>=1) then 899 899 C + ********** 900 callSISVAT_zSn900 CALL SISVAT_zSn 901 901 C + ********** 902 902 endif 903 903 else 904 904 C + ********** 905 callSISVAT_zSn905 CALL SISVAT_zSn 906 906 C + ********** 907 907 endif 908 908 909 909 C + ********** 910 ! #ve callSISVAT_wEq('_zSn ',0)910 ! #ve CALL SISVAT_wEq('_zSn ',0) 911 911 C + ********** 912 912 … … 981 981 982 982 ! ****** 983 callSnOptP(jjtime)983 CALL SnOptP(jjtime) 984 984 ! ****** 985 985 … … 996 996 ! Soil optical properties 997 997 ! ============================= 998 !Etienne: as in inlandis we do not callvgopt, we need to define998 !Etienne: as in inlandis we do not CALL vgopt, we need to define 999 999 !the albedo alb_SV and to calculate the 1000 1000 !absorbed Solar Radiation by Surfac (Normaliz)[-] SoSosv … … 1025 1025 1026 1026 ! ********** 1027 ! #ve callSISVAT_wEq('SnOptP',0)1027 ! #ve CALL SISVAT_wEq('SnOptP',0) 1028 1028 ! ********** 1029 1029 … … 1083 1083 if (iflag_temp_inlandsis == 0) then 1084 1084 1085 callSISVAT_TSo1085 CALL SISVAT_TSo 1086 1086 1087 1087 else … … 1090 1090 END DO 1091 1091 1092 callSISVAT_TS21092 CALL SISVAT_TS2 1093 1093 1094 1094 end if … … 1096 1096 1097 1097 ! ********** 1098 ! #ve callSISVAT_wEq('_TSo ',0)1098 ! #ve CALL SISVAT_wEq('_TSo ',0) 1099 1099 ! ********** 1100 1100 … … 1127 1127 1128 1128 ! ********** 1129 callSISVAT_qSn1129 CALL SISVAT_qSn 1130 1130 ! ********** 1131 1131 1132 1132 ! ********** 1133 ! #ve callSISVAT_wEq('_qSn ',0)1133 ! #ve CALL SISVAT_wEq('_qSn ',0) 1134 1134 ! ********** 1135 1135 … … 1168 1168 1169 1169 ! ********** 1170 callSISVAT_qSo1170 CALL SISVAT_qSo 1171 1171 ! #m0. (Wats_0,Wats_1,Wats_d) 1172 1172 ! ********** … … 1239 1239 if(isnoSV(1)>=1) then 1240 1240 C + ********** 1241 callSISVAT_GSn1241 CALL SISVAT_GSn 1242 1242 C + ********** 1243 1243 endif 1244 1244 else 1245 1245 C + ********** 1246 callSISVAT_GSn1246 CALL SISVAT_GSn 1247 1247 C + ********** 1248 1248 endif -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_bsn.F
r3900 r5103 1 1 2 2 3 subroutineSISVAT_BSn3 SUBROUTINE SISVAT_BSn 4 4 5 5 C +------------------------------------------------------------------------+ -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_gsn.F
r5082 r5103 1 1 2 subroutineSISVAT_GSn2 SUBROUTINE SISVAT_GSn 3 3 4 4 C +------------------------------------------------------------------------+ … … 296 296 C + ==== 297 297 298 data vector/. true./ ! Vectorization Switch298 data vector/.TRUE./ ! Vectorization Switch 299 299 data vdent1/ 0.5e8/ ! Wet Snow Metamorphism 300 300 cXF tuned for Greenland (2.e8=old value) -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_qsn.F
r5082 r5103 1 1 2 2 3 subroutineSISVAT_qSn3 SUBROUTINE SISVAT_qSn 4 4 . ( 5 5 ! #e1. EqSn_0,EqSn_1,EqSn_d -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_qso.F
r3792 r5103 1 1 2 2 3 subroutineSISVAT_qSo3 SUBROUTINE SISVAT_qSo 4 4 ! #m0. (Wats_0,Wats_1,Wats_d) 5 5 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_sno_albedo.F
r5082 r5103 1 subroutineSnOptP(jjtime)1 SUBROUTINE SnOptP(jjtime) 2 2 3 3 C +------------------------------------------------------------------------+ … … 464 464 ! prescription for each time step with NEMO values 465 465 466 c #AO if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. . true.) then466 c #AO if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. .TRUE.) then 467 467 c #AO if (AOmask(ikl) .eq. 0) then 468 468 c #AO albisv(ikl) = (1.-AOmask(ikl))* albAOsisv(ikl) -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_ts2.F
r5099 r5103 1 subroutine SISVAT_TS21 SUBROUTINE SISVAT_TS2 2 2 c #ES. (ETSo_0,ETSo_1,ETSo_d) 3 3 … … 475 475 ENDDO 476 476 477 end ! subroutine SISVAT_TS2477 end ! SUBROUTINE SISVAT_TS2 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.F
r5099 r5103 3 3 4 4 5 subroutineSISVAT_TSo5 SUBROUTINE SISVAT_TSo 6 6 ! #e1. (ETSo_0,ETSo_1,ETSo_d) 7 7 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_weq.F
r5082 r5103 1 1 2 2 3 subroutineSISVAT_wEq( labWEq ,istart)3 SUBROUTINE SISVAT_wEq( labWEq ,istart) 4 4 5 5 C +------------------------------------------------------------------------+ … … 48 48 49 49 IF (.NOT.logWEq) THEN 50 logWEq = . true.50 logWEq = .TRUE. 51 51 open(unit=45,status='unknown',file='SISVAT_wEq.ve') 52 52 rewind 45 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zag.F
r3792 r5103 1 1 2 2 3 subroutineSISVAT_zAg3 SUBROUTINE SISVAT_zAg 4 4 . (isagra,isagrb,WEagra 5 5 . ,dzagra,dzagrb,T_agra,T_agrb -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zcr.F
r3792 r5103 1 1 2 2 3 subroutineSISVAT_zCr3 SUBROUTINE SISVAT_zCr 4 4 C + 5 5 C +------------------------------------------------------------------------+ -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zsn.F
r5099 r5103 1 1 2 2 3 subroutineSISVAT_zSn3 SUBROUTINE SISVAT_zSn 4 4 5 5 C +------------------------------------------------------------------------+ … … 231 231 232 232 C + *************** 233 callSISVAT_zCr233 CALL SISVAT_zCr 234 234 C + *************** 235 235 … … 287 287 288 288 C + *************** 289 callSISVAT_zAg289 CALL SISVAT_zAg 290 290 . (isagr1,isagr2,WEagre 291 291 . ,dzagr1,dzagr2,T_agr1,T_agr2 … … 563 563 564 564 C + *************** 565 callSISVAT_zCr565 CALL SISVAT_zCr 566 566 C + *************** 567 567 … … 629 629 630 630 C + *************** 631 callSISVAT_zAg631 CALL SISVAT_zAg 632 632 . (isagr1,isagr2,WEagre 633 633 . ,dzagr1,dzagr2,T_agr1,T_agr2 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90
r5101 r5103 206 206 IF (debut) THEN 207 207 firstcall = .TRUE. 208 INI_SV = . false.208 INI_SV = .FALSE. 209 209 ELSE 210 firstcall = . false.211 INI_SV = . true.210 firstcall = .FALSE. 211 INI_SV = .TRUE. 212 212 END IF 213 213 … … 1047 1047 ! Objet: Lecture du fichier de conditions initiales pour SISVAT 1048 1048 !====================================================================== 1049 ! include "indicesol.h"1050 1051 ! include "dimsoil.h"1052 1049 include "clesphys.h" 1053 include "thermcell.h"1054 1050 include "compbl.h" 1055 1051 … … 1284 1280 IMPLICIT none 1285 1281 1286 ! include "indicesol.h"1287 ! include "dimsoil.h"1288 1282 include "clesphys.h" 1289 include "thermcell.h"1290 1283 include "compbl.h" 1291 1284
Note: See TracChangeset
for help on using the changeset viewer.