Changeset 481 for branches/GRISLIv3/SOURCES/main3D-0.4-40km.f90
- Timestamp:
- 02/09/24 10:58:55 (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/GRISLIv3/SOURCES/main3D-0.4-40km.f90
r476 r481 131 131 subroutine grisli_init 132 132 133 USE module3D_phy, only: i temp,icouple,isynchro,icompteur,iglen, &133 USE module3D_phy, only: isynchro,icompteur,iglen, & 134 134 ice,bm,bmelt,ablbord,ablbord_dtt,dt, & 135 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,itracebug136 use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,itracebug 137 137 use geography, only: nx,ny,geoplace 138 138 use deformation_mod_2lois, only:n1poly,n2poly … … 156 156 ! 1 ou 2 passe (se met a 0 tout seul si on prend un fichier .cptr) 157 157 158 ITEMP=0159 160 ! switch couple physique faible => CP et CT independant T161 ! 0 pas de trait. vert. A FAIRE niveau L0162 ! 1 pas de couplage , faible physique niveau L1163 ! 2 couplage, faible physique niveau L2164 ! 3 couplage, physique complete sans CBT niveau L3165 ! 4 idem 3 mais loi de def. Duval niveau L4166 ICOUPLE=4167 168 158 nt=-1 ! utilisee dans initialisation flottab 169 159 ! sortie profile tous les dtprofile … … 209 199 tcpt:if (ICOMPTEUR.eq.0) then 210 200 211 212 213 if ((ITEMP.eq.0).or.(ITEMP.eq.3)) then214 201 call masque(flot,mk,mk0,itracebug) 215 202 … … 224 211 ! ======================================================== 225 212 226 if (ITEMP.eq.0)call lineartemp()213 call lineartemp() 227 214 228 215 call bmelt_grounded … … 244 231 call diffusiv() 245 232 call SIA_velocities() 246 endif247 233 248 234 else ! tcpt on reprend un fichier compteur (ICOMPTEUR.eq.1)
Note: See TracChangeset
for help on using the changeset viewer.