- Timestamp:
- 01/19/17 15:21:53 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/New-remplimat/diagno-L2_mod.f90
r96 r98 28 28 29 29 integer :: ifail_diagno ! pour recuperation d'erreur 30 integer :: ifail_diagno_ramollo ! pour recuperation d'erreur shelf ramollo 30 31 integer :: iplus1,jplus1 31 32 integer :: ctvisco,iumax,jumax … … 140 141 !~ enddo 141 142 142 if (Schoof.eq.1 ) then ! flux grounding line Schoof143 if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof 143 144 call interpol_glflux ! calcul flux GL + interpolation sur voisins 144 145 endif … … 180 181 !nyd2=60 181 182 182 183 183 call rempli_L2(nxd1,nxd2,nyd1,nyd2,uxbar(nxd1:nxd2,nyd1:nyd2),uybar(nxd1:nxd2,nyd1:nyd2), & 184 184 uxb1(nxd1:nxd2,nyd1:nyd2),uyb1(nxd1:nxd2,nyd1:nyd2), & 185 185 imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno) 186 186 187 if (Schoof.eq.1 ) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo187 if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo 188 188 pvi_keep(:,:)=pvi(:,:) 189 where (flot(:,:)) 190 pvi(:,:)=pvimin 191 endwhere 192 189 where (flot(:,:).and.H(:,:).GT.2.) 190 pvi(:,:)=1.e5 191 ! pvi(:,:)=pvimin 192 endwhere 193 193 194 call rempli_L2(nxd1,nxd2,nyd1,nyd2,uxbar(nxd1:nxd2,nyd1:nyd2),uybar(nxd1:nxd2,nyd1:nyd2), & 194 195 uxb1ramollo(nxd1:nxd2,nyd1:nyd2),uyb1ramollo(nxd1:nxd2,nyd1:nyd2), & 195 imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno) 196 imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno_ramollo) 197 196 198 pvi(:,:)=pvi_keep(:,:) 199 197 200 where ((uxb1ramollo(:,:)**2 + uyb1ramollo(:,:)**2).GT.0.) 198 201 back_force(:,:) = sqrt(uxb1(:,:)**2+ uyb1(:,:)**2) / sqrt(uxb1ramollo(:,:)**2 + uyb1ramollo(:,:)**2) … … 200 203 back_force(:,:)=1. 201 204 endwhere 202 endif 203 205 206 if (ifail_diagno_ramollo.gt.0) then 207 write(6,*) ' Probleme resolution systeme L2. ramollo ifail=',ifail_diagno_ramollo 208 STOP 209 endif 204 210 !~ do j=1,ny 205 211 !~ do i=1,nx … … 214 220 !~ print*,'apres calcul rempli_L2' 215 221 !~ read(*,*) 222 endif 223 216 224 217 225 ! Dans rempli_L2
Note: See TracChangeset
for help on using the changeset viewer.