Changeset 1930 for LMDZ5/trunk/libf/dyn3d
- Timestamp:
- Jan 17, 2014, 5:45:09 PM (11 years ago)
- Location:
- LMDZ5/trunk/libf/dyn3d
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d/defrun.F
r1907 r1930 364 364 WRITE(tapeout,*)' La valeur de dzoomx passee par run.def est dif 365 365 *ferente de celle lue sur le fichier start ' 366 CALL ABORT 366 CALL ABORT_gcm("defrun", "", 1) 367 367 ENDIF 368 368 … … 370 370 WRITE(tapeout,*)' La valeur de dzoomy passee par run.def est dif 371 371 *ferente de celle lue sur le fichier start ' 372 CALL ABORT 372 CALL ABORT_gcm("defrun", "", 1) 373 373 ENDIF 374 374 … … 376 376 WRITE(6,*)' La valeur de taux passee par run.def est differente 377 377 * de celle lue sur le fichier start ' 378 CALL ABORT 378 CALL ABORT_gcm("defrun", "", 1) 379 379 ENDIF 380 380 … … 382 382 WRITE(6,*)' La valeur de tauy passee par run.def est differente 383 383 * de celle lue sur le fichier start ' 384 CALL ABORT 384 CALL ABORT_gcm("defrun", "", 1) 385 385 ENDIF 386 386 -
LMDZ5/trunk/libf/dyn3d/dynetat0.F
r1907 r1930 67 67 write(lunout,*)'dynetat0: Pb d''ouverture du fichier start.nc' 68 68 write(lunout,*)' ierr = ', ierr 69 CALL ABORT 69 CALL ABORT_gcm("dynetat0", "", 1) 70 70 ENDIF 71 71 … … 74 74 IF (ierr .NE. NF_NOERR) THEN 75 75 write(lunout,*)"dynetat0: Le champ <controle> est absent" 76 CALL abort76 CALL ABORT_gcm("dynetat0", "", 1) 77 77 ENDIF 78 78 ierr = nf90_get_var(nid, nvarid, tab_cntrl) 79 79 IF (ierr .NE. NF_NOERR) THEN 80 80 write(lunout,*)"dynetat0: Lecture echoue pour <controle>" 81 CALL abort81 CALL ABORT_gcm("dynetat0", "", 1) 82 82 ENDIF 83 83 … … 154 154 IF (ierr .NE. NF_NOERR) THEN 155 155 write(lunout,*)"dynetat0: Le champ <rlonu> est absent" 156 CALL abort156 CALL ABORT_gcm("dynetat0", "", 1) 157 157 ENDIF 158 158 ierr = nf90_get_var(nid, nvarid, rlonu) 159 159 IF (ierr .NE. NF_NOERR) THEN 160 160 write(lunout,*)"dynetat0: Lecture echouee pour <rlonu>" 161 CALL abort161 CALL ABORT_gcm("dynetat0", "", 1) 162 162 ENDIF 163 163 … … 165 165 IF (ierr .NE. NF_NOERR) THEN 166 166 write(lunout,*)"dynetat0: Le champ <rlatu> est absent" 167 CALL abort167 CALL ABORT_gcm("dynetat0", "", 1) 168 168 ENDIF 169 169 ierr = nf90_get_var(nid, nvarid, rlatu) 170 170 IF (ierr .NE. NF_NOERR) THEN 171 171 write(lunout,*)"dynetat0: Lecture echouee pour <rlatu>" 172 CALL abort172 CALL ABORT_gcm("dynetat0", "", 1) 173 173 ENDIF 174 174 … … 176 176 IF (ierr .NE. NF_NOERR) THEN 177 177 write(lunout,*)"dynetat0: Le champ <rlonv> est absent" 178 CALL abort178 CALL ABORT_gcm("dynetat0", "", 1) 179 179 ENDIF 180 180 ierr = nf90_get_var(nid, nvarid, rlonv) 181 181 IF (ierr .NE. NF_NOERR) THEN 182 182 write(lunout,*)"dynetat0: Lecture echouee pour <rlonv>" 183 CALL abort183 CALL ABORT_gcm("dynetat0", "", 1) 184 184 ENDIF 185 185 … … 187 187 IF (ierr .NE. NF_NOERR) THEN 188 188 write(lunout,*)"dynetat0: Le champ <rlatv> est absent" 189 CALL abort189 CALL ABORT_gcm("dynetat0", "", 1) 190 190 ENDIF 191 191 ierr = nf90_get_var(nid, nvarid, rlatv) 192 192 IF (ierr .NE. NF_NOERR) THEN 193 193 write(lunout,*)"dynetat0: Lecture echouee pour rlatv" 194 CALL abort194 CALL ABORT_gcm("dynetat0", "", 1) 195 195 ENDIF 196 196 … … 198 198 IF (ierr .NE. NF_NOERR) THEN 199 199 write(lunout,*)"dynetat0: Le champ <cu> est absent" 200 CALL abort200 CALL ABORT_gcm("dynetat0", "", 1) 201 201 ENDIF 202 202 ierr = nf90_get_var(nid, nvarid, cu) 203 203 IF (ierr .NE. NF_NOERR) THEN 204 204 write(lunout,*)"dynetat0: Lecture echouee pour <cu>" 205 CALL abort205 CALL ABORT_gcm("dynetat0", "", 1) 206 206 ENDIF 207 207 … … 209 209 IF (ierr .NE. NF_NOERR) THEN 210 210 write(lunout,*)"dynetat0: Le champ <cv> est absent" 211 CALL abort211 CALL ABORT_gcm("dynetat0", "", 1) 212 212 ENDIF 213 213 ierr = nf90_get_var(nid, nvarid, cv) 214 214 IF (ierr .NE. NF_NOERR) THEN 215 215 write(lunout,*)"dynetat0: Lecture echouee pour <cv>" 216 CALL abort216 CALL ABORT_gcm("dynetat0", "", 1) 217 217 ENDIF 218 218 … … 220 220 IF (ierr .NE. NF_NOERR) THEN 221 221 write(lunout,*)"dynetat0: Le champ <aire> est absent" 222 CALL abort222 CALL ABORT_gcm("dynetat0", "", 1) 223 223 ENDIF 224 224 ierr = nf90_get_var(nid, nvarid, aire) 225 225 IF (ierr .NE. NF_NOERR) THEN 226 226 write(lunout,*)"dynetat0: Lecture echouee pour <aire>" 227 CALL abort227 CALL ABORT_gcm("dynetat0", "", 1) 228 228 ENDIF 229 229 … … 231 231 IF (ierr .NE. NF_NOERR) THEN 232 232 write(lunout,*)"dynetat0: Le champ <phisinit> est absent" 233 CALL abort233 CALL ABORT_gcm("dynetat0", "", 1) 234 234 ENDIF 235 235 ierr = nf90_get_var(nid, nvarid, phis) 236 236 IF (ierr .NE. NF_NOERR) THEN 237 237 write(lunout,*)"dynetat0: Lecture echouee pour <phisinit>" 238 CALL abort238 CALL ABORT_gcm("dynetat0", "", 1) 239 239 ENDIF 240 240 … … 246 246 IF (ierr .NE. NF_NOERR) THEN 247 247 write(lunout,*)"dynetat0: Le champ <Time> est absent" 248 CALL abort248 CALL ABORT_gcm("dynetat0", "", 1) 249 249 ENDIF 250 250 ENDIF … … 252 252 IF (ierr .NE. NF_NOERR) THEN 253 253 write(lunout,*)"dynetat0: Lecture echouee <temps>" 254 CALL abort254 CALL ABORT_gcm("dynetat0", "", 1) 255 255 ENDIF 256 256 … … 258 258 IF (ierr .NE. NF_NOERR) THEN 259 259 write(lunout,*)"dynetat0: Le champ <ucov> est absent" 260 CALL abort260 CALL ABORT_gcm("dynetat0", "", 1) 261 261 ENDIF 262 262 ierr = nf90_get_var(nid, nvarid, ucov) 263 263 IF (ierr .NE. NF_NOERR) THEN 264 264 write(lunout,*)"dynetat0: Lecture echouee pour <ucov>" 265 CALL abort265 CALL ABORT_gcm("dynetat0", "", 1) 266 266 ENDIF 267 267 … … 269 269 IF (ierr .NE. NF_NOERR) THEN 270 270 write(lunout,*)"dynetat0: Le champ <vcov> est absent" 271 CALL abort271 CALL ABORT_gcm("dynetat0", "", 1) 272 272 ENDIF 273 273 ierr = nf90_get_var(nid, nvarid, vcov) 274 274 IF (ierr .NE. NF_NOERR) THEN 275 275 write(lunout,*)"dynetat0: Lecture echouee pour <vcov>" 276 CALL abort276 CALL ABORT_gcm("dynetat0", "", 1) 277 277 ENDIF 278 278 … … 280 280 IF (ierr .NE. NF_NOERR) THEN 281 281 write(lunout,*)"dynetat0: Le champ <teta> est absent" 282 CALL abort282 CALL ABORT_gcm("dynetat0", "", 1) 283 283 ENDIF 284 284 ierr = nf90_get_var(nid, nvarid, teta) 285 285 IF (ierr .NE. NF_NOERR) THEN 286 286 write(lunout,*)"dynetat0: Lecture echouee pour <teta>" 287 CALL abort287 CALL ABORT_gcm("dynetat0", "", 1) 288 288 ENDIF 289 289 … … 301 301 IF (ierr .NE. NF_NOERR) THEN 302 302 write(lunout,*)"dynetat0: Lecture echouee pour "//tname(iq) 303 CALL abort303 CALL ABORT_gcm("dynetat0", "", 1) 304 304 ENDIF 305 305 ENDIF … … 310 310 IF (ierr .NE. NF_NOERR) THEN 311 311 write(lunout,*)"dynetat0: Le champ <masse> est absent" 312 CALL abort312 CALL ABORT_gcm("dynetat0", "", 1) 313 313 ENDIF 314 314 ierr = nf90_get_var(nid, nvarid, masse) 315 315 IF (ierr .NE. NF_NOERR) THEN 316 316 write(lunout,*)"dynetat0: Lecture echouee pour <masse>" 317 CALL abort317 CALL ABORT_gcm("dynetat0", "", 1) 318 318 ENDIF 319 319 … … 321 321 IF (ierr .NE. NF_NOERR) THEN 322 322 write(lunout,*)"dynetat0: Le champ <ps> est absent" 323 CALL abort323 CALL ABORT_gcm("dynetat0", "", 1) 324 324 ENDIF 325 325 ierr = nf90_get_var(nid, nvarid, ps) 326 326 IF (ierr .NE. NF_NOERR) THEN 327 327 write(lunout,*)"dynetat0: Lecture echouee pour <ps>" 328 CALL abort328 CALL ABORT_gcm("dynetat0", "", 1) 329 329 ENDIF 330 330 -
LMDZ5/trunk/libf/dyn3d/dynredem.F
r1907 r1930 133 133 & //trim(fichnom) 134 134 write(lunout,*)' ierr = ', ierr 135 CALL ABORT 135 CALL ABORT_GCM("DYNREDEM0", "", 1) 136 136 ENDIF 137 137 c … … 512 512 IF (ierr .NE. NF_NOERR) THEN 513 513 write(lunout,*)"dynredem1: Pb. d ouverture "//trim(fichnom) 514 CALL abort514 call abort_gcm("dynredem1", "", 1) 515 515 ENDIF 516 516 -
LMDZ5/trunk/libf/dyn3d/fxhyp.F
r1907 r1930 178 178 WRITE(6,*)'Modifier les valeurs de grossismx ,tau ou dzoomx ', 179 179 , ' et relancer ! *** ' 180 CALL ABORT 180 CALL ABORT_GCM("FXHYP", "", 1) 181 181 ENDIF 182 182 c -
LMDZ5/trunk/libf/dyn3d/fyhyp.F
r1907 r1930 162 162 WRITE(6,*)'Modifier les valeurs de grossismy ,tauy ou dzoomy', 163 163 , ' et relancer ! *** ' 164 CALL ABORT 164 CALL ABORT_GCM("FYHYP", "", 1) 165 165 166 166 ENDIF -
LMDZ5/trunk/libf/dyn3d/gcm.F
r1907 r1930 327 327 start_time = starttime 328 328 ELSE 329 WRITE(lunout,*)'Je m''arrete' 330 CALL abort 329 call abort_gcm("gcm", "'Je m''arrete'", 1) 331 330 ENDIF 332 331 ENDIF -
LMDZ5/trunk/libf/dyn3d/grid_atob.F
r1907 r1930 55 55 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN 56 56 PRINT*, 'imar ou jmar trop grand', imar, jmar 57 CALL ABORT 57 CALL ABORT_GCM("", "", 1) 58 58 ENDIF 59 59 c … … 119 119 ELSE 120 120 PRINT*, 'probleme,i,j=', i,j 121 ccc CALL ABORT 121 ccc CALL ABORT_GCM("", "", 1) 122 122 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) 123 123 #ifdef CRAY … … 193 193 IF (imar.GT.400 .OR. jmar.GT.200) THEN 194 194 PRINT*, 'imar ou jmar trop grand', imar, jmar 195 CALL ABORT 195 CALL ABORT_GCM("", "", 1) 196 196 ENDIF 197 197 c 198 198 IF (imdep.GT.400 .OR. jmdep.GT.200) THEN 199 199 PRINT*, 'imdep ou jmdep trop grand', imdep, jmdep 200 CALL ABORT 200 CALL ABORT_GCM("", "", 1) 201 201 ENDIF 202 202 c … … 258 258 PRINT*, 'Probleme grave,i,j,indx,indy=', 259 259 . i,j,indx(i,j),indy(i,j) 260 CALL abort260 call abort_gcm("", "", 1) 261 261 ENDIF 262 262 ENDDO … … 309 309 ELSE 310 310 PRINT*, 'probleme,i,j=', i,j 311 CALL ABORT 311 CALL ABORT_GCM("", "", 1) 312 312 ENDIF 313 313 ENDDO … … 345 345 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN 346 346 PRINT*, 'imar ou jmar trop grand', imar, jmar 347 CALL ABORT 347 CALL ABORT_GCM("", "", 1) 348 348 ENDIF 349 349 c … … 409 409 ELSE 410 410 PRINT*, 'probleme,i,j=', i,j 411 CALL ABORT 411 CALL ABORT_GCM("", "", 1) 412 412 ENDIF 413 413 ENDDO … … 452 452 IF (imar.GT.400 .OR. jmar.GT.400) THEN 453 453 PRINT*, 'imar ou jmar trop grand', imar, jmar 454 CALL ABORT 454 CALL ABORT_GCM("", "", 1) 455 455 ENDIF 456 456 c … … 513 513 ELSE 514 514 PRINT*, 'probleme,i,j=', i,j 515 ccc CALL ABORT 515 ccc CALL ABORT_GCM("", "", 1) 516 516 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) 517 517 #ifdef CRAY … … 577 577 IF (imar.GT.400 .OR. jmar.GT.400) THEN 578 578 PRINT*, 'imar ou jmar trop grand', imar, jmar 579 CALL ABORT 579 CALL ABORT_GCM("", "", 1) 580 580 ENDIF 581 581 c … … 642 642 ELSE 643 643 PRINT*, 'probleme,i,j=', i,j 644 ccc CALL ABORT 644 ccc CALL ABORT_GCM("", "", 1) 645 645 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) 646 646 #ifdef CRAY … … 712 712 IF (immod.GT.2200 .OR. jmmod.GT.1100) THEN 713 713 PRINT*, 'immod ou jmmod trop grand', immod, jmmod 714 CALL ABORT 714 CALL ABORT_GCM("", "", 1) 715 715 ENDIF 716 716 c … … 801 801 ELSE 802 802 PRINT*, 'probleme,i,j=', i,j 803 CALL ABORT 803 CALL ABORT_GCM("", "", 1) 804 804 ENDIF 805 805 ENDDO … … 875 875 ELSE 876 876 PRINT*, 'probleme,i,j=', i,j 877 ccc CALL ABORT 877 ccc CALL ABORT_GCM("", "", 1) 878 878 CALL dist_sphe(xmod(i),ymod(j),xtmp,ytmp,imtmp,jmtmp,distans) 879 879 #ifdef CRAY -
LMDZ5/trunk/libf/dyn3d/grid_noro.F
r1907 r1930 111 111 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN 112 112 PRINT*, 'imar or jmar too big', imar, jmar 113 CALL ABORT 113 CALL ABORT_GCM("GRID_NORO", "", 1) 114 114 ENDIF 115 115 116 116 IF(imdep.ne.iusn.or.jmdep.ne.jusn)then 117 117 print *,' imdep or jmdep bad dimensions:',imdep,jmdep 118 call abort 118 call abort_gcm("grid_noro", "", 1) 119 119 ENDIF 120 120 121 121 IF(imar+1.ne.iim+1.or.jmar.ne.jjm+1)THEN 122 122 print *,' imar or jmar bad dimensions:',imar,jmar 123 call abort 123 call abort_gcm("grid_noro", "", 1) 124 124 ENDIF 125 125 -
LMDZ5/trunk/libf/dyn3d/juldate.F
r1907 r1930 9 9 implicit real (a-h,o-z) 10 10 frac=((os/60.+om)/60.+oh)/24. 11 ojou=d float(ijou)+frac12 year=d float(ian)13 rmon=d float(imoi)11 ojou=dble(ijou)+frac 12 year=dble(ian) 13 rmon=dble(imoi) 14 14 if (imoi .le. 2) then 15 15 year=year-1. -
LMDZ5/trunk/libf/dyn3d/ran1.F
r1907 r1930 28 28 IX3=MOD(IA3*IX3+IC3,M3) 29 29 J=1+(97*IX3)/M3 30 IF(J.GT.97.OR.J.LT.1) PAUSE30 IF(J.GT.97.OR.J.LT.1) stop 1 31 31 RAN1=R(J) 32 32 R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
Note: See TracChangeset
for help on using the changeset viewer.