Changeset 68
- Timestamp:
- 06/03/16 16:06:35 (8 years ago)
- Location:
- trunk/SOURCES
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/ANT15-LBq_files/module_choix_antar15_LBq.f90
r4 r68 36 36 ! use lect_topo_anteis ! pour l'antarctique 40 km 37 37 38 use lect_topo_ant_gen ! read antarctic topo with .dat files 38 !use lect_topo_ant_gen ! read antarctic topo with .dat files Catherine 39 use lect_topo_anteis ! pour l'antarctique 40 km ou 15km Christophe 39 40 !use lect_topo_ant_CISM_gen 40 41 !use lect_topo_nord40 ! pour l'hemisphere nord 40 km … … 45 46 46 47 ! use lect_clim_act_ant_gen 47 use lect_clim_acc_T_ant_gen 48 ! use lect_clim_acc_T_ant_gen ! module Catherine climat Antarctique 49 50 !--------------Module climat --------------- 51 !use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 52 !use climat_Grice2sea_years_mod ! climat force par fichier SMB directement (grice2sea) 53 use climat_Grice2sea_years_perturb_mod ! forcage avec SMB type MAR et evolutionn temporelle (futur ou carotte) 48 54 49 55 ! use lect_clim_act_nord40 ! pour l'hemisphere nord et l'eurasie … … 51 57 52 58 !--------------Lecture climat forcage--------------- 53 use climat_perturb_mod59 !use climat_perturb_mod 54 60 !use climat_forcage_mod 55 61 !use climat_synthes_mod 56 62 !use climat_profil_mod 57 63 !use climat_regions_delta 64 65 !use ablation_mod ! calcul de l'ablation (PDD ou autre methode) 66 use no_ablation ! pas de calcul de l'ablation => lecture fichier SMB (necessaire avec climat_Grice2sea_years_mod) 58 67 59 68 ! pas de lacs proglaciaires … … 64 73 ! attention tracer fait appel a un modele climat qui doit etre compatible avec celui choisi 65 74 ! dans modile choix. 66 67 75 use notracer_mod 68 76 … … 75 83 !use deformation_mod ! module concernant les lois de deformation 76 84 use deformation_mod_2lois ! module concernant les lois de deformation 77 78 !--------------Module propriete thermique de la glace79 use prop_thermiques_ice80 ! use prop_therm_ice_heino81 85 82 86 ! eau basale … … 104 108 ! use dragging_beta_buoy ! beta fonction de la flottaison 105 109 106 use dragging_beta_nolin ! beta non lineaire 110 !use dragging_beta_nolin ! beta non lineaire 111 112 use dragging_neff_slope ! dragging Aurelien et Christophe 107 113 108 114 !use beta_iter_vitbil_mod … … 111 117 !use dragging_plastic_LGM 112 118 !------------ spinup ----------------------------------- 113 119 use no_spinup ! spinup=0 ou spinup=1 114 120 ! use spinup_vitbil ! plusieurs variantes de spinup : inclue le dragging 115 121 ! enlever le use dragging dans ce cas là … … 133 139 !use bmelt_ant_regions ! pour l'Antarctique avec régions 134 140 135 use bmelt_seuil_prof 141 !use bmelt_seuil_prof 142 use bmelt_ant_regions ! pour l'Antarctique avec régions 136 143 137 ! USEBMELT_NOR_REGIONS ! pour le nord avec régions144 ! use BMELT_NOR_REGIONS ! pour le nord avec régions 138 145 ! use BMELT_NOR_DEPTH ! pour le nord avec profondeur d'eau 139 146 -
trunk/SOURCES/ANT15-LBq_files/paradim-ant15_LBq_mod.f90
r4 r68 16 16 ! character 17 17 character(len=7), parameter :: geoplace='A-LBq15' 18 character(len=2 3), parameter :: dirnameinp='../INPUT/ANT-LBq-15km/' !< input directory19 character(len= 17), parameter :: dirforcage='../INPUT/Forcage/' !< input directory18 character(len=25), parameter :: dirnameinp='../../INPUT/ANT-LBq-15km/' !< input directory 19 character(len=20), parameter :: dirforcage='../../INPUT/Forcage/' !< input directory 20 20 21 21 ! dimensionnement grilles -
trunk/SOURCES/Ant40_files/bmelt-ant-regions_mod.f90
r49 r68 14 14 MODULE BMELT_ANT_REGIONS ! cat juillet 2005 15 15 16 USE module3D_phy 17 16 use module3D_phy 18 17 19 18 implicit none … … 26 25 !< utilises pour moduler la fusion sous le shelf 27 26 integer, dimension(10) :: region !< pour écrire dans le fichier param 28 character(len=30),dimension(10) :: regname !< nom des régions 27 character(len=30),dimension(10) :: regname !< nom des régions 28 character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves 29 29 real :: bsupshelf 30 30 integer :: grd !< pour une sortie … … 51 51 bmelt_Amery,bmgrz_Amery,bmelt_PIG,bmgrz_PIG, & 52 52 bmelt_Pen,bmgrz_Pen,bmelt_other,bmgrz_other, & 53 bmelt_talus,bmgrz_talus,bmelt_coef 53 bmelt_talus,bmgrz_talus,bmelt_coef,file_number_shelves 54 54 55 55 … … 78 78 write(num_rep_42,*) 'bmgrz_talus = ',bmgrz_talus 79 79 write(num_rep_42,*) 'bmelt_coef = ',bmelt_coef 80 write(num_rep_42,*) 'file_numer_shelves = ',file_number_shelves 80 81 write(num_rep_42,*)'/' 81 82 write(num_rep_42,428) '! bmelt_Ross & bmgrz_Ross : fonte basale secteur Ross' … … 87 88 write(num_rep_42,428) '! bmelt_talus & bmgrz_talus : fonte basale apres talus cont' 88 89 write(num_rep_42,428) '! bmelt_coef : coef fonte (1 pour conserver val)' 90 write(num_rep_42,428) '! file_numer_ice-shelves : fichier zones ice shelves' 89 91 write(num_rep_42,*) 90 92 91 93 92 ! ecriture dans le fichier parametres93 write(num_rep_42,*)'fusion basale sous les ice shelves : bmelt-ant-regions_mod'94 write(num_rep_42,*)'----------------------------------------------------------'95 96 97 ! lecture du fichier contenant les distances au talu continental (m)98 open(88,file=TRIM(DIRNAMEINP)//'distance_talu-40km.xy')99 100 do j=1,ny101 do i=1,nx102 read(88,'(i3,1x,i3,1x,f10.2)') k,k,dist_talu(i,j)103 enddo104 enddo105 close(88)106 107 ! les bords de la grille sont mis a dist_talus=0 de force108 109 110 dist_talu(1,:)=0. !a gauche111 dist_talu(2,:)=0.112 113 dist_talu(nx-1,:)=0. ! a droite114 dist_talu(nx,:)=0.115 116 dist_talu(:,1)=0. ! en bas117 dist_talu(:,2)=0.118 119 dist_talu(:,ny-1)=0. ! en haut120 dist_talu(:,ny)=0.121 122 123 debug_3D(:,:,32)=dist_talu(:,:)124 94 125 95 ! lecture du fichier contenant les types de shelves Ronne-Flichner 126 96 ! ->1, Ross -> 2 , Amery -> 3, PIG-> 4, les petits shelves au dessus de PIG -> 5 127 128 open(88,file=TRIM(DIRNAMEINP)//'numer-ice-shelves-juil07.dat',status='OLD')129 97 typeshelf(:,:)=100 130 do k=1,nx*ny 131 read(88,*,end=36) i,j,typeshelf(i,j) 132 end do 133 36 close(88) 134 135 98 file_number_shelves=TRIM(DIRNAMEINP)//trim(file_number_shelves) 99 ! call lect_input(1,'z',1,typeshelf,file_number_shelves,file_ncdf) 100 ! open(88,file=TRIM(DIRNAMEINP)//'numer-ice-shelves-juil07.dat',status='OLD') 101 open(88,file=TRIM(file_number_shelves),status='OLD') 102 do j=1,ny 103 do i=1,nx 104 read(88,*) typeshelf(i,j) 105 enddo 106 enddo 107 close(88) 136 108 137 109 region(:)=0 … … 147 119 do i=1,nx 148 120 149 talus: if ( dist_talu(i,j).GT.5000.) then121 talus: if (Bsoc(i,j).LT.2500.) then 150 122 if (nint(typeshelf(i,j)).eq.1) then ! Ronne-Filchner FRis 151 123 bmshelf(i,j)=bmelt_FRis -
trunk/SOURCES/Ant40_files/module_choix-antar40-0.4.f90
r62 r68 47 47 !use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 48 48 !use climat_perturb_mod ! climat perturbe a reverifier Dec 2015 49 use climat_Grice2sea_years_mod ! climat force par fichier SMB directement (grice2sea)50 !use climat_Grice2sea_years_perturb_mod ! climat force par fichier SMB directement (grice2sea) + index temperature carotte de glace49 !use climat_Grice2sea_years_mod ! climat force par fichier SMB directement (grice2sea) 50 use climat_Grice2sea_years_perturb_mod ! climat force par fichier SMB directement (grice2sea) + index temperature carotte de glace 51 51 52 52 … … 80 80 ! sliding- dragging 81 81 ! use module sliding_vitbal 82 use sliding_Bindschadler 82 ! use sliding_Bindschadler 83 use no_sliding 83 84 ! use sliding_dragging_heino ! loi de glissement heino : mis dans diffusiv : 84 85 -
trunk/SOURCES/Antarctique_general_files/output_anta_mod-0.4.f90
r4 r68 123 123 124 124 filout = runname//'_vol_regions.dat' 125 filout = TRIM(DIRNAMEOUT)// 'time-series/'//TRIM(filout)125 filout = TRIM(DIRNAMEOUT)//TRIM(filout) 126 126 open(nvol,file=filout) 127 127 128 128 filout = runname// '_volbuoy_regions.dat' 129 filout = TRIM(DIRNAMEOUT)// 'time-series/'//TRIM(filout)129 filout = TRIM(DIRNAMEOUT)//TRIM(filout) 130 130 open(nvolbuoy,file=filout) 131 131 132 132 filout = runname//'_meanhdot_regions.dat' 133 filout = TRIM(DIRNAMEOUT)// 'time-series/'//TRIM(filout)133 filout = TRIM(DIRNAMEOUT)//TRIM(filout) 134 134 open(nhdot,file=filout) 135 135 136 136 filout = runname//'_sigmahdot_regions.dat' 137 filout = TRIM(DIRNAMEOUT)// 'time-series/'//TRIM(filout)137 filout = TRIM(DIRNAMEOUT)//TRIM(filout) 138 138 open(nsigma_hdot,file=filout) 139 139 140 140 filout = runname//'_nbpoints_regions.dat' 141 filout = TRIM(DIRNAMEOUT)// 'time-series/'//TRIM(filout)141 filout = TRIM(DIRNAMEOUT)//TRIM(filout) 142 142 open(npoints,file=filout) 143 143 -
trunk/SOURCES/Fichiers-parametres/A-LBq15_LISTE-VAR-NETCDF.dat
r4 r68 82 82 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ fin du header commentaire 83 83 ============================================ classe 1 84 1 100 4 ! class 1 nombre de snapmax 100 4 type de noeud84 1 220 4 ! class 1 nombre de snapmax 100 4 type de noeud 85 85 'o' '>' '^' 'x' 86 86 ----------------------------------------------------------- … … 114 114 Bm 115 115 10 116 1 0 -1116 1 1 1 117 117 ----------------------------------------------------------- 118 118 Acc 119 119 11 120 0 4 0120 1 1 1 121 121 ----------------------------------------------------------- 122 122 Abl 123 123 12 124 1 0 -1124 1 1 1 125 125 ----------------------------------------------------------- 126 126 Calv 127 127 13 128 0 1 -1128 0 1 1 129 129 ----------------------------------------------------------- 130 130 dHdt … … 134 134 Tann 135 135 18 136 1 0 -1136 1 1 1 137 137 ----------------------------------------------------------- 138 138 Tjja 139 139 19 140 1 0 -1140 1 1 1 141 141 ----------------------------------------------------------- 142 142 Tb … … 144 144 1 1 1 145 145 ----------------------------------------------------------- 146 ghf 147 23 148 1 1 -1 149 ----------------------------------------------------------- 146 150 phid 147 151 24 148 1 0 -1152 1 1 1 149 153 ----------------------------------------------------------- 150 154 Bmelt 151 155 25 152 1 0 -1156 1 1 1 153 157 ----------------------------------------------------------- 154 158 Uxbar_maj 155 159 30 156 01 1160 1 1 1 157 161 ----------------------------------------------------------- 158 162 Uybar_maj 159 163 31 160 01 1164 1 1 1 161 165 ----------------------------------------------------------- 162 166 Umag … … 168 172 1 1 1 169 173 ----------------------------------------------------------- 170 171 174 bm_ref_t 175 197 176 1 1 1 177 ----------------------------------------------------------- 178 grad_bm 179 198 180 1 1 1 181 ----------------------------------------------------------- 182 bm_time 183 199 184 1 1 1 185 ----------------------------------------------------------- 172 186 ============================================ classe 2 173 2 100 3 ! class 2 nombre de snapmax 100 3 type de noeud a sortir187 2 220 3 ! class 2 nombre de snapmax 100 3 type de noeud a sortir 174 188 '>' '^' 'o' 175 189 ----------------------------------------------------------- … … 204 218 frtmx 205 219 40 206 0 0 -1220 1 1 1 207 221 ----------------------------------------------------------- 208 222 frtmy ! sur noeud mineur 209 223 41 210 0 0 -1224 1 1 1 211 225 ----------------------------------------------------------- 212 226 tobmx ! sur noeud mineur 213 227 42 214 01 1228 1 1 1 215 229 ----------------------------------------------------------- 216 230 tobmy 217 231 43 218 01 1232 1 1 1 219 233 ----------------------------------------------------------- 220 234 toshe … … 224 238 epsxx !'o' 225 239 45 226 1 0-1240 1 3 -1 227 241 ----------------------------------------------------------- 228 242 epsyy !'o' 229 243 46 230 1 0-1244 1 3 -1 231 245 ----------------------------------------------------------- 232 246 epsxy !'x' 233 247 47 234 1 01248 1 3 1 235 249 ----------------------------------------------------------- 236 250 eps 'o' 237 251 48 238 1 0 1 239 ----------------------------------------------------------- 240 Abar 'o' 241 49 242 1 1 1 252 1 3 1 243 253 ----------------------------------------------------------- 244 254 pvi ! 'o' … … 264 274 Hwat 265 275 60 266 01 1276 1 1 1 267 277 ----------------------------------------------------------- 268 278 Hdw … … 292 302 Nefmx 293 303 68 294 01 1304 1 1 1 295 305 ----------------------------------------------------------- 296 306 Nefmy 297 307 69 298 02 1308 1 2 1 299 309 ----------------------------------------------------------- 300 310 posx … … 316 326 H_uw_x 317 327 192 318 01 1328 1 1 1 319 329 ----------------------------------------------------------- 320 330 H_uw_y 321 331 193 322 0 1 1323 -----------------------------------------------------------324 time_float325 194326 332 1 1 1 327 333 ----------------------------------------------------------- -
trunk/SOURCES/Fichiers-parametres/A-LBq15_TEMPS-NETCDF.dat
r4 r68 15 15 ----------------------------------------------------------------------- 16 16 4 ndtsortie lecture de dtncdf de sortie 17 10.17 5000. 18 18 10000. 19 19 20000. 20 20 1.e10 21 21 ------------------------------------------------------------------------------------------- 22 24 npredeft lecture de predef_tsort 23 1990 24 1991 25 1992 26 1993 27 1994 28 1995 29 1996 30 1997 31 1998 32 1999 22 9 npredeft lecture de predef_tsort 23 -21000 24 10 25 100 26 200 27 500 28 1000 33 29 2000 34 200135 200236 200337 200438 200539 200640 200741 200842 200843 201044 201145 201246 201547 248 649 750 851 952 10.53 1154 1255 1356 1457 1658 1759 1860 1961 15.62 20.63 30.64 40.65 50.66 15067 20.68 100.69 500.70 5.71 10.72 =======73 20.74 50.75 100.76 200.77 500.78 1000.79 10000.80 20000.81 70.82 80.83 90.84 100085 1000086 5000087 2588 3589 6690 10091 15092 20093 25094 30095 35096 40097 45098 50099 300100 400101 3000102 4000103 30 5000 104 5010105 5020106 5050107 5100108 5200109 5300110 5400111 5500112 5600113 5700114 5800115 5900116 6000117 6500118 31 7000 119 7500120 8000121 122 32 ------------------------------------------------------------------------------------------- 123 -
trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90
r67 r68 421 421 end if 422 422 423 coefbmshelf=(1.+tafor/10.) 423 ! coefmshelf est un coefficient qui fait varier bmgrz et bmshelf en fonction de tafor 424 coefbmshelf=(1.+tafor/10.) ! coefbmshelf=0 pour tafor=-10deg 424 425 coefbmshelf=max(coefbmshelf,0.) 425 426 coefbmshelf=min(coefbmshelf,2.) -
trunk/SOURCES/Makefile.grisli.inc
r65 r68 172 172 173 173 174 Liste_ANT15-LBq = lect-Ant_clim_Acc-T_gen_dat.o output_anta_mod-0.4.o \ 175 dragging-vit_bil_LBq_gen_mod.o dragging_LGM_mod.o dragging_plastic_LGM_mod.o \ 176 dragging_prescr_beta_mod.o dragging_prescr_beta_buoyency_mod.o dragging_prescr_beta_nolin_mod.o \ 177 fake-routines-ant_mod.o beta_iter_vitbil_mod.o \ 178 lect-Ant_gen2010_dat.o module_choix_antar15_LBq.o \ 179 massb-ant_perturb_Tparam.o track_ant40_mod.o 180 181 174 #Liste_ANT15-LBq = lect-Ant_clim_Acc-T_gen_dat.o output_anta_mod-0.4.o \ 175 # dragging-vit_bil_LBq_gen_mod.o dragging_LGM_mod.o dragging_plastic_LGM_mod.o \ 176 # dragging_prescr_beta_mod.o dragging_prescr_beta_buoyency_mod.o dragging_prescr_beta_nolin_mod.o \ 177 # fake-routines-ant_mod.o beta_iter_vitbil_mod.o \ 178 # lect-Ant_gen2010_dat.o module_choix_antar15_LBq.o \ 179 # massb-ant_perturb_Tparam.o track_ant40_mod.o 180 181 Liste_ANT15-LBq = output_anta_mod-0.4.o \ 182 lect-Ant_gen2010_dat.o \ 183 lect-anteis_mod.o \ 184 bmelt-ant-regions_mod.o \ 185 fake-routines-ant_mod.o \ 186 module_choix_antar15_LBq.o \ 187 track_ant40_mod.o \ 182 188 183 189 Liste_hudson = sedim_declar_hudson_mod.o climat-hudson_mod.o \ … … 361 367 $(FT) $(NETCDFINCLUDE) -c Ant40_files/$*.f90 362 368 369 # ANT15-LBq_files 370 %.o: ANT15-LBq_files/%.f90 371 $(FT) $(NETCDFINCLUDE) -c ANT15-LBq_files/$*.f90 372 363 373 # Greeneem_files/Greeneem15_files 364 374 %.o: Greeneem_files/Greeneem15_files/%.f90 … … 410 420 #--------------------------------------- 411 421 412 LBqAntar15 : $(Dim_ANT15-LBq) $(mod_dim_communs) \ 413 $(Proto_recul) \ 414 $(mod_communs) \ 415 $(Proto_recul) \ 416 $(mod_clim_perturb) \ 417 $(mod_tracers) \ 422 Ant-15 : $(Dim_ANT15-LBq) $(mod_dim_communs) \ 423 $(toy_recul) \ 424 $(mod_communs) \ 425 $(mod_clim_tof) \ 426 $(mod_no_tracers) \ 418 427 $(mod_ell) $(Liste_ANT15-LBq) \ 419 428 $(diagnoshelf) \ 420 429 $(Liste_Netcdf) \ 421 $(routines_communes) steps_time_loop.o $(routine_elliptiques) 422 423 $(LK) -o ../bin/LBqAntar-15 $(Dim_ANT15-LBq) $(mod_dim_communs) \ 424 $(Proto_recul) \ 425 $(mod_communs) \ 426 $(mod_clim_perturb) \ 427 $(mod_tracers) \ 428 $(mod_ell) \ 429 $(Liste_ANT15-LBq) \ 430 $(diagnoshelf) $(Liste_Netcdf) \ 431 $(routines_communes) steps_time_loop.o $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) 430 $(routines_communes) steps_time_loop.o \ 431 $(routine_elliptiques) \ 432 $(Liste_BLAS) 433 434 $(LK) -o ../bin/Ant-15 \ 435 $(Dim_ANT15-LBq) $(mod_dim_communs) \ 436 $(toy_recul) \ 437 $(mod_communs) \ 438 $(mod_clim_tof) \ 439 $(mod_no_tracers) \ 440 $(mod_ell) $(Liste_ANT15-LBq) \ 441 $(diagnoshelf) \ 442 $(Liste_Netcdf) \ 443 $(routines_communes) steps_time_loop.o \ 444 $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) $(Liste_BLAS) 432 445 433 446 Recul_ice2sea : $(Dim_ANT15-LBq) $(mod_dim_communs) \ -
trunk/SOURCES/diffusiv-polyn-0.6.f90
r65 r68 66 66 call slope_surf 67 67 68 !call sliding ! au sens vitesse de glissement, afq, commented 21 apr 16 68 call sliding ! au sens vitesse de glissement 69 69 70 70 ! le glissement est maintenant dans un module a part choisi dans le module choix -
trunk/SOURCES/initial-0.3.f90
r65 r68 120 120 call init_diagno ! initialisation de la resolution equation elliptique vitesses 121 121 122 !call init_sliding ! initialisation du glissement, afq, commented 21 apr 16122 call init_sliding ! initialisation du glissement 123 123 124 124 call init_spinup ! initialisation du spinup
Note: See TracChangeset
for help on using the changeset viewer.