- Timestamp:
- 01/31/24 16:56:43 (5 months ago)
- Location:
- branches/GRISLIv3
- Files:
-
- 6 deleted
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/GRISLIv3/SOURCES/3D-physique-gen_mod.f90
r470 r474 127 127 ! dimensionnements qui etaient avant dans H_prescribe 128 128 real,dimension(nx,ny) :: Hp !< H value if prescribed 129 real,dimension(nx,ny) :: Hp0 !< H value if prescribed (reference value)130 129 integer,dimension(nx,ny) :: i_Hp !< 1 if H is prescribed on this node, else 0 131 integer,dimension(nx,ny) :: i_Hp0 !< i_hp mask reference value does not change with time132 130 integer, dimension(nx,ny) :: imx_diag !< masque pour eq elliptique 133 131 integer, dimension(nx,ny) :: imy_diag !< masque pour eq elliptique … … 228 226 logical,dimension(nx,ny) :: FLOTMX !< if flotting (one neighboord floats) '>' 229 227 logical,dimension(nx,ny) :: FLOTMY !< if flotting (one neighboord floats) '^' 230 logical,dimension(nx,ny) :: FRONT_RESC!< comme front (mais boolean) pour remplimat rescue231 228 logical,dimension(nx,ny) :: FBM !< 232 229 logical,dimension(nx,ny) :: GZMX !< point stream ">" … … 248 245 integer :: num_ritz = 35 !< Id of ritz 249 246 integer :: num_rep_42 = 42 !< Id of reponse 250 integer :: num_trace = 81 !< Id of trace-reboot251 247 integer :: num_file = 88 !< Id of 252 248 integer :: num_kelvin = 90 !< Id of kelvin.res -
branches/GRISLIv3/SOURCES/Alps_files/lect-alps_mod.f90
r356 r474 12 12 character(len=100) :: ghf_fich ! fichier grille 13 13 character(len=80) :: filin 14 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure15 14 character(len=100) :: file_ncdf !< fichier netcdf issue des fichiers .dat 16 15 -
branches/GRISLIv3/SOURCES/Alps_files/module_choix-alps.f90
r450 r474 23 23 !! @note - use bmelt_nor_regions 24 24 !! @note - use fake_nor 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_hemin40_mod 29 26 !! @note - use printtable … … 133 130 !-------------- Outputs----------------------------- 134 131 ! 135 use out_profile !profils136 !use out_cptr ! fichier de reprise137 !use out_hz ! module de sorties snapshot dans un plan horizontal138 139 132 !-----Suivi temporel 140 133 !use output_hemin40_mod -
branches/GRISLIv3/SOURCES/Ant16_files/module_choix-ant16km.f90
r452 r474 23 23 !! @note - use bmelt_ant_regions 24 24 !! @note - use fake_ant 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_antarcti_mod 29 26 !! @note - use printtable … … 126 123 !-------------- Outputs----------------------------- 127 124 ! 128 use out_profile !profils129 130 125 !-----Suivi temporel 131 126 ! use output_heminord_mod -
branches/GRISLIv3/SOURCES/Ant40_files/lect-anteis_mod.f90
r465 r474 33 33 ! character(len=100) :: file2 34 34 35 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure36 35 ! afq 8/1/19: real :: sealevel0 37 36 … … 69 68 70 69 71 !!$! lecture de la topo actuelle72 !!$! ---------------------------73 !!$ open (20,file=TRIM(DIRNAMEINP)//file2,status='old')74 !!$75 !!$ read(20,'(A80)') TITRE76 !!$ read(20,*) NI,NJ,NXX,NYY,STEP77 !!$ read(20,*)78 !!$ do J=1,ny79 !!$ do I=1,nx80 !!$ read (20,*) S0(I,J),H0(I,J),BSOC0(I,J)81 !!$ S0(i,j)=max(S0(i,j),0.)82 !!$ end do83 !!$ end do84 !!$ close(20)85 !!$86 !!$87 !!$! lecture de la topo de depart88 !!$! ---------------------------89 !!$ open (20,file=TRIM(DIRNAMEINP)//file1,status='old')90 !!$! open (20,file='../INPUT-DATA/hemin.g50')91 !!$ read(20,'(A80)') TITRE92 !!$ read(20,*) NI,NJ,NXX,NYY,STEP93 !!$ read(20,*)94 !!$ do J=1,ny95 !!$ do I=1,nx96 !!$ read (20,*) S(I,J),H(I,J),BSOC(I,J)97 !!$ end do98 !!$ end do99 !!$ close(20)100 101 ! lecture de la topo de référence102 ! -------------------------------103 ! Cette topo sert a calculer le socle de reference pour l'isostasie104 ! voir init_iso et a avoir une surface de reference pour les temperatures105 ! lecture adaptee aux fichiers ZBL.dat ou netcdf ou grd106 ! call lect_input(1,'Bsoc',1,Bsoc0,topo_ref,file_ncdf) ! socle107 ! call lect_input(1,'S',1,S0,topo_ref,file_ncdf) ! surface108 ! call lect_input(1,'H',1,H0,topo_ref,file_ncdf) ! epaisseur109 110 70 ! lecture pour eviter plantage avec compile -O0 111 71 call Read_Ncdf_var('Bsoc',topo_ref,tab) … … 115 75 call Read_Ncdf_var('H',topo_ref,tab) 116 76 H0(:,:) = tab(:,:) 117 !~ !cdc correction de 3 pts qui posent probleme (pente tres forte):118 !~ S0(38,53)=1400.119 !~ Bsoc0(38,53)=-1000.120 !~ H0(38,53)=S0(38,53)-Bsoc0(38,53)121 !~122 !~ S0(35,53)=1300.123 !~ Bsoc0(35,53)=-1000.124 !~ H0(35,53)=S0(35,53)-Bsoc0(35,53)125 !~126 !~ S0(35,56)=1200.127 !~ Bsoc0(35,56)=-1000.128 !~ H0(35,56)=S0(35,56)-Bsoc0(35,56)129 130 !cdc correction point pole sud :131 ! S0(71,71)=(S0(71,70)+S0(71,72)+S0(70,71)+S0(72,71))/4.132 ! Bsoc0(71,71)=(Bsoc0(71,70)+Bsoc0(71,72)+Bsoc0(70,71)+Bsoc0(72,71))/4.133 134 77 135 78 ! where (S0(:,:).GT.0) … … 150 93 ! lecture de la topo de depart 151 94 ! --------------------------- 152 ! lecture adaptee aux fichiers ZBL.dat ou netcdf ou grd153 ! call lect_input(1,'Bsoc',1,Bsoc,topo_dep,file_ncdf) ! socle154 ! call lect_input(1,'S',1,S,topo_dep,file_ncdf) ! surface155 ! call lect_input(1,'H',1,H,topo_dep,file_ncdf) ! epaisseur156 95 ! lecture pour eviter plantage avec compile -O0 157 96 call Read_Ncdf_var('Bsoc',topo_dep,tab) … … 162 101 H(:,:) = tab(:,:) 163 102 164 !~ !cdc correction de 3 pts qui posent probleme (pente tres forte):165 !~ S(39,54)=1400.166 !~ Bsoc(39,54)=-1000.167 !~ H(39,54)=S(39,54)-Bsoc(39,54)168 !~169 !~ S(36,54)=1300.170 !~ Bsoc(36,54)=-1000.171 !~ H(36,54)=S(36,54)-Bsoc(36,54)172 !~173 !~ S(36,57)=1200.174 !~ Bsoc(36,57)=-1000.175 !~ H(36,57)=S(36,57)-Bsoc(36,57)176 177 ! S(71,71)=(S(71,70)+S(71,72)+S(70,71)+S(72,71))/4.178 ! Bsoc(71,71)=(Bsoc(71,70)+Bsoc(71,72)+Bsoc(70,71)+Bsoc(72,71))/4.179 103 ! where (S(:,:).GT.0) 180 104 ! H(:,:)=S(:,:)-BSOC(:,:) … … 184 108 ! afq 8/1/19: S(:,:)=max(S(:,:),0.) ! pour etre au niveau des mers : ATTENTION si SEALEV <0 185 109 H(:,:)=max(H(:,:),0.) ! pas d epaisseurs negatives 186 187 188 !!$189 !!$ ! socle190 !!$ filin='bedelev-2000-40km.dat'191 !!$ call lect_eis(nx,ny,BSOC,filin,DIRNAMEINP)192 !!$ write(num_rep_42,*) 'fichier BSOC (socle) : ', filin193 !!$ Bsoc0(:,:)=Bsoc(:,:)194 !!$ ! surface195 !!$ filin='surface-2000-40km.dat'196 !!$ call lect_eis(nx,ny,S,filin,DIRNAMEINP)197 !!$ write(num_rep_42,*) 'fichier surface : ', filin198 !!$199 !!$ ! epaisseur200 !!$ filin='icethic-2000-40km.dat'201 !!$ call lect_eis(nx,ny,H,filin,DIRNAMEINP)202 !!$ write(num_rep_42,*) 'fichier epaisseur : ', filin203 !!$ H0(:,:)=H(:,:)204 !!$205 !!$206 !!$ ! masque207 !!$ filin='maskHUY40km.dat'208 !!$ call lect_ieis(nx,ny,MK0,filin,DIRNAMEINP)209 !!$ write(num_rep_42,*) 'fichier masque : ', filin210 !!$211 !!$ ! enlever les epaisseurs fictives de Philippe212 !!$ do I=1,NX213 !!$ do J=1,NY214 !!$ ! attention les valeurs de masque de Philippe sont differentes215 !!$ if ((MK0(I,J).gt.1).and.(H(i,j).le.218.00)) then ! mer libre216 !!$ H(i,j)=1.217 !!$ S(I,J)=H(i,j)-(RO/ROW)*H(i,j)218 !!$ endif219 !!$ end do220 !!$ end do221 222 110 223 111 ! S0(:,:)=S(:,:) … … 236 124 enddo 237 125 enddo 238 239 240 ! calcul des courbures du socle241 242 ! call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &243 ! bidon(:,:,4),socle_cry,bidon(:,:,5))244 ! socle_cry(:,:)=socle_cry(:,:)*dx*dx245 246 ! lecture des coordonnées geographiques247 ! filin=TRIM(DIRNAMEINP)//'coord-Ant-40km.dat'248 126 249 127 ! les coordonnees sont calculees en °dec avec GMT, … … 264 142 open(88,file=ghf_fich) 265 143 266 ! write(num_rep_42,*) 'flux geothermique Shapiro : ', TRIM(DIRNAMEINP)//'ijphi-40km-ant.dat'267 ! do k=1,nx*ny268 ! read(88,*) i,j,ghf(i,j)269 ! end do270 ! close(88)271 272 ! call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf)273 144 ! pour eviter plantage -O0 274 145 call Read_Ncdf_var('ghf',ghf_fich,tab) -
branches/GRISLIv3/SOURCES/Ant40_files/module_choix-antar40-0.4.f90
r465 r474 23 23 !! @note - use bmelt_ant_regions 24 24 !! @note - use fake_ant 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_antarcti_mod 29 26 !< … … 125 122 !-------------- Outputs----------------------------- 126 123 ! 127 use out_profile !profils128 129 124 !-----Suivi temporel 130 125 ! use output_heminord_mod -
branches/GRISLIv3/SOURCES/Eurasie_files/lect-eurasie_mod.f90
r297 r474 12 12 character(len=100) :: ghf_fich ! fichier grille 13 13 character(len=80) :: filin 14 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure15 14 character(len=100) :: file_ncdf !< fichier netcdf issue des fichiers .dat 16 15 … … 117 116 H(:,:)=varloc(:,:) 118 117 119 120 ! calcul des courbures du socle121 ! call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &122 ! bidon(:,:,4),socle_cry,bidon(:,:,5))123 ! socle_cry(:,:)=socle_cry(:,:)*dx*dx124 125 ! lecture des coordonnées geographiques126 127 ! filin=TRIM(DIRNAMEINP)//grid_topo128 129 118 ! les coordonnees sont calculees en °dec avec GMT, 130 119 ! les longitudes sont comprises entre -180 et +180 (negative a l'Ouest de … … 140 129 xmax=xcc(nx,ny)/1000. 141 130 ymax=ycc(nx,ny)/1000. 142 143 ! lecture du flux geothermique de Shapiro144 ! open(88,file=TRIM(DIRNAMEINP)//'ijphi_hemin40.dat')145 ! open(88,file=ghf_fich)146 131 147 !!$ write(42,*) 'flux geothermique Shapiro : ',TRIM(DIRNAMEINP)//'ijphi_hemin40.dat'148 !!$ do k=1,nx*ny149 !!$ read(88,*) i,j,ghf(i,j)150 !!$ end do151 !!$ close(88)152 132 153 133 call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf) … … 158 138 ! ghf(:,:)=-SECYEAR/1000.*55. !B6norcg2 159 139 160 ! print*,'lect topo'161 ! print*,'shb',S(101,91),H(101,91),B(101,91)162 ! print*,'shb0',S0(101,91),H0(101,91),BSOC0(101,91)163 140 ! Initialisation du Masque 164 141 !------------------------------------------------ -
branches/GRISLIv3/SOURCES/Eurasie_files/module_choix-eurasie.f90
r450 r474 23 23 !! @note - use bmelt_nor_regions 24 24 !! @note - use fake_nor 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_hemin40_mod 29 26 !! @note - use printtable … … 131 128 !-------------- Outputs----------------------------- 132 129 ! 133 use out_profile !profils134 !use out_cptr ! fichier de reprise135 !use out_hz ! module de sorties snapshot dans un plan horizontal136 137 130 !-----Suivi temporel 138 131 !use output_hemin40_mod -
branches/GRISLIv3/SOURCES/Grismip6_files/module_choix-grismip6.f90
r452 r474 23 23 !! @note - use bmelt_ant_regions 24 24 !! @note - use fake_ant 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_antarcti_mod 29 26 !! @note - use printtable … … 127 124 !-------------- Outputs----------------------------- 128 125 ! 129 use out_profile !profils130 131 126 !-----Suivi temporel 132 127 ! use output_heminord_mod -
branches/GRISLIv3/SOURCES/Hemin15_files/lect-hemin15_mod.f90
r224 r474 11 11 character(len=100) :: grid_topo ! fichier grille 12 12 character(len=100) :: ghf_fich ! fichier grille 13 character(len=80) :: filin 14 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure 13 character(len=80) :: fili 15 14 character(len=100) :: file_ncdf !< fichier netcdf issue des fichiers .dat 16 15 … … 55 54 56 55 57 ! calcul des courbures du socle58 ! call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &59 ! bidon(:,:,4),socle_cry,bidon(:,:,5))60 ! socle_cry(:,:)=socle_cry(:,:)*dx*dx61 62 56 ! lecture des coordonnées geographiques 63 57 ! les coordonnees sont calculees en °dec avec GMT, -
branches/GRISLIv3/SOURCES/Hemin15_files/module_choix-hemin15.f90
r450 r474 23 23 !! @note - use equat_adv_diff_2D 24 24 !! @note - use fake_nor 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_hemin40_mod 29 26 !! @note - use printtable … … 121 118 !-------------- Outputs----------------------------- 122 119 ! 123 use out_profile !profils124 125 120 !-----Suivi temporel 126 121 use output_hemin40_mod -
branches/GRISLIv3/SOURCES/Hemin40_files/lect-hemin40_mod.f90
r465 r474 14 14 character(len=100) :: ghf_fich ! fichier grille 15 15 character(len=80) :: filin 16 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure17 16 character(len=100) :: file_ncdf !< fichier netcdf issue des fichiers .dat 18 17 … … 120 119 121 120 122 ! calcul des courbures du socle123 ! call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &124 ! bidon(:,:,4),socle_cry,bidon(:,:,5))125 ! socle_cry(:,:)=socle_cry(:,:)*dx*dx126 127 ! lecture des coordonnées geographiques128 129 ! filin=TRIM(DIRNAMEINP)//grid_topo130 131 121 ! les coordonnees sont calculees en °dec avec GMT, 132 122 ! les longitudes sont comprises entre -180 et +180 (negative a l'Ouest de … … 142 132 xmax=xcc(nx,ny)/1000. 143 133 ymax=ycc(nx,ny)/1000. 144 145 ! lecture du flux geothermique de Shapiro146 ! open(88,file=TRIM(DIRNAMEINP)//'ijphi_hemin40.dat')147 ! open(88,file=ghf_fich)148 134 149 !!$ write(42,*) 'flux geothermique Shapiro : ',TRIM(DIRNAMEINP)//'ijphi_hemin40.dat'150 !!$ do k=1,nx*ny151 !!$ read(88,*) i,j,ghf(i,j)152 !!$ end do153 !!$ close(88)154 135 155 136 call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf) … … 160 141 ! ghf(:,:)=-SECYEAR/1000.*55. !B6norcg2 161 142 162 ! print*,'lect topo'163 ! print*,'shb',S(101,91),H(101,91),B(101,91)164 ! print*,'shb0',S0(101,91),H0(101,91),BSOC0(101,91)165 143 ! Initialisation du Masque 166 144 !------------------------------------------------ -
branches/GRISLIv3/SOURCES/Hemin40_files/module_choix-hemin40-0.4.f90
r473 r474 23 23 !! @note - use equat_adv_diff_2D 24 24 !! @note - use fake_nor 25 !! @note - use out_profile26 !! @note - use out_cptr27 !! @note - use out_hz28 25 !! @note - use output_hemin40_mod 29 26 !! @note - use printtable … … 130 127 !-------------- Outputs----------------------------- 131 128 ! 132 use out_profile !profils133 !use out_cptr ! fichier de reprise134 !use out_hz ! module de sorties snapshot dans un plan horizontal135 136 129 !-----Suivi temporel 137 130 use output_hemin40_mod -
branches/GRISLIv3/SOURCES/Laure16_files/lect-laurentide_mod.f90
r465 r474 12 12 character(len=100) :: ghf_fich ! fichier grille 13 13 character(len=80) :: filin 14 real, dimension(nx,ny,5) :: bidon ! pour l'appel a courbure15 14 character(len=100) :: file_ncdf !< fichier netcdf issue des fichiers .dat 16 15 … … 120 119 H(:,:)=varloc(:,:) 121 120 122 123 ! calcul des courbures du socle124 ! call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &125 ! bidon(:,:,4),socle_cry,bidon(:,:,5))126 ! socle_cry(:,:)=socle_cry(:,:)*dx*dx127 128 ! lecture des coordonnées geographiques129 130 ! filin=TRIM(DIRNAMEINP)//grid_topo131 132 121 ! les coordonnees sont calculees en °dec avec GMT, 133 122 ! les longitudes sont comprises entre -180 et +180 (negative a l'Ouest de … … 145 134 xmax=xcc(nx,ny)/1000. 146 135 ymax=ycc(nx,ny)/1000. 147 148 ! lecture du flux geothermique de Shapiro149 ! open(88,file=TRIM(DIRNAMEINP)//'ijphi_hemin40.dat')150 ! open(88,file=ghf_fich)151 136 152 !!$ write(42,*) 'flux geothermique Shapiro : ',TRIM(DIRNAMEINP)//'ijphi_hemin40.dat'153 !!$ do k=1,nx*ny154 !!$ read(88,*) i,j,ghf(i,j)155 !!$ end do156 !!$ close(88)157 137 158 138 call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf) … … 163 143 ! ghf(:,:)=-SECYEAR/1000.*55. !B6norcg2 164 144 165 ! print*,'lect topo'166 ! print*,'shb',S(101,91),H(101,91),B(101,91)167 ! print*,'shb0',S0(101,91),H0(101,91),BSOC0(101,91)168 145 ! Initialisation du Masque 169 146 !------------------------------------------------ -
branches/GRISLIv3/SOURCES/Laure16_files/module_choix-laurentide.f90
r450 r474 24 24 !! @note - use bmelt_nor_regions 25 25 !! @note - use fake_nor 26 !! @note - use out_profile27 !! @note - use out_cptr28 !! @note - use out_hz29 26 !! @note - use output_hemin40_mod 30 27 !! @note - use printtable … … 134 131 !-------------- Outputs----------------------------- 135 132 ! 136 use out_profile !profils137 !use out_cptr ! fichier de reprise138 !use out_hz ! module de sorties snapshot dans un plan horizontal139 140 133 !-----Suivi temporel 141 134 !use output_hemin40_mod -
branches/GRISLIv3/SOURCES/MISMIP3D_files/module_choix_mismip3d.f90
r450 r474 24 24 !! @note - use bmelt_seuil_prof 25 25 !! @note - use fake_ant 26 !! @note - use out_profile27 !! @note - use out_cptr28 !! @note - use out_hz29 26 !! @note - use output_antarcti_mod 30 27 !! @note - use printtable … … 148 145 !-------------- Outputs----------------------------- 149 146 ! 150 use out_profile !profils151 !use out_cptr ! fichier de reprise152 use out_hz ! module de sorties snapshot dans un plan horizontal153 154 147 !-----Suivi temporel 155 148 ! use output_heminord_mod -
branches/GRISLIv3/SOURCES/Makefile.grisli.inc
r453 r474 62 62 interface_input.o \ 63 63 iso_declar_mod-0.3.o isostasie_mod-0.3.o noisostasie_mod-0.3.o \ 64 deformation_mod_2lois.o deformation_mod_2lois_isotherme.o \ 65 no_sliding.o \ 64 deformation_mod_2lois.o no_sliding.o \ 66 65 fake_beta_iter_vitbil_mod.o \ 67 66 dragging_param_beta_mod.o dragging_param_beta_sedim_mod.o \ … … 70 69 calving_frange_ISMIP_glaciers.o calving_frange_ISMIP_fracture.o \ 71 70 no_calving.o no_rsl.o imposed_regions_rsl.o lake_rsl_mod.o \ 72 out_profile_mod.omix-SIA-L1_mod.o \71 mix-SIA-L1_mod.o \ 73 72 furst_schoof_mod.o bilan_flux_output_mod.o \ 74 73 relaxation_water_diffusion.o \ 75 prescribe-H-i2s_mod.o \76 74 resol_adv_diff_2D-sept2009.o \ 77 75 conserv-mass-adv-diff_sept2009_mod.o \ … … 122 120 out_cptr_mod.o lineartemp-0.2.o litho-0.4.o main3D-0.4-40km.o masque.o \ 123 121 moy_mxmy.o neffect-0.4.o next-time.o \ 124 no-icethick.o printdebug.o \125 122 slope_surf.o \ 126 123 strain_rate.o \ 127 124 tab-litho-0.3.o taubed-0.3.o \ 128 125 velocities-polyn-0.3.o \ 129 courbures.osteps_time_loop.o126 steps_time_loop.o 130 127 131 128 -
branches/GRISLIv3/SOURCES/conserv-mass-adv-diff_sept2009_mod.f90
r468 r474 24 24 use runparam, only: itracebug,num_tracebug 25 25 use reso_adv_diff_2D_vect, only: init_reso_adv_diff_2D,resol_adv_diff_2D_vect 26 ! use prescribe_H, only:prescribe_fixed_points,prescribe_paleo_gl_shelf,prescribe_present_H_gl,&27 ! prescribe_present_H_gl_bmelt,init_prescribe_H28 ! use toy_retreat_mod, only:time_step_recul29 26 30 27 implicit none … … 72 69 73 70 call init_reso_adv_diff_2D 74 ! call init_prescribe_H ! initialize grounding line mask75 71 76 72 Dx1=1./Dx -
branches/GRISLIv3/SOURCES/initial-0.3.f90
r446 r474 56 56 !initialisation du niveau marin local 57 57 call input_rsl 58 59 60 ! call input_profile() ! Les fichiers profils61 !62 58 63 59 if (icompteur.ne.0) then ! reprise d un fichier cptr -
branches/GRISLIv3/SOURCES/steps_time_loop.f90
r468 r474 106 106 ! outputs 107 107 !==================================================================== 108 109 ! vertical plan snapshots (profiles)110 !------------------------------------------------------------------111 112 !if ((NT.eq.1) & ! .or.(NT.eq.2)113 ! .or.(mod(abs(dble(TIME)),dble(DTPROFILE)).lt.dble(dtmin)) &114 ! .or.(abs(TIME-TEND).lt.dtmin)) then115 ! call sortieprofile()116 !endif117 108 118 109 ! sorties netcdf Hassine (aout 2010) (2D and 3D) -
branches/GRISLIv3/launching_area/Master-files/module_choix-antar40-0.4-cycles.f90
r473 r474 24 24 !! @note - use bmelt_ant_regions 25 25 !! @note - use fake_ant 26 !! @note - use out_profile27 !! @note - use out_cptr28 !! @note - use out_hz29 26 !! @note - use output_antarcti_mod 30 27 !< … … 127 124 !-------------- Outputs----------------------------- 128 125 ! 129 use out_profile !profils130 131 126 !-----Suivi temporel 132 127 ! use output_heminord_mod
Note: See TracChangeset
for help on using the changeset viewer.