- Timestamp:
- 01/19/16 17:01:36 (8 years ago)
- Location:
- branches/iLoveclim
- Files:
-
- 12 edited
- 2 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/iLoveclim
-
branches/iLoveclim/SOURCES/Ant40_files/lect-anteis_mod.f90
r27 r35 115 115 ! endwhere 116 116 117 118 where (BSOC(:,:).LT.-9999.) 119 BSOC(:,:)=-9999. 120 endwhere 121 117 122 sealevel0=0. ! voir a passer dans le fichier parametre 118 123 S0(:,:)=max(S0(:,:),sealevel0) ! pour etre au niveau des mers : ATTENTION si SEALEV <0 … … 231 236 open(88,file=ghf_fich) 232 237 233 write(num_rep_42,*) 'flux geothermique Shapiro : ', TRIM(DIRNAMEINP)//'ijphi-40km-ant.dat' 234 235 do k=1,nx*ny 236 read(88,*) i,j,ghf(i,j) 237 end do 238 close(88) 238 ! write(num_rep_42,*) 'flux geothermique Shapiro : ', TRIM(DIRNAMEINP)//'ijphi-40km-ant.dat' 239 ! do k=1,nx*ny 240 ! read(88,*) i,j,ghf(i,j) 241 ! end do 242 ! close(88) 243 244 call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf) ! socle 245 246 239 247 240 248 ! pour passer les flux des mW/m2 au J/m2/an -
branches/iLoveclim/SOURCES/Ant40_files/module_choix-antar40-0.4.f90
r27 r35 45 45 46 46 !--------------Lecture climat forcage--------------- 47 !use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 48 use climat_perturb_mod ! climat perturbe a reverifier Dec 2015 47 use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 48 !use climat_perturb_mod ! climat perturbe a reverifier Dec 2015 49 !use climat_Grice2sea_years_mod ! climat force par fichier SMB directement (grice2sea) 49 50 50 51 !use climat_forcage_mod … … 105 106 106 107 !--------------Fusion basale------------------------ 107 use bmelt_ant_regions ! pour l'Antarctique avec régions 108 !use bmelt_ant_regions ! pour l'Antarctique avec régions 109 use bmelt_ant_regions_oce ! pour Ant region et temp ocean 108 110 !use bmelt_nor_regions ! pour le nord avec régions 109 111 !use bmelt_nor_depth ! pour le nord avec profondeur d'eau -
branches/iLoveclim/SOURCES/Fichiers-parametres/anteis1_LISTE-VAR-NETCDF.dat
r27 r35 144 144 1 1 1 145 145 ----------------------------------------------------------- 146 ghf 147 23 148 1 1 -1 149 ----------------------------------------------------------- 146 150 phid 147 151 24 -
branches/iLoveclim/SOURCES/Fichiers-parametres/greeneem15_LISTE-VAR-NETCDF.dat
r19 r35 144 144 1 1 1 145 145 ----------------------------------------------------------- 146 ghf 147 23 148 1 1 -1 149 ----------------------------------------------------------- 146 150 phid 147 151 24 -
branches/iLoveclim/SOURCES/Fichiers-parametres/hemin15_LISTE-VAR-NETCDF.dat
r14 r35 144 144 1 1 1 145 145 ----------------------------------------------------------- 146 ghf 147 23 148 1 1 -1 149 ----------------------------------------------------------- 146 150 phid 147 151 24 -
branches/iLoveclim/SOURCES/Fichiers-parametres/hemin40_LISTE-VAR-NETCDF.dat
r12 r35 144 144 1 1 1 145 145 ----------------------------------------------------------- 146 ghf 147 23 148 1 1 -1 149 ----------------------------------------------------------- 146 150 phid 147 151 24 -
branches/iLoveclim/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90
r31 r35 16 16 17 17 18 use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,dirnameinp 18 use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,dirnameinp,coefbmshelf 19 19 !use lect_climref_Ice2sea 20 20 use netcdf … … 86 86 87 87 namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file 88 ! dmr --- the namelist should be at the header position (moved up from below) 89 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time 88 90 89 91 428 format(A) … … 102 104 smb_file = trim(dirnameinp)//trim(smb_file) 103 105 104 call Read_Ncdf_var(' z',smb_file,tab)106 call Read_Ncdf_var('smb',smb_file,tab) 105 107 106 108 ! call lect_input(3,'smb',1,bm,smb_file,trim(dirnameinp)//trim(runname)//'.nc') … … 130 132 temp_annual_file = trim(dirnameinp)//trim(temp_annual_file) 131 133 132 call Read_Ncdf_var(' z',temp_annual_file,tab)134 call Read_Ncdf_var('Tann',temp_annual_file,tab) 133 135 Tann(:,:) = tab(:,:) 134 136 ! call lect_input(3,'Tann',1,Tann,temp_annual_file,trim(dirnameinp)//trim(runname)//'.nc') … … 145 147 !lecture namelist 146 148 147 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time 149 ! dmr --- namelist should be in the header part in the standard ... 150 ! namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time 148 151 149 152 rewind(num_param) ! pour revenir au debut du fichier param_list.dat … … 281 284 call init_grad_smb 282 285 endif 286 ! coefbmshelf coefficient pour la fusion basale sous les ice shelves 287 ! ici pas d'evolution temporelle donc coefbmshelf est fixe 288 coefbmshelf=1. 283 289 284 290 return -
branches/iLoveclim/SOURCES/Makefile.grisli.inc
r31 r35 214 214 lect-anteis_mod.o \ 215 215 bmelt-ant-regions_mod.o \ 216 bmelt-ant-regions-oce_mod.o \ 216 217 fake-routines-ant_mod.o module_choix-antar40-0.4.o \ 217 218 track_ant40_mod.o … … 515 516 516 517 bmelt-ant-regions_mod.o : Ant40_files/bmelt-ant-regions_mod.f90 517 $(FT) Ant40_files/bmelt-ant-regions_mod.f90 518 $(FT) Ant40_files/bmelt-ant-regions_mod.f90 519 520 bmelt-ant-regions-oce_mod.o : Ant40_files/bmelt-ant-regions-oce_mod.f90 521 $(FT) Ant40_files/bmelt-ant-regions-oce_mod.f90 518 522 519 523 bmelt-nor-regions_mod.o : Heminord_files/bmelt-nor-regions_mod.f90 -
branches/iLoveclim/SOURCES/Netcdf-routines/Description_Variables.dat
r4 r35 225 225 "celsius" 226 226 "Basal temperature " 227 ------------------------------------------------------ 228 "xxx" 229 23 "ghf" 1 230 "o" 231 "geothermal_heat_flux" 232 "geothermal_heat_flux" 233 "m.w/m²" 234 "Geothermal heat flux" 227 235 ------------------------------------------------------ 228 236 "xxx" -
branches/iLoveclim/SOURCES/initial-0.3.f90
r32 r35 73 73 endif 74 74 75 !cdccall sortie_ncdf_cat75 call sortie_ncdf_cat 76 76 77 77 ! calcul de Hmx et Hmy -> shift=-1, dim=1 -> H(i-1,j) -
branches/iLoveclim/SOURCES/lect-eis.f90
r31 r35 18 18 19 19 implicit none 20 character(len= 17) :: DIRNAMEINP20 character(len=20) :: DIRNAMEINP 21 21 integer nx,ny,i,j,kdum 22 22 real data(nx,ny)
Note: See TracChangeset
for help on using the changeset viewer.