Changeset 1939 for LMDZ5/trunk
- Timestamp:
- Jan 21, 2014, 3:23:17 PM (11 years ago)
- Location:
- LMDZ5/trunk/libf
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d/fxhyp.F
r1930 r1939 307 307 308 308 309 310 309 IF(ik.EQ.1.and.grossism.EQ.1.) THEN 311 310 xvrai(1) = xvrai(iip1)-depi -
LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.F
r1908 r1939 68 68 & 'dynetat0_loc: Pb d''ouverture du fichier start.nc' 69 69 write(lunout,*)' ierr = ', ierr 70 CALL ABORT 70 CALL ABORT_GCM("DYNETAT0", "", 1) 71 71 ENDIF 72 72 … … 75 75 IF (ierr .NE. NF_NOERR) THEN 76 76 write(lunout,*)"dynetat0_loc: Le champ <controle> est absent" 77 CALL abort 77 CALL abort_gcm("dynetat0", "", 1) 78 78 ENDIF 79 79 #ifdef NC_DOUBLE … … 84 84 IF (ierr .NE. NF_NOERR) THEN 85 85 write(lunout,*)"dynetat0_loc: Lecture echoue pour <controle>" 86 CALL abort 86 CALL abort_gcm("dynetat0", "", 1) 87 87 ENDIF 88 88 … … 158 158 IF (ierr .NE. NF_NOERR) THEN 159 159 write(lunout,*)"dynetat0_loc: Le champ <rlonu> est absent" 160 CALL abort 160 CALL abort_gcm("dynetat0", "", 1) 161 161 ENDIF 162 162 #ifdef NC_DOUBLE … … 167 167 IF (ierr .NE. NF_NOERR) THEN 168 168 write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlonu>" 169 CALL abort 169 CALL abort_gcm("dynetat0", "", 1) 170 170 ENDIF 171 171 … … 173 173 IF (ierr .NE. NF_NOERR) THEN 174 174 write(lunout,*)"dynetat0_loc: Le champ <rlatu> est absent" 175 CALL abort 175 CALL abort_gcm("dynetat0", "", 1) 176 176 ENDIF 177 177 #ifdef NC_DOUBLE … … 182 182 IF (ierr .NE. NF_NOERR) THEN 183 183 write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlatu>" 184 CALL abort 184 CALL abort_gcm("dynetat0", "", 1) 185 185 ENDIF 186 186 … … 188 188 IF (ierr .NE. NF_NOERR) THEN 189 189 write(lunout,*)"dynetat0_loc: Le champ <rlonv> est absent" 190 CALL abort 190 CALL abort_gcm("dynetat0", "", 1) 191 191 ENDIF 192 192 #ifdef NC_DOUBLE … … 197 197 IF (ierr .NE. NF_NOERR) THEN 198 198 write(lunout,*)"dynetat0_loc: Lecture echouee pour <rlonv>" 199 CALL abort 199 CALL abort_gcm("dynetat0", "", 1) 200 200 ENDIF 201 201 … … 203 203 IF (ierr .NE. NF_NOERR) THEN 204 204 write(lunout,*)"dynetat0_loc: Le champ <rlatv> est absent" 205 CALL abort 205 CALL abort_gcm("dynetat0", "", 1) 206 206 ENDIF 207 207 #ifdef NC_DOUBLE … … 212 212 IF (ierr .NE. NF_NOERR) THEN 213 213 write(lunout,*)"dynetat0_loc: Lecture echouee pour rlatv" 214 CALL abort 214 CALL abort_gcm("dynetat0", "", 1) 215 215 ENDIF 216 216 … … 218 218 IF (ierr .NE. NF_NOERR) THEN 219 219 write(lunout,*)"dynetat0_loc: Le champ <cu> est absent" 220 CALL abort 220 CALL abort_gcm("dynetat0", "", 1) 221 221 ENDIF 222 222 #ifdef NC_DOUBLE … … 227 227 IF (ierr .NE. NF_NOERR) THEN 228 228 write(lunout,*)"dynetat0_loc: Lecture echouee pour <cu>" 229 CALL abort 229 CALL abort_gcm("dynetat0", "", 1) 230 230 ENDIF 231 231 … … 233 233 IF (ierr .NE. NF_NOERR) THEN 234 234 write(lunout,*)"dynetat0_loc: Le champ <cv> est absent" 235 CALL abort 235 CALL abort_gcm("dynetat0", "", 1) 236 236 ENDIF 237 237 #ifdef NC_DOUBLE … … 242 242 IF (ierr .NE. NF_NOERR) THEN 243 243 write(lunout,*)"dynetat0_loc: Lecture echouee pour <cv>" 244 CALL abort 244 CALL abort_gcm("dynetat0", "", 1) 245 245 ENDIF 246 246 … … 248 248 IF (ierr .NE. NF_NOERR) THEN 249 249 write(lunout,*)"dynetat0_loc: Le champ <aire> est absent" 250 CALL abort 250 CALL abort_gcm("dynetat0", "", 1) 251 251 ENDIF 252 252 #ifdef NC_DOUBLE … … 257 257 IF (ierr .NE. NF_NOERR) THEN 258 258 write(lunout,*)"dynetat0_loc: Lecture echouee pour <aire>" 259 CALL abort 259 CALL abort_gcm("dynetat0", "", 1) 260 260 ENDIF 261 261 … … 265 265 IF (ierr .NE. NF_NOERR) THEN 266 266 write(lunout,*)"dynetat0_loc: Le champ <phisinit> est absent" 267 CALL abort 267 CALL abort_gcm("dynetat0", "", 1) 268 268 ENDIF 269 269 #ifdef NC_DOUBLE … … 274 274 IF (ierr .NE. NF_NOERR) THEN 275 275 write(lunout,*)"dynetat0_loc: Lecture echouee pour <phisinit>" 276 CALL abort 276 CALL abort_gcm("dynetat0", "", 1) 277 277 ENDIF 278 278 phis(ijb_u:ije_u)=phis_glo(ijb_u:ije_u) … … 286 286 IF (ierr .NE. NF_NOERR) THEN 287 287 write(lunout,*)"dynetat0: Le champ <Time> est absent" 288 CALL abort 288 CALL abort_gcm("dynetat0", "", 1) 289 289 ENDIF 290 290 ENDIF … … 296 296 IF (ierr .NE. NF_NOERR) THEN 297 297 write(lunout,*)"dynetat0_loc: Lecture echouee <temps>" 298 CALL abort 298 CALL abort_gcm("dynetat0", "", 1) 299 299 ENDIF 300 300 … … 302 302 IF (ierr .NE. NF_NOERR) THEN 303 303 write(lunout,*)"dynetat0_loc: Le champ <ucov> est absent" 304 CALL abort 304 CALL abort_gcm("dynetat0", "", 1) 305 305 ENDIF 306 306 … … 314 314 IF (ierr .NE. NF_NOERR) THEN 315 315 write(lunout,*)"dynetat0_loc: Lecture echouee pour <ucov>" 316 CALL abort 316 CALL abort_gcm("dynetat0", "", 1) 317 317 ENDIF 318 318 … … 324 324 IF (ierr .NE. NF_NOERR) THEN 325 325 write(lunout,*)"dynetat0_loc: Le champ <vcov> est absent" 326 CALL abort 326 CALL abort_gcm("dynetat0", "", 1) 327 327 ENDIF 328 328 #ifdef NC_DOUBLE … … 333 333 IF (ierr .NE. NF_NOERR) THEN 334 334 write(lunout,*)"dynetat0_loc: Lecture echouee pour <vcov>" 335 CALL abort 335 CALL abort_gcm("dynetat0", "", 1) 336 336 ENDIF 337 337 vcov(ijb_v:ije_v,:)=vcov_glo(ijb_v:ije_v,:) … … 342 342 IF (ierr .NE. NF_NOERR) THEN 343 343 write(lunout,*)"dynetat0_loc: Le champ <teta> est absent" 344 CALL abort 344 CALL abort_gcm("dynetat0", "", 1) 345 345 ENDIF 346 346 #ifdef NC_DOUBLE … … 351 351 IF (ierr .NE. NF_NOERR) THEN 352 352 write(lunout,*)"dynetat0_loc: Lecture echouee pour <teta>" 353 CALL abort 353 CALL abort_gcm("dynetat0", "", 1) 354 354 ENDIF 355 355 … … 375 375 write(lunout,*) 376 376 & "dynetat0_loc: Lecture echouee pour "//tname(iq) 377 CALL abort 377 CALL abort_gcm("dynetat0", "", 1) 378 378 ENDIF 379 379 q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:) … … 388 388 IF (ierr .NE. NF_NOERR) THEN 389 389 write(lunout,*)"dynetat0_loc: Le champ <masse> est absent" 390 CALL abort 390 CALL abort_gcm("dynetat0", "", 1) 391 391 ENDIF 392 392 #ifdef NC_DOUBLE … … 397 397 IF (ierr .NE. NF_NOERR) THEN 398 398 write(lunout,*)"dynetat0_loc: Lecture echouee pour <masse>" 399 CALL abort 399 CALL abort_gcm("dynetat0", "", 1) 400 400 ENDIF 401 401 masse(ijb_u:ije_u,:)=masse_glo(ijb_u:ije_u,:) … … 406 406 IF (ierr .NE. NF_NOERR) THEN 407 407 write(lunout,*)"dynetat0_loc: Le champ <ps> est absent" 408 CALL abort 408 CALL abort_gcm("dynetat0", "", 1) 409 409 ENDIF 410 410 #ifdef NC_DOUBLE … … 415 415 IF (ierr .NE. NF_NOERR) THEN 416 416 write(lunout,*)"dynetat0_loc: Lecture echouee pour <ps>" 417 CALL abort 417 CALL abort_gcm("dynetat0", "", 1) 418 418 ENDIF 419 419 -
LMDZ5/trunk/libf/dyn3dmem/dynredem_loc.F
r1907 r1939 138 138 & //trim(fichnom) 139 139 write(lunout,*)' ierr = ', ierr 140 CALL ABORT 140 CALL ABORT_GCM("DYNREDEM0", "", 1) 141 141 ENDIF 142 142 c … … 588 588 IF (ierr .NE. NF_NOERR) THEN 589 589 write(lunout,*)"dynredem1: Pb. d ouverture "//trim(fichnom) 590 CALL abort 590 CALL abort_gcm("dynredem1", "", 1) 591 591 ENDIF 592 592 -
LMDZ5/trunk/libf/dyn3dmem/fxhyp.F
r1907 r1939 48 48 c 49 49 REAL dzoom 50 REAL *8xlon(iip1),xprimm(iip1),xuv51 REAL *8xtild(0:nmax2)52 REAL *8fhyp(0:nmax2),ffdx,beta,Xprimt(0:nmax2)53 REAL *8Xf(0:nmax2),xxpr(0:nmax2)54 REAL *8xvrai(iip1),xxprim(iip1)55 REAL *8pi,depi,epsilon,xzoom,fa,fb56 REAL *8Xf1, Xfi , a0,a1,a2,a3,xi250 REAL(KIND=8) xlon(iip1),xprimm(iip1),xuv 51 REAL(KIND=8) xtild(0:nmax2) 52 REAL(KIND=8) fhyp(0:nmax2),ffdx,beta,Xprimt(0:nmax2) 53 REAL(KIND=8) Xf(0:nmax2),xxpr(0:nmax2) 54 REAL(KIND=8) xvrai(iip1),xxprim(iip1) 55 REAL(KIND=8) pi,depi,epsilon,xzoom,fa,fb 56 REAL(KIND=8) Xf1, Xfi , a0,a1,a2,a3,xi2 57 57 INTEGER i,it,ik,iter,ii,idif,ii1,ii2 58 REAL *8xi,xo1,xmoy,xlon2,fxm,Xprimin59 REAL *8champmin,champmax,decalx58 REAL(KIND=8) xi,xo1,xmoy,xlon2,fxm,Xprimin 59 REAL(KIND=8) champmin,champmax,decalx 60 60 INTEGER is2 61 61 SAVE is2 62 62 63 REAL *8heavyside63 REAL(KIND=8) heavyside 64 64 65 65 pi = 2. * ASIN(1.) … … 157 157 WRITE(6,*)'Modifier les valeurs de grossismx ,tau ou dzoomx ', 158 158 , ' et relancer ! *** ' 159 CALL ABORT 159 CALL ABORT_GCM("FXHYP", "", 1) 160 160 ENDIF 161 161 c … … 290 290 xxprim(1) = xxprim(iip1) 291 291 ENDIF 292 292 293 DO i = 1 , iim 293 294 xlon(i) = xvrai(i) -
LMDZ5/trunk/libf/dyn3dpar/fxhyp.F
r1907 r1939 69 69 c 70 70 if (iim==1) then 71 72 print*,'Longitudes calculees a la main pour iim=1'73 71 74 72 rlonm025(1)=-pi/2. … … 180 178 WRITE(6,*)'Modifier les valeurs de grossismx ,tau ou dzoomx ', 181 179 , ' et relancer ! *** ' 182 CALL ABORT 180 CALL ABORT_GCM("FXHYP", "", 1) 183 181 ENDIF 184 182 c … … 307 305 308 306 1500 CONTINUE 309 310 307 311 308 -
LMDZ5/trunk/libf/dyn3dpar/gcm.F
r1907 r1939 332 332 start_time = starttime 333 333 ELSE 334 WRITE(lunout,*)'Je m''arrete' 335 CALL abort 334 call abort_gcm("gcm", "'Je m''arrete'", 1) 336 335 ENDIF 337 336 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.