- Timestamp:
- 12/17/15 17:53:06 (9 years ago)
- Location:
- branches/iLoveclim/SOURCES
- Files:
-
- 6 deleted
- 12 edited
- 3 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/iLoveclim/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90
r22 r31 87 87 namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file 88 88 89 ! dmr moved up since namelist declaration should be in the header90 91 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time92 93 89 428 format(A) 94 90 rewind(num_param) ! pour revenir au debut du fichier param_list.dat … … 149 145 !lecture namelist 150 146 147 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time 151 148 152 149 rewind(num_param) ! pour revenir au debut du fichier param_list.dat -
branches/iLoveclim/SOURCES/Hemin15_files/module_choix-hemin15.f90
r25 r31 77 77 use deformation_mod_2lois ! module concernant les lois de deformation 78 78 79 !--------------Module propritete thermique de la glace80 ! use prop_therm_ice_heino81 82 79 ! eau basale 83 80 use eau_basale -
branches/iLoveclim/SOURCES/Makefile
r26 r31 10 10 11 11 # Choice for compiler (ifort = 1, gfortran = 0) 12 ifort ?= 012 ifort ?= 1 13 13 14 14 # Choice for librairies (mkl_c = 1 = MKL, = 0 = BLAS) -
branches/iLoveclim/SOURCES/Makefile.grisli-gfortran.inc
r26 r31 43 43 mod_clim_tof = climat_forcage_mois_mod.o climat-forcage-insolation_mod_oneway.o \ 44 44 climat-forcage-insolation_mod.o climat_GrIce2sea_years_mod.o \ 45 climat-perturb_mod-0.4.o \ 45 46 ablation_mod.o no_ablation_mod.o 46 47 … … 210 211 track_greeneem_mod.o 211 212 213 Liste_Ant40 = output_anta40_mod-0.4.o \ 214 lect-anteis_mod.o \ 215 bmelt-ant-regions_mod.o \ 216 fake-routines-ant_mod.o module_choix-antar40-0.4.o \ 217 track_ant40_mod.o 218 212 219 Liste_eura40 = output_eurasie40_mod-0.1.o \ 213 220 lect-eurasie_mod.o lakes-prescribed_mod-0.1.o \ … … 269 276 Dim_hemin15 = paradim-hemin15_mod.o geography-hemin15_mod.o 270 277 278 Dim_Ant40 = paradim-ant40_mod.o geography-ant40_mod.o 279 271 280 Dim_eura40 = paradim-euras40_mod.o geography-euras40_mod.o 272 281 … … 282 291 283 292 Dim_gsearise05 = paradim-gsearise05_mod.o geography-gsearise05_mod.o 284 285 293 286 294 # Liste_BLAS = band.o reduc.o \ … … 289 297 # BLAS/strsm.o BLAS/xerbla.o 290 298 291 Liste_BLAS = band.oreduc.o \299 Liste_BLAS = LAPACK/band.o LAPACK/reduc.o \ 292 300 BLAS/isamax.o BLAS/lsame.o BLAS/scopy.o \ 293 301 BLAS/slamch.o BLAS/slabad.o BLAS/slange.o \ … … 604 612 $(F_NETCDF) climat-forcage-insolation_mod.f90 605 613 614 climat-perturb_mod-0.4.o : climat-perturb_mod-0.4.f90 615 $(F_NETCDF) climat-perturb_mod-0.4.f90 616 606 617 ablation_mod.o : ablation_mod.f90 607 618 $(FT) ablation_mod.f90 … … 648 659 $(FT) accum_month_lapsecouche.f90 649 660 650 climat-perturb_mod-0.4.o : climat-perturb_mod-0.4.f90651 $(FT) climat-perturb_mod-0.4.f90652 653 661 climat-perturb-mois_mod.o : climat-perturb-mois_mod.f90 654 662 $(FT) climat-perturb-mois_mod.f90 … … 817 825 $(FT) Antarctique_general_files/output_anta_mod-0.4.f90 818 826 827 output_anta40_mod-0.4.o : Ant40_files/output_anta40_mod-0.4.f90 828 $(FT) Ant40_files/output_anta40_mod-0.4.f90 829 819 830 output_global_mod.o :output_global_mod.f90 820 831 $(FT) output_global_mod.f90 … … 1023 1034 $(FT) lastoutput-0.2.f 1024 1035 1025 lect-eis.o : lect-eis. F1026 $(FT) lect-eis. F1036 lect-eis.o : lect-eis.f90 1037 $(FT) lect-eis.f90 1027 1038 1028 1039 lect_datfile.o : lect_datfile.f90 … … 1035 1046 $(FT) limit_file.f90 1036 1047 1037 lineartemp-0.2.o : lineartemp-0.2. F#module3d.mod1038 $(FT) lineartemp-0.2. F1048 lineartemp-0.2.o : lineartemp-0.2.f90 1049 $(FT) lineartemp-0.2.f90 1039 1050 1040 1051 litho-0.3.o : litho-0.3.f90#module3d.mod iso_declar.mod … … 1340 1351 $(FT) New-remplimat/graph-sgbsv.f90 1341 1352 1342 # Routines venant de bibliotheques exterieures 1343 reduc.o : reduc.f#module3d.mod 1344 $(FT) $(F_90) reduc.f 1345 1346 # band vient de Lapack 1347 band.o: band.f 1348 $(FT) $(F_90) band.f 1353 # LAPACK 1354 reduc.o : LAPACK/reduc.f 1355 $(FT) $(F_90) LAPACK/reduc.f 1356 1357 band.o: LAPACK/band.f 1358 $(FT) $(F_90) LAPACK/band.f 1349 1359 1350 1360 #BLAS … … 1678 1688 $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) 1679 1689 1690 Ant-40 : $(Dim_Ant40) $(mod_dim_communs) \ 1691 $(toy_recul) \ 1692 $(mod_communs) \ 1693 $(mod_clim_tof) \ 1694 $(mod_no_tracers) \ 1695 $(mod_ell) $(Liste_Ant40) \ 1696 $(diagnoshelf) \ 1697 $(Liste_Netcdf) \ 1698 $(routines_communes) steps_time_loop.o \ 1699 $(routine_elliptiques) \ 1700 $(Liste_BLAS) 1701 1702 $(LK) -o ../bin/Ant-40 \ 1703 $(Dim_Ant40) $(mod_dim_communs) \ 1704 $(toy_recul) \ 1705 $(mod_communs) \ 1706 $(mod_clim_tof) \ 1707 $(mod_no_tracers) \ 1708 $(mod_ell) $(Liste_Ant40) \ 1709 $(diagnoshelf) \ 1710 $(Liste_Netcdf) \ 1711 $(routines_communes) steps_time_loop.o \ 1712 $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) $(Liste_BLAS) 1713 1680 1714 1681 1715 -
branches/iLoveclim/SOURCES/Makefile.grisli.inc
r27 r31 297 297 # BLAS/strsm.o BLAS/xerbla.o 298 298 ifeq ($(mkl_c), 0) 299 Liste_BLAS = band.oreduc.o \299 Liste_BLAS = LAPACK/band.o LAPACK/reduc.o \ 300 300 BLAS/isamax.o BLAS/lsame.o BLAS/scopy.o \ 301 301 BLAS/slamch.o BLAS/slabad.o BLAS/slange.o \ … … 1034 1034 $(FT) lastoutput-0.2.f 1035 1035 1036 lect-eis.o : lect-eis. F1037 $(FT) lect-eis. F1036 lect-eis.o : lect-eis.f90 1037 $(FT) lect-eis.f90 1038 1038 1039 1039 lect_datfile.o : lect_datfile.f90 … … 1046 1046 $(FT) limit_file.f90 1047 1047 1048 lineartemp-0.2.o : lineartemp-0.2. F#module3d.mod1049 $(FT) lineartemp-0.2. F1048 lineartemp-0.2.o : lineartemp-0.2.f90 1049 $(FT) lineartemp-0.2.f90 1050 1050 1051 1051 litho-0.3.o : litho-0.3.f90#module3d.mod iso_declar.mod … … 1351 1351 $(FT) New-remplimat/graph-sgbsv.f90 1352 1352 1353 # Routines venant de bibliotheques exterieures 1354 reduc.o : reduc.f#module3d.mod 1355 $(FT) $(F_90) reduc.f 1356 1357 # band vient de Lapack 1358 band.o: band.f 1359 $(FT) $(F_90) band.f 1353 # LAPACK 1354 reduc.o : LAPACK/reduc.f 1355 $(FT) $(F_90) LAPACK/reduc.f 1356 1357 band.o: LAPACK/band.f 1358 $(FT) $(F_90) LAPACK/band.f 1360 1359 1361 1360 ifeq ($(mkl_c), 0) -
branches/iLoveclim/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90
r23 r31 114 114 integer,dimension(:),allocatable :: num_ncdf_file !< compteur des fichiers netcdf par class 115 115 116 real*8, dimension(:,:), pointer :: tab 117 real*8, dimension(:,:,:), pointer :: tab1 118 real*8, dimension(:,:,:), pointer :: tab1T 116 real*8, dimension(:,:), pointer :: tab => null() !< tableau 2d real ecrit dans le fichier 117 real*8, dimension(:,:,:), pointer :: tab1 => null() !< tableau 3d real 118 real*8, dimension(:,:,:), pointer :: tab1T => null() !< tableau 3d real pour la temperature 119 119 120 120 … … 653 653 real (kind=kind(0.d0)) :: timetmp !< variable intermediaire 654 654 character(len=20) :: nametmp !< nom intermediaire 655 real*8,pointer,dimension(:) :: liste_time => null() !< liste des snapshot des variables ecrites en netcdf 656 real*8,pointer,dimension(:) :: x,y,x1,y1,z,nzzm 657 real*8,pointer,dimension(:,:):: lat,lon 655 real*8,pointer,dimension(:) :: liste_time => null() !< liste des snapshot des variables ecrites en netcdf 656 real*8,pointer,dimension(:) :: x,y,x1,y1,z,nzzm 657 real*8,pointer,dimension(:,:):: lat,lon => null() 658 658 integer :: i,j,l,k,p 659 659 logical :: fait … … 674 674 nametmp = 'f_' 675 675 endif 676 677 allocate(tab(nx,ny),tab1(nx,ny,nz),tab1T(nx,ny,nz+nzm)) 676 677 if (.not.associated(tab)) allocate(tab(nx,ny)) 678 if (.not.associated(tab1)) allocate(tab1(nx,ny,nz)) 679 if (.not.associated(tab1T)) allocate(tab1T(nx,ny,nz+nzm)) 678 680 679 681 if (.not.associated(liste_time)) then -
branches/iLoveclim/SOURCES/Temperature-routines/icetemp_declar_mod.f90
r25 r31 55 55 ! ___________________________ 56 56 Real,allocatable,Dimension(:,:):: Tbmer !< Temperature De La Mer A La Base De L'Ice Shelf 57 ! Real,allocatable,Dimension(:,:) :: Alpha !< Pente Locale Sur Noeud Majeur58 ! Real,allocatable,Dimension(:,:) :: Ubar !< Vitesse Sur Noeud Majeur59 57 Real,allocatable,Dimension(:,:) :: Chalglissx,Chalglissy !< Chaleur De Glissement 60 58 Integer,allocatable,Dimension(:,:) :: Iadvec_w,Iadvec_e,Iadvec_s,Iadvec_n -
branches/iLoveclim/SOURCES/Temperature-routines/icetemp_mod.f90
r25 r31 85 85 !Tab Nx,Ny 86 86 Allocate(Tbmer(Nx,Ny),& 87 ! Alpha(Nx,Ny),&88 ! Ubar(Nx,Ny),&89 87 Chalglissx(Nx,Ny),& 90 88 Chalglissy(Nx,Ny),& … … 118 116 Rr=0. 119 117 Hh=0. 120 Tdot=0 118 Tdot=0 121 119 Chal2_x=0. 122 120 Chal2_y=0. … … 298 296 299 297 Deallocate (Aa,Bb,Cc,Rr,Hh,Tdot,Abis,Bbis,Cbis,Rbis,Hbis,Ee,Tbmer, & 300 Chalglissx,Chalglissy,Iadvec_w,Iadvec_e, & !Alpha,Ubar298 Chalglissx,Chalglissy,Iadvec_w,Iadvec_e, & 301 299 Iadvec_s,Iadvec_n,Chal2_x,Chal2_y,Chal2_z,Chal2_xy,Chaldef_maj, & 302 300 Advecx,Advecy,Advec,Chalx,Chaly,Ffx,Ffy,T3d_new) -
branches/iLoveclim/SOURCES/ablation_mod.f90
r17 r31 31 31 real, parameter :: PYG=2*PI_L/NYEAR ! ct for PDD calculation (PYG remplace PY qui existe dans CLIMBER 32 32 REAL, dimension(nx,ny) :: PDDCT ! ct for PDD calculation 33 REAL, dimension(nx,ny) :: PDDCT2 ! ct for PDD calculation 33 REAL, dimension(nx,ny) :: PDDCT2 ! ct for PDD calculation 34 integer :: methode_abl=0 ! selection methode pdd (0) ou van den Berg 2008 (1) 34 35 35 36 contains … … 83 84 84 85 85 USE module3d_phy,only:Tjuly,Tann,Tmois,acc,pdd,TS,Tshelf,precip,BM, S,dice,cl86 USE module3d_phy,only:Tjuly,Tann,Tmois,acc,pdd,TS,Tshelf,precip,BM,Abl,S,dice,cl 86 87 87 88 IMPLICIT NONE … … 128 129 ENDWHERE 129 130 130 131 ELSEIF (pdd_type.EQ.2) THEN ! pdd Tarasov 131 132 132 133 where (TJULY(:,:).gt.10.) 133 134 134 Cice_2D(:,:) = 8.3*1e-3 135 Csnow_2D(:,:) = 4.3*1e-3 135 136 endwhere 136 137 where(TJULY(:,:).gt.-1..and.TJULY(:,:).lt.10.) 137 138 138 Cice_2D(:,:) = 1e-3*(8.3+0.0067*(10.-TJULY(:,:))**3) 139 Csnow_2D(:,:) = 1e-3*(2.8+0.15*TJULY(:,:)) 139 140 endwhere 140 141 where(TJULY(:,:).le.-1.) 141 142 142 Cice_2D(:,:) = 17.22*1e-3 143 Csnow_2D(:,:) = 2.65*1e-3 143 144 endwhere 144 145 145 146 sigma_ice_2D(:,:)=sigma_ice 146 147 ENDIF 147 148 148 149 … … 187 188 188 189 ! calcul du Bilan de masse 189 IF (pdd_type.EQ.1) THEN 190 PDS(:,:)=ACC(:,:)/Csnow_2D(:,:) 191 SIMAX(:,:)=ACC(:,:)*CSI_2D(:,:) 192 PDSI(:,:)=SIMAX(:,:)/Cice_2D(:,:) 190 IF (methode_abl.EQ.0) THEN 191 IF (pdd_type.EQ.1) THEN 192 PDS(:,:)=ACC(:,:)/Csnow_2D(:,:) 193 SIMAX(:,:)=ACC(:,:)*CSI_2D(:,:) 194 PDSI(:,:)=SIMAX(:,:)/Cice_2D(:,:) 193 195 ! avec regel de 60% puis fonte (2 premiers where) : 194 196 ! WHERE (PDD(:,:).LE.CSI_2D*PDS(:,:)) … … 200 202 ! SIF(:,:)=SIMAX(:,:) 201 203 ! endwhere 202 WHERE (PDD(:,:).LE.PDS(:,:)) ! test avec regel de 60% progressif203 BM(:,:)=ACC(:,:)-PDD(:,:)*Csnow_2D*(1-CSI_2D(:,:))204 SIF(:,:)=PDD(:,:)*Csnow_2D(:,:)*(1-CSI_2D(:,:))205 endwhere206 WHERE ((PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:)+PDSI(:,:)))207 BM(:,:)=SIMAX(:,:)-(PDD(:,:)-PDS(:,:))*Cice_2D208 SIF(:,:)=SIMAX(:,:)209 endwhere210 WHERE (PDS(:,:)+PDSI(:,:).LE.PDD(:,:))211 BM(:,:)=(PDS(:,:)+PDSI(:,:)-PDD(:,:))*Cice_2D212 SIF(:,:)=SIMAX(:,:)213 endwhere214 ELSEIF (PDD_type.EQ.2) THEN ! PDD Tarasov204 WHERE (PDD(:,:).LE.PDS(:,:)) ! test avec regel de 60% progressif 205 BM(:,:)=ACC(:,:)-PDD(:,:)*Csnow_2D*(1-CSI_2D(:,:)) 206 SIF(:,:)=PDD(:,:)*Csnow_2D(:,:)*(1-CSI_2D(:,:)) 207 endwhere 208 WHERE ((PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:)+PDSI(:,:))) 209 BM(:,:)=SIMAX(:,:)-(PDD(:,:)-PDS(:,:))*Cice_2D 210 SIF(:,:)=SIMAX(:,:) 211 endwhere 212 WHERE (PDS(:,:)+PDSI(:,:).LE.PDD(:,:)) 213 BM(:,:)=(PDS(:,:)+PDSI(:,:)-PDD(:,:))*Cice_2D 214 SIF(:,:)=SIMAX(:,:) 215 endwhere 216 ELSEIF (PDD_type.EQ.2) THEN ! PDD Tarasov 215 217 PDS(:,:)=ACC(:,:)/Csnow_2D(:,:) 216 218 pr_ice_eq(:,:) = amax1(0.,((PRECIP(:,:)/DICE)-ACC(:,:))) ! precipe liquide (ice equivalent) … … 231 233 232 234 ! where(snowmelt(:,:).lt.ACC(:,:)) 233 234 235 refr2(:,:) = 2.2*(ACC(:,:)-snowmelt(:,:))-(cpsurf(:,:)/CL)*amin1(TANN(:,:),0.) 236 refreezed_ice(:,:) = amin1(pr_ice_eq(:,:)+snowmelt(:,:),refr2(:,:)) 235 237 ! elsewhere 236 238 ! refr2(:,:) = -(cpsurf(:,:)/CL)*amin1(TANN(:,:),0.) … … 252 254 endwhere 253 255 254 ELSE ! pdd standard reeh256 ELSE ! pdd standard reeh 255 257 ! (* Positive degrees required to melt the snow layer *) 256 258 PDS(:,:)=ACC(:,:)/Csnow_2D(:,:) … … 260 262 PDSI(:,:)=SIMAX(:,:)/Cice_2D(:,:) 261 263 ! avec regel de 60% puis fonte (2 premiers where) : 262 WHERE (PDD(:,:).LE.CSI_2D*PDS(:,:))263 BM(:,:)=ACC(:,:)264 SIF(:,:)=PDD(:,:)*Csnow_2D(:,:)265 endwhere266 WHERE ((CSI_2D*PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:)))267 BM(:,:)=ACC(:,:)+SIMAX(:,:)-PDD(:,:)*Csnow_2D268 SIF(:,:)=SIMAX(:,:)269 endwhere264 WHERE (PDD(:,:).LE.CSI_2D*PDS(:,:)) 265 BM(:,:)=ACC(:,:) 266 SIF(:,:)=PDD(:,:)*Csnow_2D(:,:) 267 endwhere 268 WHERE ((CSI_2D*PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:))) 269 BM(:,:)=ACC(:,:)+SIMAX(:,:)-PDD(:,:)*Csnow_2D 270 SIF(:,:)=SIMAX(:,:) 271 endwhere 270 272 ! WHERE (PDD(:,:).LE.PDS(:,:)) ! test avec regel de 60% progressif| 271 273 ! BM(:,:)=ACC(:,:)-PDD(:,:)*Csnow_2D*(1-CSI_2D) ! remplace les 2 premiers where | 272 274 ! SIF(:,:)=PDD(:,:)*Csnow_2D(:,:)*(1-CSI_2D(:,:)) !----------------------------------| 273 275 ! endwhere 274 WHERE ((PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:)+PDSI(:,:))) 275 BM(:,:)=SIMAX(:,:)-(PDD(:,:)-PDS(:,:))*Cice_2D 276 SIF(:,:)=SIMAX(:,:) 277 endwhere 278 WHERE (PDS(:,:)+PDSI(:,:).LE.PDD(:,:)) 279 BM(:,:)=(PDS(:,:)+PDSI(:,:)-PDD(:,:))*Cice_2D 280 SIF(:,:)=SIMAX(:,:) 281 endwhere 282 276 WHERE ((PDS(:,:).LT.PDD(:,:)).AND.(PDD(:,:).LE.PDS(:,:)+PDSI(:,:))) 277 BM(:,:)=SIMAX(:,:)-(PDD(:,:)-PDS(:,:))*Cice_2D 278 SIF(:,:)=SIMAX(:,:) 279 endwhere 280 WHERE (PDS(:,:)+PDSI(:,:).LE.PDD(:,:)) 281 BM(:,:)=(PDS(:,:)+PDSI(:,:)-PDD(:,:))*Cice_2D 282 SIF(:,:)=SIMAX(:,:) 283 endwhere 284 ENDIF 285 ELSEIF ( methode_abl.EQ.1 ) THEN ! Insolation Temperature Melt equation, van den Berg, 2008 286 SIMAX(:,:)=ACC(:,:)*CSI 287 SIF(:,:)=SIMAX(:,:) 288 ELSE 289 print*,'ablation.f90 : pb methode calcul BM' 290 print*,'ATTENTION methode_abl > 1 NON COMPATIBLE AVEC iLOVECLIM' 291 stop 283 292 ENDIF 284 293 285 294 ! calcul de la temperature de surface (utilisee dans icetemp) : 286 TS(:,:)=(TANN(:,:)+26.6*SIF(:,:)) 287 TS(:,:)=min(0.0,TS(:,:)) 288 tshelf(:,:)=TS(:,:) 289 290 !!$! (*** ablation ***) 291 !!$ DO i=1,nx 292 !!$ DO j=1,ny 293 !!$! (* positive degrees required to melt the snow layer *) 294 !!$ pds=acc(i,j)/csnow 295 !!$! (* maximum amount of super. ice that can be formed *) 296 !!$ simax=acc(i,j)*csi 297 !!$! (* pos. degrees required to melt the superimposed ice *) 298 !!$ pdsi=simax/cice 299 !!$ IF (pdd(i,j).LE.csi*pds) THEN 300 !!$ sif=pdd(i,j)*csnow 301 !!$ ELSE 302 !!$ sif=simax 303 !!$ ENDIF 304 !!$ 305 !!$! surface temperature 306 !!$ ts(i,j)=(tann(i,j)+26.6*sif) 307 !!$ ts(i,j)=MIN(0.0,ts(i,j)) 308 !!$ tshelf(i,j)=ts(i,j) 309 !!$ 310 !!$! mass balance 311 !!$ IF (pdd(i,j).LE.csi*pds) bm(i,j)=acc(i,j) ! toute la fonte est recuperee en regel 312 !!$ IF ((csi*pds.LT.pdd(i,j)).AND.(pdd(i,j).LE.pds)) & ! fonte de la neige uniquement 313 !!$ bm(i,j)=acc(i,j)+simax-pdd(i,j)*csnow 314 !!$ IF ((pds.LT.pdd(i,j)).AND.(pdd(i,j).LE.pds+pdsi)) & ! fonte toute neige et une partie glace regel 315 !!$ bm(i,j)=simax-(pdd(i,j)-pds)*cice 316 !!$ IF (pds+pdsi.LE.pdd(i,j)) bm(i,j)=(pds+pdsi-pdd(i,j))*cice !fonte de tt neige tt regel + glace annees precedentes 317 318 319 !940 format('%%%% ',a,' time=',f8.0,' %%%%') 295 TS(:,:)=(TANN(:,:)+26.6*SIF(:,:)) 296 TS(:,:)=min(0.0,TS(:,:)) 297 tshelf(:,:)=TS(:,:) 298 Abl(:,:)=BM(:,:)-Acc(:,:) 299 320 300 END SUBROUTINE ABLATION 321 301 -
branches/iLoveclim/SOURCES/initial-0.3.f90
r24 r31 71 71 endif 72 72 73 call sortie_ncdf_cat73 !cdc call sortie_ncdf_cat 74 74 75 75 ! calcul de Hmx et Hmy -> shift=-1, dim=1 -> H(i-1,j) -
branches/iLoveclim/SOURCES/spinup_mod.f90
r22 r31 209 209 coef_defmx(i,j) = 1. 210 210 Uxbar(i,j) = Vcol_x(i,j) 211 flgzmx(i,j) = transfer(Vcol_x(i,j),flgzmx(i,j)) 211 ! dmr below is a cast from a real*4 to logical*4 212 ! dmr cannot be implicit in gfortran 213 ! dmr flgzmx(i,j) = Vcol_x(i,j) 214 flgzmx(i,j) = (nint(Vcol_x(i,j)).ne.0) 212 215 uxdef(i,j) = 0. 213 216 Ubx(i,j) = Vcol_x(i,j) … … 250 253 coef_defmy(i,j) = 1. 251 254 Uybar(i,j) = Vcol_y(i,j) 252 flgzmy(i,j) = transfer(Vcol_y(i,j),flgzmy(i,j)) 255 ! dmr below is a cast from a real*4 to logical*4 256 ! dmr cannot be implicit in gfortran 257 ! dmr flgzmy(i,j) = Vcol_y(i,j) 258 flgzmy(i,j) = (nint(Vcol_y(i,j)).ne.0) 259 253 260 uydef(i,j) = 0. 254 261 Uby(i,j) = Vcol_y(i,j) -
branches/iLoveclim/SOURCES/util_recovery.f90
r10 r31 47 47 read(num_file,*) !saut de la ligne "-------" 48 48 read(num_file,*) i 49 allocate(tab_time(i)) 49 if (.not. allocated(tab_time)) then ! aurel dec15 50 allocate(tab_time(i)) 51 end if 50 52 do j=1,i 51 53 read(num_file,*) tab_time(j)
Note: See TracChangeset
for help on using the changeset viewer.