| 1 | ! ********************************************************************** |
|---|
| 2 | ! * GRISLI Grenoble Ice Shelves-Land Ice |
|---|
| 3 | ! ********************************************************************** |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | ! Ont participe a l'ecriture de ce modele : |
|---|
| 7 | ! |
|---|
| 8 | ! Catherine Ritz (tout du long) |
|---|
| 9 | ! Adeline Fabre (la partie Gremlins) |
|---|
| 10 | ! Vincent Rommelaere (ice shelves et ice streams) |
|---|
| 11 | ! Christophe Dumas (debut f90, (Antarctique) |
|---|
| 12 | ! Vincent Peyaud (portage HN,calving, front, hydrol) |
|---|
| 13 | ! Cyril Mazauric (AGRIF) |
|---|
| 14 | ! Hassine Baya (netcdf, doxygen, icetemp,...) |
|---|
| 15 | ! |
|---|
| 16 | ! catritz@lgge.obs.ujf-grenoble.fr |
|---|
| 17 | ! |
|---|
| 18 | ! ********************************************************************** |
|---|
| 19 | |
|---|
| 20 | |
|---|
| 21 | |
|---|
| 22 | !> \mainpage GRISLI Modele 3D De Calotte Glaciaire |
|---|
| 23 | !! |
|---|
| 24 | !! \section start Pour commencer |
|---|
| 25 | !! Le programme principal est dans le module main3D. |
|---|
| 26 | !! Ce module est dans le fichier main3D-0.4-40km.f90. |
|---|
| 27 | !! |
|---|
| 28 | !! \section tree Arbre d'appel |
|---|
| 29 | !! |
|---|
| 30 | !! - call grisli_init() |
|---|
| 31 | !! - step_grisli() |
|---|
| 32 | !! - sortie_ncdf_cat() |
|---|
| 33 | !! - testsort_time_ncdf() |
|---|
| 34 | !! - initial() |
|---|
| 35 | !! - sortie_hz_multi() |
|---|
| 36 | !! - call step_grisli1() |
|---|
| 37 | !! |
|---|
| 38 | !< |
|---|
| 39 | |
|---|
| 40 | |
|---|
| 41 | !> \file main3D-0.4-40km.f90 GRISLI Modele 3D De Calotte Glaciaire |
|---|
| 42 | !! programme principal |
|---|
| 43 | !! (voir l'\ref tree) |
|---|
| 44 | !! |
|---|
| 45 | !! @brief modele flow line d'evolution de calotte |
|---|
| 46 | !! @authors Catherine Ritz catritz@lgge.obs.ujf-grenoble.fr (tout du long) |
|---|
| 47 | !! @authors Adeline Fabre (la partie Gremlins) |
|---|
| 48 | !! @authors Vincent Rommelaere (ice shelves et ice streams) |
|---|
| 49 | !! @authors Christophe Dumas (debut f90, (Antarctique) |
|---|
| 50 | !! @authors Vincent Peyaud (portage HN,calving, front, hydrol) |
|---|
| 51 | !! @authors Cyril Mazauric (AGRIF) |
|---|
| 52 | !! |
|---|
| 53 | !! |
|---|
| 54 | !! @note use module3D_phy |
|---|
| 55 | !! @note use module_choix |
|---|
| 56 | !! @note use flottab_mod |
|---|
| 57 | !! @note use icetempmod |
|---|
| 58 | !! @note use sorties_ncdf_grisli |
|---|
| 59 | !! @note use diagno_mod |
|---|
| 60 | !! @note use resolmeca_SIA_L1 |
|---|
| 61 | !! |
|---|
| 62 | !! |
|---|
| 63 | !! Ce module appelle les routines suivantes : |
|---|
| 64 | !! - grisli_init() |
|---|
| 65 | !! - step_grisli1() |
|---|
| 66 | !! - step_output() |
|---|
| 67 | !! |
|---|
| 68 | !< |
|---|
| 69 | |
|---|
| 70 | !> \namespace main3D GRISLI Modele 3D De Calotte Glaciaire |
|---|
| 71 | !! programme principal |
|---|
| 72 | !! (voir l'\ref tree) |
|---|
| 73 | !! |
|---|
| 74 | !! |
|---|
| 75 | !! @brief modele flow line d'evolution de calotte |
|---|
| 76 | !! @authors Catherine Ritz catritz@lgge.obs.ujf-grenoble.fr (tout du long) |
|---|
| 77 | !! @authors Adeline Fabre (la partie Gremlins) |
|---|
| 78 | !! @authors Vincent Rommelaere (ice shelves et ice streams) |
|---|
| 79 | !! @authors Christophe Dumas (debut f90, (Antarctique) |
|---|
| 80 | !! @authors Vincent Peyaud (portage HN,calving, front, hydrol) |
|---|
| 81 | !! @authors Cyril Mazauric (AGRIF) |
|---|
| 82 | !! |
|---|
| 83 | !! |
|---|
| 84 | !! @note use module3D_phy |
|---|
| 85 | !! @note use module_choix |
|---|
| 86 | !! @note use flottab_mod |
|---|
| 87 | !! @note use icetempmod |
|---|
| 88 | !! @note use sorties_ncdf_grisli |
|---|
| 89 | !! @note use diagno_mod |
|---|
| 90 | !! @note use resolmeca_SIA_L1 |
|---|
| 91 | !! |
|---|
| 92 | !! |
|---|
| 93 | !! @todo itracebug : faire une routine |
|---|
| 94 | !! |
|---|
| 95 | !! Ce module appelle les routines suivantes : |
|---|
| 96 | !! - grisli_init() |
|---|
| 97 | !! - step_grisli1() |
|---|
| 98 | !! |
|---|
| 99 | !! Defined in file main3D-0.4-40km.f90 |
|---|
| 100 | !< |
|---|
| 101 | |
|---|
| 102 | program main3D |
|---|
| 103 | |
|---|
| 104 | use module3D_phy, only: time,iout,ntmax |
|---|
| 105 | use runparam, only: tend,itracebug,nt |
|---|
| 106 | |
|---|
| 107 | implicit none |
|---|
| 108 | |
|---|
| 109 | |
|---|
| 110 | ! good luck |
|---|
| 111 | |
|---|
| 112 | call grisli_init ! Initializations |
|---|
| 113 | |
|---|
| 114 | time_loop: do nt=1,ntmax !____________________________ debut boucle temporelle |
|---|
| 115 | |
|---|
| 116 | if (time.ge.tend) exit |
|---|
| 117 | if (time.gt.10) itracebug = 0 |
|---|
| 118 | call step_time_loop() |
|---|
| 119 | |
|---|
| 120 | end do time_loop |
|---|
| 121 | if (itracebug.eq.1) call tracebug('dans main avant call out_recovery ') |
|---|
| 122 | call out_recovery(iout) |
|---|
| 123 | |
|---|
| 124 | write(6,*) "end of the run at time = ",time |
|---|
| 125 | write(6,*) "_____________________________________________________________________" |
|---|
| 126 | |
|---|
| 127 | end program main3D |
|---|
| 128 | |
|---|
| 129 | |
|---|
| 130 | !--------------------------------------------------------------------------------------- |
|---|
| 131 | subroutine grisli_init |
|---|
| 132 | |
|---|
| 133 | USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen,num_sealevel, & |
|---|
| 134 | num_ts_ritz,num_ic_vo,num_ic_by,num_ic_dm,num_ic_dc,num_ic_df, & |
|---|
| 135 | s,h,b,bsoc,flot,mk,mk0,uxbar,uybar,hwater,time,timemax,ndebug,ndebug_max |
|---|
| 136 | use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,runname,itracebug |
|---|
| 137 | use geography, only: nx,ny,geoplace |
|---|
| 138 | use deformation_mod_2lois, only:n1poly,n2poly |
|---|
| 139 | use bilan_eau_mod, only: init_bilan_eau |
|---|
| 140 | use module_choix, only: forclim,ablation,bmeltshelf,calving,flow_general,flowlaw |
|---|
| 141 | ! module_choix donne acces a tous les modules |
|---|
| 142 | ! de declaration des packages |
|---|
| 143 | use flottab_mod, only: flottab |
|---|
| 144 | use sorties_ncdf_grisli, only: iglob_ncdf,testsort_time_ncdf,init_sortie_ncdf,testsort_time_ncdf, & |
|---|
| 145 | sortie_ncdf_cat |
|---|
| 146 | use util_recovery, only: dtout |
|---|
| 147 | |
|---|
| 148 | ! use track_debug |
|---|
| 149 | |
|---|
| 150 | implicit none |
|---|
| 151 | |
|---|
| 152 | integer :: i,j |
|---|
| 153 | |
|---|
| 154 | if (itracebug.eq.1) call tracebug(' Entree dans routine grisli_init') |
|---|
| 155 | ! switch pour passer ou non par T lliboutry calcule => 0, ne passe pas, |
|---|
| 156 | ! 1 ou 2 passe (se met a 0 tout seul si on prend un fichier .cptr) |
|---|
| 157 | |
|---|
| 158 | ITEMP=0 |
|---|
| 159 | |
|---|
| 160 | ! switch couple physique faible => CP et CT independant T |
|---|
| 161 | ! 0 pas de trait. vert. A FAIRE niveau L0 |
|---|
| 162 | ! 1 pas de couplage , faible physique niveau L1 |
|---|
| 163 | ! 2 couplage, faible physique niveau L2 |
|---|
| 164 | ! 3 couplage, physique complete sans CBT niveau L3 |
|---|
| 165 | ! 4 idem 3 mais loi de def. Duval niveau L4 |
|---|
| 166 | ICOUPLE=4 |
|---|
| 167 | |
|---|
| 168 | nt=-1 ! utilisee dans initialisation flottab |
|---|
| 169 | ! sortie profile tous les dtprofile |
|---|
| 170 | DTPROFILE=50000. |
|---|
| 171 | ! ----------------------------------fin des modifs run les plus usuelles |
|---|
| 172 | DIRNAMEOUT='../RESULTATS/' |
|---|
| 173 | !DIRNAMEOUT='./' |
|---|
| 174 | |
|---|
| 175 | call initial ! routine qui appel toutes les routines d'initialisation |
|---|
| 176 | |
|---|
| 177 | |
|---|
| 178 | ! call init_sortie_ncdf |
|---|
| 179 | ! call sortie_ncdf_cat |
|---|
| 180 | ! STOP |
|---|
| 181 | |
|---|
| 182 | ! compteur tous les DTCPT |
|---|
| 183 | DTCPT=dtout |
|---|
| 184 | |
|---|
| 185 | |
|---|
| 186 | ! ************ OPEN FILES.RITZ **************** |
|---|
| 187 | |
|---|
| 188 | if ((geoplace.eq.'anteis1').or.(geoplace.eq.'ant20km')) then |
|---|
| 189 | ! fichier de reference pour le niveau des mers |
|---|
| 190 | open(num_sealevel,file=TRIM(DIRNAMEOUT)//'sealevel'//runname//'.ritz',position="append") |
|---|
| 191 | open(num_ts_ritz,file=TRIM(DIRNAMEOUT)//'ts_'//runname//'.ritz',position="append") |
|---|
| 192 | open(num_ic_vo,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'vo.ritz',position="append") |
|---|
| 193 | open(num_ic_by,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'by.ritz',position="append") |
|---|
| 194 | open(num_ic_dm,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'dm.ritz',position="append") |
|---|
| 195 | open(num_ic_dc,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'dc.ritz',position="append") |
|---|
| 196 | open(num_ic_df,file=TRIM(DIRNAMEOUT)//'ic_'//runname//'df.ritz',position="append") |
|---|
| 197 | endif |
|---|
| 198 | |
|---|
| 199 | !------------------------------ INITIALISATION ---------------------------- |
|---|
| 200 | ! |
|---|
| 201 | ! ecriture netcdf apres initialisation |
|---|
| 202 | |
|---|
| 203 | |
|---|
| 204 | |
|---|
| 205 | call testsort_time_ncdf(dble(tbegin)) |
|---|
| 206 | if (iglob_ncdf .EQ. 1) call sortie_ncdf_cat |
|---|
| 207 | |
|---|
| 208 | ! test vincent car certains H(i,j)=0 dans fichier de reprise |
|---|
| 209 | do j=1,ny |
|---|
| 210 | do i=1,nx |
|---|
| 211 | H(i,j)=max(0.,H(i,j)) |
|---|
| 212 | enddo |
|---|
| 213 | enddo |
|---|
| 214 | |
|---|
| 215 | |
|---|
| 216 | call forclim ! initialisation BM et TS |
|---|
| 217 | call ablation |
|---|
| 218 | |
|---|
| 219 | |
|---|
| 220 | |
|---|
| 221 | ! ----------- CALCULATION OF INITIAL TEMPERATURES |
|---|
| 222 | |
|---|
| 223 | tcpt:if (ICOMPTEUR.eq.0) then |
|---|
| 224 | |
|---|
| 225 | |
|---|
| 226 | |
|---|
| 227 | if ((ITEMP.eq.0).or.(ITEMP.eq.3)) then |
|---|
| 228 | call masque(flot,mk,mk0,itracebug) |
|---|
| 229 | |
|---|
| 230 | call Neffect() |
|---|
| 231 | |
|---|
| 232 | call flottab() |
|---|
| 233 | |
|---|
| 234 | call Neffect() |
|---|
| 235 | |
|---|
| 236 | |
|---|
| 237 | ! call vitbilan_lect ! routine de lecture des vitesses de bilan |
|---|
| 238 | ! ======================================================== |
|---|
| 239 | |
|---|
| 240 | if (ITEMP.eq.0) call lineartemp() |
|---|
| 241 | |
|---|
| 242 | call bmelt_grounded |
|---|
| 243 | call bmeltshelf |
|---|
| 244 | |
|---|
| 245 | |
|---|
| 246 | call flow_general |
|---|
| 247 | |
|---|
| 248 | do iglen=n1poly,n2poly |
|---|
| 249 | call flowlaw(iglen) |
|---|
| 250 | end do |
|---|
| 251 | |
|---|
| 252 | call Neffect() |
|---|
| 253 | call flottab() |
|---|
| 254 | call calving |
|---|
| 255 | call ablation_bord |
|---|
| 256 | call flottab |
|---|
| 257 | call Neffect() |
|---|
| 258 | call diffusiv() |
|---|
| 259 | call SIA_velocities() |
|---|
| 260 | endif |
|---|
| 261 | |
|---|
| 262 | else ! tcpt on reprend un fichier compteur (ICOMPTEUR.eq.1) |
|---|
| 263 | |
|---|
| 264 | time=tbegin ! prend le temps du compteur |
|---|
| 265 | |
|---|
| 266 | |
|---|
| 267 | call masque(flot,mk,mk0,itracebug) |
|---|
| 268 | call flottab() |
|---|
| 269 | call neffect() |
|---|
| 270 | call flottab() |
|---|
| 271 | call masque(flot,mk,mk0,itracebug) |
|---|
| 272 | |
|---|
| 273 | do i=1,nx |
|---|
| 274 | do j=1,ny |
|---|
| 275 | if (S(i,j).lt.0) then |
|---|
| 276 | print*,i,j,S(i,j) |
|---|
| 277 | goto 11115 |
|---|
| 278 | endif |
|---|
| 279 | enddo |
|---|
| 280 | enddo |
|---|
| 281 | 11115 continue |
|---|
| 282 | |
|---|
| 283 | call bmeltshelf ! afq -- |
|---|
| 284 | |
|---|
| 285 | ! ======================================================== |
|---|
| 286 | call flow_general |
|---|
| 287 | |
|---|
| 288 | do iglen=n1poly,n2poly |
|---|
| 289 | call flowlaw(iglen) |
|---|
| 290 | end do |
|---|
| 291 | |
|---|
| 292 | |
|---|
| 293 | call Neffect() |
|---|
| 294 | call flottab |
|---|
| 295 | call diffusiv() |
|---|
| 296 | call SIA_velocities() |
|---|
| 297 | call strain_rate |
|---|
| 298 | |
|---|
| 299 | endif tcpt |
|---|
| 300 | ! fin du test sur icompteur |
|---|
| 301 | |
|---|
| 302 | ! call init_sortie_ncdf |
|---|
| 303 | ! call sortie_ncdf_cat |
|---|
| 304 | |
|---|
| 305 | call flottab() |
|---|
| 306 | call Neffect() |
|---|
| 307 | call flottab() |
|---|
| 308 | |
|---|
| 309 | if (icompteur.eq.0) then |
|---|
| 310 | do i=1,nx |
|---|
| 311 | do j=1,ny |
|---|
| 312 | if (.not.flot(i,j)) then |
|---|
| 313 | B(i,j) = Bsoc(i,j) |
|---|
| 314 | Uxbar(i,j) = 0. |
|---|
| 315 | Uybar(i,j) = 0. |
|---|
| 316 | end if |
|---|
| 317 | end do |
|---|
| 318 | end do |
|---|
| 319 | endif |
|---|
| 320 | |
|---|
| 321 | do i=2,nx-1 |
|---|
| 322 | do j=2,ny-1 |
|---|
| 323 | hwater(i,j)=max(hwater(i,j),0.) |
|---|
| 324 | enddo |
|---|
| 325 | enddo |
|---|
| 326 | timemax=time |
|---|
| 327 | isynchro=1 |
|---|
| 328 | ndebug=0 |
|---|
| 329 | ndebug_max=9 |
|---|
| 330 | |
|---|
| 331 | call init_bilan_eau |
|---|
| 332 | call step_thermomeca() ! un tour dans la boucle temporelle, partie avant icethick |
|---|
| 333 | call init_sortie_ncdf |
|---|
| 334 | if (itracebug.eq.1) call tracebug(' fin routine grisli_init') |
|---|
| 335 | call testsort_time_ncdf(dble(tbegin)) |
|---|
| 336 | |
|---|
| 337 | if (iglob_ncdf .EQ. 1) call sortie_ncdf_cat |
|---|
| 338 | |
|---|
| 339 | return |
|---|
| 340 | end subroutine grisli_init |
|---|
![(please configure the [header_logo] section in trac.ini)](/grisli/chrome/site/Laurentide-16km-3D-Umag-water-bordure-GRISLI-small.png)
