Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (7 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARdSV.F90
r5087 r5099 97 97 INTEGER :: iso1 !rajout hjp for ini - check if ivg,iso is possible 98 98 99 !100 101 99 END MODULE VARdSV -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARphy.F90
r5087 r5099 8 8 ! Etienne, 08112020: merging YOMCST_SISVAT, PHY_SV and VARphy 9 9 ! but redundancies remain 10 ! 10 11 11 ! ATTENTION!!!!: ce fichier include est compatible format fixe/format libre 12 12 ! veillez à n'utiliser que des ! pour les commentaires 13 13 ! et à bien positionner les & des lignes de continuation 14 14 ! (les placer en colonne 6 et en colonne 73) 15 ! 16 ! 15 16 17 17 ! A1.0 Fundamental constants 18 18 !---------------------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARxSV.F90
r3900 r5099 201 201 INTEGER,DIMENSION(nb_wri),SAVE :: lwriSV ! OUTPUT point vec Index 202 202 !$OMP THREADPRIVATE(lwriSV) 203 ! 203 204 204 INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: ii__SV ! WORK point i Coordinate 205 205 !$OMP THREADPRIVATE(ii__SV) … … 329 329 ! +--SISVAT OUTPUT Variables 330 330 ! + ----------------------------- 331 332 ! 331 333 332 ALLOCATE(ii__SV(klonv)) ! WORK point i Coordinate 334 333 ALLOCATE(jj__SV(klonv)) ! WORK point j Coordinate -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARySV.F90
r3900 r5099 86 86 ALLOCATE(alb6sv(klonv,6))! 6-band Albedo 87 87 88 !89 88 ALLOCATE(albssv(klonv)) ! Soil Albedo [-] 90 89 ALLOCATE(SoSosv(klonv)) !? Surface Absorbed Solar Radiat. -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/inlandsis.F
r5082 r5099 1031 1031 ! Surface Emissivity (Etienne: simplified calculation for landice) 1032 1032 ! ============================================================= 1033 ! 1033 1034 1034 DO ikl=1,knonv 1035 1035 LSnMsk = min( 1,isnoSV(ikl)) -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_ts2.F
r5082 r5099 126 126 127 127 ! write(*,*)'T check' 128 ! 128 129 129 ! DO ig = 1,knonv 130 130 ! DO jk = 1,isnoSV(ig) !nsno … … 132 132 ! TsisSV(ig,jk) = TsisSV(ig,isnoSV(ig)) 133 133 ! ENDIF 134 ! 134 135 135 ! IF (TsisSV(ig,jk) <= 1.) THEN !hj check 136 136 ! TsisSV(ig,jk) = 273.15 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.F
r5093 r5099 467 467 468 468 if(ist__s==1) then ! to reduce computer time 469 ! 469 470 470 DO it_srf=1,nt_srf ! 471 471 dt_ver = dt_ver +dt_srf ! … … 481 481 DO ikl=1,knonv 482 482 ist = max(0,isotSV(ikl)-100*isnoSV(ikl)) ! 0 if H2O 483 ! 483 484 484 Psi(ikl) = ! 485 485 . psidSV(ist) ! DR97, Eqn 3.34 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zsn.F
r5082 r5099 162 162 . isnoSV(ikl)-isn +1 )) ! Snow Pack 163 163 . * min(unun, ! 164 ! 164 165 165 ! 1st additional Condition to accept OKthin 166 166 . max(zero, ! combination 167 167 . sign(unun,G1snSV(ikl, isn ) ! G1 with same 168 168 . *G1snSV(ikl,max(1,isn-1)))) ! sign => OK 169 ! 169 170 170 ! 2nd additional Condition to accept OKthin 171 171 . + max(zero, ! G1>0 172 172 . sign(unun,G1snSV(ikl, isn ))) ! =>OK 173 ! 173 174 174 ! 3rd additional Condition to accept OKthin 175 175 . + max(zero, ! dz too small -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90
r5087 r5099 593 593 END DO 594 594 595 !596 595 ! + +++ READ FORCINGS: END +++ 597 596 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ … … 751 750 ! to LMDZ. In LMDZ it is done in soil.F, which is not used in the case 752 751 ! of SISVAT, therefore it's needed here. 753 ! 752 754 753 USE mod_phys_lmdz_mpi_data, ONLY : is_mpi_root 755 754 USE mod_phys_lmdz_para … … 862 861 !C | | 863 862 !C +------------------------------------------------------------------------+ 864 ! 865 ! 863 866 864 867 865 !C +--Global Variables … … 933 931 !c #kd END DO 934 932 !c #kd END IF 935 ! 933 936 934 ! IF (nsol.ne.4) THEN 937 935 ! DO isl= 0,-nsol,-1 … … 941 939 ! . + misl_2 * 0.003) * 10**(nisl_2)) * 4. 942 940 !C +... dz_dSV(0) = Hapex-Sahel Calibration: 4 mm 943 ! 941 944 942 !c +SH dz_dSV(isl)=(((1-misl_2) * 0.001 945 943 !c +SH. + misl_2 * 0.003) * 10**(nisl_2)) * 1. 946 ! 944 947 945 !c #05 dz_dSV(isl)=(((1-misl_2) * 0.001 948 946 !c #05. + misl_2 * 0.008) * 10**(nisl_2)) * 0.5
Note: See TracChangeset
for help on using the changeset viewer.