- Timestamp:
- 01/19/16 17:01:36 (8 years ago)
- Location:
- branches/iLoveclim
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/iLoveclim
-
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
Note: See TracChangeset
for help on using the changeset viewer.