Changeset 3893 for trunk/LMDZ.GENERIC/libf/dynphy_lonlat
- Timestamp:
- Aug 13, 2025, 7:53:49 PM (4 months ago)
- Location:
- trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd
- Files:
-
- 8 edited
-
datareadnc.F (modified) (1 diff)
-
defrun_new.F (modified) (16 diffs)
-
grid_noro1.F (modified) (1 diff)
-
inichim_newstart.F90 (modified) (4 diffs)
-
lect_start_archive.F (modified) (33 diffs)
-
newstart.F (modified) (12 diffs)
-
start2archive.F (modified) (7 diffs)
-
write_archive.F (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/datareadnc.F
r3713 r3893 127 127 write(*,*)' http://www.lmd.jussieu.fr/~lmdz/planets/'// 128 128 & 'generic/datagcm/' 129 STOP 129 call abort_physic("datareadnc", 130 & "cannot open file surface file",1) 130 131 ENDIF 131 132 -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/defrun_new.F
r2354 r3893 249 249 250 250 IF( ABS(clon - clonn).GE. 0.001 ) THEN 251 PRINT *,' La valeur de clon passee par run.def est ' 252 * ,'differente de celle lue sur le fichier start ' 253 STOP 251 call abort_physic("defrun_new", 252 & "Value mismatch of clon in run.def and start file",1) 254 253 ENDIF 255 254 c … … 259 258 260 259 IF( ABS(clat - clatt).GE. 0.001 ) THEN 261 PRINT *,' La valeur de clat passee par run.def est ' 262 * ,'differente de celle lue sur le fichier start ' 263 STOP 260 call abort_physic("defrun_new", 261 & "Value mismatch of clat in run.def and start file",1) 264 262 ENDIF 265 263 … … 275 273 write(*,*)'place de alphax, alphay. cf. dyn3d). ' 276 274 write(*,*) 277 stop 275 call abort_physic("defrun_new", 276 & "ERREUR : dans run.def, grossismx =0",1) 278 277 end if 279 278 280 279 IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN 281 280 PRINT *,' La valeur de grossismx passee par run.def est ' 282 * ,'differente de celle lue sur le fichier 281 * ,'differente de celle lue sur le fichier start =', 283 282 * grossismx 284 283 if (grossismx.eq.0) then … … 287 286 grossismx=grossismxx 288 287 else 289 STOP 288 call abort_physic("defrun_new", 289 & "Value mismatch of grossismx in run.def and start file",1) 290 290 endif 291 291 ENDIF … … 303 303 grossismy=grossismyy 304 304 else 305 STOP 305 call abort_physic("defrun_new", 306 & "Value mismatch of grossismy in run.def and start file",1) 306 307 endif 307 308 ENDIF … … 309 310 310 311 IF( grossismx.LT.1. ) THEN 311 PRINT *,' *** ATTENTION !! grossismx < 1 . *** '312 STOP312 call abort_physic("defrun_new", 313 & "ATTENTION !! grossismx < 1",1) 313 314 ELSE 314 315 alphax = 1. - 1./ grossismx … … 316 317 317 318 IF( grossismy.LT.1. ) THEN 318 PRINT *,' *** ATTENTION !! grossismy < 1 . *** '319 STOP319 call abort_physic("defrun_new", 320 & "ATTENTION !! grossismy < 1",1) 320 321 ELSE 321 322 alphay = 1. - 1./ grossismy … … 334 335 PRINT *,' *** fxyhypb lu sur le fichier start est F ', 335 336 * 'alors qu il est T sur run.def ***' 336 STOP 337 call abort_physic("defrun_new", 338 & "fxyhypb=F dans start alors que =T dans run.def",1) 337 339 ENDIF 338 340 ELSE … … 341 343 PRINT *,' *** fxyhypb lu sur le fichier start est T ', 342 344 * 'alors qu il est F sur run.def **** ' 343 STOP 345 call abort_physic("defrun_new", 346 & "fxyhypb=T dans start alors que =F dans run.def",1) 344 347 ENDIF 345 348 ENDIF … … 348 351 349 352 IF( fxyhypb ) THEN 350 IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN 351 PRINT *,' La valeur de dzoomx passee par run.def est ' 352 * ,'differente de celle lue sur le fichier start ' 353 STOP 354 ENDIF 353 IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN 354 call abort_physic("defrun_new", 355 & "Value mismatch of dzoomx between run.def and startfile",1) 356 ENDIF 355 357 ENDIF 356 358 … … 359 361 360 362 IF( fxyhypb ) THEN 361 IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN 362 PRINT *,' La valeur de dzoomy passee par run.def est ' 363 * ,'differente de celle lue sur le fichier start ' 364 STOP 365 ENDIF 363 IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN 364 call abort_physic("defrun_new", 365 & "Value mismatch of dzoomy between run.def and startfile",1) 366 ENDIF 366 367 ENDIF 367 368 … … 374 375 IF( fxyhypb ) THEN 375 376 IF( ABS(taux - tauxx).GE. 0.001 ) THEN 376 377 WRITE(6,*)' La valeur de taux passee par run.def est', 377 378 * 'differente de celle lue sur le fichier start ' 378 CALL ABORT 379 call abort_physic("defrun_new", 380 & "Value mismatch of taux between run.def and startfile",1) 379 381 ENDIF 380 382 381 383 IF( ABS(tauy - tauyy).GE. 0.001 ) THEN 382 384 WRITE(6,*)' La valeur de tauy passee par run.def est', 383 385 * 'differente de celle lue sur le fichier start ' 384 CALL ABORT 386 call abort_physic("defrun_new", 387 & "Value mismatch of tauy between run.def and startfile",1) 385 388 ENDIF 386 389 ENDIF … … 416 419 417 420 IF( grossismx.LT.1. ) THEN 418 PRINT *,' *** ATTENTION !! grossismx < 1 . *** '419 STOP421 call abort_physic("defrun_new", 422 & "ATTENTION !! grossismx < 1",1) 420 423 ELSE 421 424 alphax = 1. - 1./ grossismx … … 423 426 424 427 IF( grossismy.LT.1. ) THEN 425 PRINT *,' *** ATTENTION !! grossismy < 1 . *** '426 STOP428 call abort_physic("defrun_new", 429 & "ATTENTION !! grossismy < 1",1) 427 430 ELSE 428 431 alphay = 1. - 1./ grossismy … … 513 516 ELSE 514 517 write(tapeerr,*) ' WHERE IS run.def ? WE NEED IT !!!!!!!!!!!!!!' 515 stop518 call abort_physic("defrun_new","Unable to find run.def",1) 516 519 ENDIF ! of IF(ierr.eq.0) 517 520 … … 525 528 else 526 529 c Avec Zoom 527 if (.not.fxyhypb) stop 'With zoom, fxyhypb should be set to T528 & in run.def for this version... -> STOP ! '530 if (.not.fxyhypb) call abort_physic("defrun_new", 531 & "With zoom, fxyhypb should be T in run.def for this version",1) 529 532 end if 530 533 #else 531 write(*,*) "defrun_new should not be used in parallel mode!"532 stop534 call abort_physic("defrun_new", 535 & "defrun_new should not be used in parallel mode",1) 533 536 #endif 534 537 ! of #ifndef CPP_PARA -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/grid_noro1.F
r1422 r3893 109 109 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN 110 110 PRINT*, 'imar ou jmar trop grand', imar, jmar 111 CALL ABORT111 CALL abort_physic("grid_noro1.F","imar ou jmar trop grand",1) 112 112 ENDIF 113 113 114 114 IF(imdep.ne.iusn.or.jmdep.ne.jusn)then 115 115 print *,' imdep ou jmdep mal dimensionnes:',imdep,jmdep 116 call abort 116 call abort_physic("grid_noro1.F", 117 & "imdep ou jmdep mal dimensionnes",1) 117 118 ENDIF 118 119 119 120 IF(imar+1.gt.iim+1.or.jmar.gt.jjm+1)THEN 120 121 print *,' imar ou jmar mal dimensionnes:',imar,jmar 121 call abort 122 call abort_physic("grid_noro1.F", 123 & "imar ou jmar mal dimensionnes",1) 122 124 ENDIF 123 125 c -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/inichim_newstart.F90
r3073 r3893 94 94 write(*,*) "calchim: error reading line of tracers" 95 95 write(*,*) " (first lines of traceur.def) " 96 stop96 call abort_physic("traceur.def","error reading line of tracers",1) 97 97 ENDIF 98 98 ENDDO 99 99 ENDIF ! if modern or standard traceur.def 100 100 else 101 write(*,*) "calchim: error opening traceur.def in inichim_newstart" 102 stop 101 call abort_physic("inichim_newstart","Unable to open traceur.def",1) 103 102 endif 104 103 … … 165 164 else 166 165 write(*,*) 'inichim_newstart: error opening ', fil 167 stop166 call abort_physic("inichim_newstart", "Unable to open chemimal profile file",1) 168 167 endif 169 168 close(90) … … 181 180 pqx(ilon,ilat,ilay,iq) = 1-sum(pqx(ilon,ilat,ilay,:)) 182 181 if (pqx(ilon,ilat,ilay,iq)<=0.) then 183 write(*,*) 'inichim_newstart: vmr tot > 1 not possible' 184 stop 182 call abort_physic("inichim_newstart","vmr tot > 1 not possible",1) 185 183 end if 186 184 end do … … 192 190 end do 193 191 if (.not.foundback) then 194 write(*,*) 'inichim_newstart: you need to set a background gas' 195 write(*,*) ' by qx=1. in traceur.def' 196 stop 192 call abort_physic("inichim_newstart","you need to set a background gas by qx=1. in traceur.def",1) 197 193 end if 198 194 -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F
r3423 r3893 319 319 ierr = NF_INQ_VARID (nid, "controle", nvarid) 320 320 IF (ierr .NE. NF_NOERR) THEN 321 PRINT*, "Lect_start_archive: champ <controle> not found"322 CALL abort321 CALL abort_physic("lect_start_archive", 322 & "Field <controle> not found",1) 323 323 ENDIF 324 324 #ifdef NC_DOUBLE … … 328 328 #endif 329 329 IF (ierr .NE. NF_NOERR) THEN 330 PRINT*, "lect_start_archive: Lecture echoue pour <controle>"331 CALL abort330 CALL abort_physic("lect_start_archive", 331 & "Lecture echoue pour <controle>",1) 332 332 ENDIF 333 333 c … … 340 340 ierr = NF_INQ_VARID (nid, "rlonv", nvarid) 341 341 IF (ierr .NE. NF_NOERR) THEN 342 PRINT*, "lect_start_archive: Field <rlonv> not found"343 CALL abort342 CALL abort_physic("lect_start_archive", 343 & "Field <rlonv> not found",1) 344 344 ENDIF 345 345 #ifdef NC_DOUBLE … … 349 349 #endif 350 350 IF (ierr .NE. NF_NOERR) THEN 351 PRINT*, "lect_start_archive: Failed loading <rlonv>"352 CALL abort351 CALL abort_physic("lect_start_archive", 352 & "Failed loading <rlonv>",1) 353 353 ENDIF 354 354 c 355 355 ierr = NF_INQ_VARID (nid, "rlatu", nvarid) 356 356 IF (ierr .NE. NF_NOERR) THEN 357 PRINT*, "lect_start_archive: Field <rlatu> not found"358 CALL abort357 CALL abort_physic("lect_start_archive", 358 & "Field <rlatu> not found",1) 359 359 ENDIF 360 360 #ifdef NC_DOUBLE … … 364 364 #endif 365 365 IF (ierr .NE. NF_NOERR) THEN 366 PRINT*, "lect_start_archive: Failed loading <rlatu>"367 CALL abort366 CALL abort_physic("lect_start_archive", 367 & "Failed loading <rlatu>",1) 368 368 ENDIF 369 369 c 370 370 ierr = NF_INQ_VARID (nid, "rlonu", nvarid) 371 371 IF (ierr .NE. NF_NOERR) THEN 372 PRINT*, "lect_start_archive: Field <rlonu> not found"373 CALL abort372 CALL abort_physic("lect_start_archive", 373 & "Field <rlonu> not found",1) 374 374 ENDIF 375 375 #ifdef NC_DOUBLE … … 379 379 #endif 380 380 IF (ierr .NE. NF_NOERR) THEN 381 PRINT*, "lect_start_archive: Failed loading <rlonu>"382 CALL abort381 CALL abort_physic("lect_start_archive", 382 & "Failed loading <rlonu>",1) 383 383 ENDIF 384 384 c 385 385 ierr = NF_INQ_VARID (nid, "rlatv", nvarid) 386 386 IF (ierr .NE. NF_NOERR) THEN 387 PRINT*, "lect_start_archive: Field <rlatv> not found"388 CALL abort387 CALL abort_physic("lect_start_archive", 388 & "Field <rlatv> not found",1) 389 389 ENDIF 390 390 #ifdef NC_DOUBLE … … 394 394 #endif 395 395 IF (ierr .NE. NF_NOERR) THEN 396 PRINT*, "lect_start_archive: Failed loading <rlatv>"397 CALL abort396 CALL abort_physic("lect_start_archive", 397 & "Failed loading <rlatv>",1) 398 398 ENDIF 399 399 c … … 426 426 IF (ierr .NE. NF_NOERR) THEN 427 427 PRINT*, "Nothing to do..." 428 CALL abort 428 CALL abort_physic("lect_start_archive", 429 & "Field <bps> or <sig_s> not found",1) 429 430 ENDIF 430 431 ENDIF … … 435 436 #endif 436 437 IF (ierr .NE. NF_NOERR) THEN 437 PRINT*, "lect_start_archive: Failed loading <bps>"438 CALL abort438 CALL abort_physic("lect_start_archive", 439 & "Failed loading <bps>",1) 439 440 END IF 440 441 … … 460 461 #endif 461 462 if (ierr .NE. NF_NOERR) then 462 PRINT*, "lect_start_archive: Failed reading <soildepth>"463 CALL abort463 CALL abort_physic("lect_start_archive", 464 & "Failed reading <soildepth>",1) 464 465 endif 465 466 … … 488 489 #endif 489 490 if (ierr .NE. NF_NOERR) then 490 PRINT*, "lect_start_archive: Failed reading <soildepth>"491 CALL abort491 CALL abort_physic("lect_start_archive", 492 & "Failed reading <soildepth>",1) 492 493 endif 493 494 endif ! of if (depthinterpol) … … 513 514 #endif 514 515 if (ierr .NE. NF_NOERR) then 515 PRINT*, "lect_start_archive: Failed reading <inertiedat>"516 CALL abort516 CALL abort_physic("lect_start_archive", 517 & "Failed reading <inertiedat>",1) 517 518 endif 518 519 endif … … 524 525 ierr = NF_INQ_VARID (nid, "phisinit", nvarid) 525 526 IF (ierr .NE. NF_NOERR) THEN 526 PRINT*, "lect_start_archive: Field <phisinit> not found"527 CALL abort527 CALL abort_physic("lect_start_archive", 528 & "Field <phisinit> not found",1) 528 529 ENDIF 529 530 #ifdef NC_DOUBLE … … 533 534 #endif 534 535 IF (ierr .NE. NF_NOERR) THEN 535 PRINT*, "lect_start_archive: Failed loading <phisinit>"536 CALL abort536 CALL abort_physic("lect_start_archive", 537 & "Failed loading <phisinit>",1) 537 538 ENDIF 538 539 … … 553 554 ierr = NF_INQ_DIMID (nid, "temps", nvarid) 554 555 IF (ierr .NE. NF_NOERR) THEN 555 PRINT*, "lect_start_archive: Field <Time> not found"556 CALL abort556 CALL abort_physic("lect_start_archive", 557 & "Field <Time> not found",1) 557 558 endif 558 559 ENDIF … … 568 569 #endif 569 570 IF (ierr .NE. NF_NOERR) THEN 570 PRINT*, "lect_start_archive: Failed loading <Time>"571 CALL abort571 CALL abort_physic("lect_start_archive", 572 & "Failed loading <Time>",1) 572 573 ENDIF 573 574 c … … 626 627 ierr = NF_INQ_VARID (nid, "emis", nvarid) 627 628 IF (ierr .NE. NF_NOERR) THEN 628 PRINT*, "lect_start_archive: Field <emis> not found"629 CALL abort629 CALL abort_physic("lect_start_archive", 630 & "Field <emis> not found",1) 630 631 ENDIF 631 632 #ifdef NC_DOUBLE … … 635 636 #endif 636 637 IF (ierr .NE. NF_NOERR) THEN 637 PRINT*, "lect_start_archive: Failed loading <emis>"638 CALL abort638 CALL abort_physic("lect_start_archive", 639 & "Failed loading <emis>",1) 639 640 ENDIF 640 641 c 641 642 ierr = NF_INQ_VARID (nid, "ps", nvarid) 642 643 IF (ierr .NE. NF_NOERR) THEN 643 PRINT*, "lect_start_archive: Field <ps> not found"644 CALL abort644 CALL abort_physic("lect_start_archive", 645 & "Field <ps> not found",1) 645 646 ENDIF 646 647 #ifdef NC_DOUBLE … … 650 651 #endif 651 652 IF (ierr .NE. NF_NOERR) THEN 652 PRINT*, "lect_start_archive: Failed loading <ps>"653 CALL abort653 CALL abort_physic("lect_start_archive", 654 & "Failed loading <ps>",1) 654 655 ENDIF 655 656 c 656 657 ierr = NF_INQ_VARID (nid, "tsurf", nvarid) 657 658 IF (ierr .NE. NF_NOERR) THEN 658 PRINT*, "lect_start_archive: Field <tsurf> not found"659 CALL abort659 CALL abort_physic("lect_start_archive", 660 & "Field <tsurf> not found",1) 660 661 ENDIF 661 662 #ifdef NC_DOUBLE … … 665 666 #endif 666 667 IF (ierr .NE. NF_NOERR) THEN 667 PRINT*, "lect_start_archive: Failed loading <tsurf>"668 CALL abort668 CALL abort_physic("lect_start_archive", 669 & "Failed loading <tsurf>",1) 669 670 ENDIF 670 671 c 671 672 ierr = NF_INQ_VARID (nid, "q2surf", nvarid) 672 673 IF (ierr .NE. NF_NOERR) THEN 673 PRINT*, "lect_start_archive: Field <q2surf> not found"674 CALL abort674 CALL abort_physic("lect_start_archive", 675 & "Field <q2surf> not found",1) 675 676 ENDIF 676 677 #ifdef NC_DOUBLE … … 680 681 #endif 681 682 IF (ierr .NE. NF_NOERR) THEN 682 PRINT*, "lect_start_archive: Failed loading <q2surf>"683 CALL abort683 CALL abort_physic("lect_start_archive", 684 & "Failed loading <q2surf>",1) 684 685 ENDIF 685 686 c … … 837 838 PRINT*, "lect_start_archive: ", 838 839 & "Field <","Tg"//str2,"> not found" 839 CALL abort 840 CALL abort_physic("lect_start_archive", 841 & "Field <Tg> not found",1) 840 842 ENDIF 841 843 #ifdef NC_DOUBLE … … 849 851 PRINT*, "lect_start_archive: ", 850 852 & "Failed reading <","Tg"//str2,">" 851 CALL abort 853 CALL abort_physic("lect_start_archive", 854 & "Failed reading <Tg>",1) 852 855 ENDIF 853 856 c … … 862 865 ierr=NF_INQ_VARID(nid,"tsoil",nvarid) 863 866 if (ierr.ne.NF_NOERR) then 864 write(*,*)"lect_start_archive: Cannot find <tsoil>" 865 call abort 867 CALL abort_physic("lect_start_archive", 868 & "Cannot find <tsoil>",1) 866 869 else 867 870 #ifdef NC_DOUBLE … … 884 887 ierr = NF_INQ_VARID (nid,"temp", nvarid) 885 888 IF (ierr .NE. NF_NOERR) THEN 886 PRINT*, "lect_start_archive: Field <temp> not found"887 CALL abort889 CALL abort_physic("lect_start_archive", 890 & "Field <temp> not found",1) 888 891 ENDIF 889 892 #ifdef NC_DOUBLE … … 893 896 #endif 894 897 IF (ierr .NE. NF_NOERR) THEN 895 PRINT*, "lect_start_archive: Failed loading <temp>"896 CALL abort898 CALL abort_physic("lect_start_archive", 899 & "Failed loading <temp>",1) 897 900 ENDIF 898 901 c 899 902 ierr = NF_INQ_VARID (nid,"u", nvarid) 900 903 IF (ierr .NE. NF_NOERR) THEN 901 PRINT*, "lect_start_archive: Field <u> not found"902 CALL abort904 CALL abort_physic("lect_start_archive", 905 & "Field <u> not found",1) 903 906 ENDIF 904 907 #ifdef NC_DOUBLE … … 908 911 #endif 909 912 IF (ierr .NE. NF_NOERR) THEN 910 PRINT*, "lect_start_archive: Failed loading <u>"911 CALL abort913 CALL abort_physic("lect_start_archive", 914 & "Failed loading <u>",1) 912 915 ENDIF 913 916 c 914 917 ierr = NF_INQ_VARID (nid,"v", nvarid) 915 918 IF (ierr .NE. NF_NOERR) THEN 916 PRINT*, "lect_start_archive: Field <v> not found"917 CALL abort919 CALL abort_physic("lect_start_archive", 920 & "Field <v> not found",1) 918 921 ENDIF 919 922 #ifdef NC_DOUBLE … … 923 926 #endif 924 927 IF (ierr .NE. NF_NOERR) THEN 925 PRINT*, "lect_start_archive: Failed loading <v>"926 CALL abort928 CALL abort_physic("lect_start_archive", 929 & "Failed loading <v>",1) 927 930 ENDIF 928 931 c 929 932 ierr = NF_INQ_VARID (nid,"q2atm", nvarid) 930 933 IF (ierr .NE. NF_NOERR) THEN 931 PRINT*, "lect_start_archive: Field <q2atm> not found"932 CALL abort934 CALL abort_physic("lect_start_archive", 935 & "Field <q2atm> not found",1) 933 936 ENDIF 934 937 #ifdef NC_DOUBLE … … 938 941 #endif 939 942 IF (ierr .NE. NF_NOERR) THEN 940 PRINT*, "lect_start_archive: Failed loading <q2atm>"941 CALL abort943 CALL abort_physic("lect_start_archive", 944 & "Failed loading <q2atm>",1) 942 945 ENDIF 943 946 c … … 989 992 #endif 990 993 IF (ierr .NE. NF_NOERR) THEN 991 PRINT*, "lect_start_archive: Failed loading <du_nonoro_gwd>"992 CALL abort994 CALL abort_physic("lect_start_archive", 995 & "Failed loading <du_nonoro_gwd>",1) 993 996 ENDIF 994 997 ENDIF … … 1007 1010 #endif 1008 1011 IF (ierr .NE. NF_NOERR) THEN 1009 PRINT*, "lect_start_archive: Failed loading <dv_nonoro_gwd>"1010 CALL abort1012 CALL abort_physic("lect_start_archive", 1013 & "Failed loading <dv_nonoro_gwd>",1) 1011 1014 ENDIF 1012 1015 ENDIF … … 1025 1028 #endif 1026 1029 IF (ierr .NE. NF_NOERR) THEN 1027 PRINT*, "lect_start_archive: Failed loading <east_gwstress>"1028 CALL abort1030 CALL abort_physic("lect_start_archive", 1031 & "Failed loading <east_gwstress>",1) 1029 1032 ENDIF 1030 1033 ENDIF … … 1043 1046 #endif 1044 1047 IF (ierr .NE. NF_NOERR) THEN 1045 PRINT*, "lect_start_archive: Failed loading <west_gwstress>"1046 CALL abort1048 CALL abort_physic("lect_start_archive", 1049 & "Failed loading <west_gwstress>",1) 1047 1050 ENDIF 1048 1051 ENDIF -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F
r3836 r3893 259 259 write(6,*)' Problem opening file:',fichnom 260 260 write(6,*)' ierr = ', ierr 261 CALL ABORT 261 CALL abort_physic("newstart", 262 & "Problem opening start_archive.nc file",1) 262 263 ENDIF 263 264 tab0 = 50 … … 275 276 write(6,*)' Problem opening file:',fichnom 276 277 write(6,*)' ierr = ', ierr 277 CALL ABORT 278 CALL abort_physic("newstart", 279 & "Problem opening start.nc file", 1) 278 280 ENDIF 279 281 … … 283 285 write(6,*)' Problem opening file:',fichnom 284 286 write(6,*)' ierr = ', ierr 285 CALL ABORT 287 CALL abort_physic("newstart", 288 & "Problem opening startfi.nc file", 1) 286 289 ENDIF 287 290 … … 1095 1098 ! check that there is indeed a water vapor tracer 1096 1099 if (igcm_h2o_vap.eq.0) then 1097 write(*,*) "No water vapour tracer! Can't use this option"1098 stop1100 call abort_physic("newstart", 1101 & "No water vapour tracer! Can't use wetstart",1) 1099 1102 endif 1100 1103 DO l=1,llm … … 1115 1118 else if (trim(modif) .eq. 'noglacier') then 1116 1119 if (igcm_h2o_ice.eq.0) then 1117 write(*,*) "No water ice tracer! Can't use this option"1118 stop1120 call abort_physic("newstart", 1121 & "No water ice tracer! Can't use noglacier",1) 1119 1122 endif 1120 1123 do ig=1,ngridmx … … 1132 1135 else if (trim(modif) .eq. 'watercapn') then 1133 1136 if (igcm_h2o_ice.eq.0) then 1134 write(*,*) "No water ice tracer! Can't use this option"1135 stop1137 call abort_physic("newstart", 1138 & "No water ice tracer! Can't use watercapn",1) 1136 1139 endif 1137 1140 … … 1158 1161 else if (trim(modif) .eq. 'watercaps') then 1159 1162 if (igcm_h2o_ice.eq.0) then 1160 write(*,*) "No water ice tracer! Can't use this option"1161 stop1163 call abort_physic("newstart", 1164 & "No water ice tracer! Can't use watercaps",1) 1162 1165 endif 1163 1166 … … 1185 1188 else if (trim(modif) .eq. 'noacglac') then 1186 1189 if (igcm_h2o_ice.eq.0) then 1187 write(*,*) "No water ice tracer! Can't use this option"1188 stop1190 call abort_physic("newstart", 1191 & "No water ice tracer! Can't use noacglac",1) 1189 1192 endif 1190 1193 DO j=1,jjp1 … … 1214 1217 else if (trim(modif) .eq. 'oborealis') then 1215 1218 if (igcm_h2o_ice.eq.0) then 1216 write(*,*) "No water ice tracer! Can't use this option"1217 stop1219 call abort_physic("newstart", 1220 & "No water ice tracer! Can't use oborealis",1) 1218 1221 endif 1219 1222 DO j=1,jjp1 … … 1246 1249 else if (trim(modif) .eq. 'iborealis') then 1247 1250 if (igcm_h2o_ice.eq.0) then 1248 write(*,*) "No water ice tracer! Can't use this option"1249 stop1251 call abort_physic("newstart", 1252 & "No water ice tracer! Can't use iborealis",1) 1250 1253 endif 1251 1254 DO j=1,jjp1 … … 1270 1273 else if (trim(modif) .eq. 'oceanball') then 1271 1274 if (igcm_h2o_ice.eq.0) then 1272 write(*,*) "No water ice tracer! Can't use this option"1273 stop1275 call abort_physic("newstart", 1276 & "No water ice tracer! Can't use oceanball",1) 1274 1277 endif 1275 1278 DO j=1,jjp1 … … 1299 1302 else if (trim(modif) .eq. 'iceball') then 1300 1303 if (igcm_h2o_ice.eq.0) then 1301 write(*,*) "No water ice tracer! Can't use this option"1302 stop1304 call abort_physic("newstart", 1305 & "No water ice tracer! Can't use iceball",1) 1303 1306 endif 1304 1307 DO j=1,jjp1 -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/start2archive.F
r3836 r3893 179 179 IF (ierr.NE.NF_NOERR) THEN 180 180 write(6,*)' Pb d''ouverture du fichier'//trim(fichnom) 181 CALL ABORT 181 CALL abort_physic("start2archive", 182 & "Pb d'ouverture du ficher start.nc",1) 182 183 ENDIF 183 184 184 185 ierr = NF_INQ_VARID (nid1, "controle", varid) 185 186 IF (ierr .NE. NF_NOERR) THEN 186 PRINT*, "start2archive: Le champ <controle> est absent"187 CALL abort187 CALL abort_physic("start2archive", 188 & "Le champ <controle> est absent",1) 188 189 ENDIF 189 190 #ifdef NC_DOUBLE … … 193 194 #endif 194 195 IF (ierr .NE. NF_NOERR) THEN 195 PRINT*, "start2archive: Lecture echoue pour <controle>"196 CALL abort196 CALL abort_physic("start2archive", 197 & "Lecture echoue pour <controle>",1) 197 198 ENDIF 198 199 … … 204 205 IF (ierr.NE.NF_NOERR) THEN 205 206 write(6,*)' Pb d''ouverture du fichier'//trim(fichnom) 206 CALL ABORT 207 CALL abort_physic("start2archive", 208 & "Pb d'ouverture du ficher startfi.nc",1) 207 209 ENDIF 208 210 ierr = NF_INQ_DIMID(nid1,"subsurface_layers",varid) 209 211 IF (ierr .NE. NF_NOERR) THEN 210 PRINT*, "start2archive: No subsurface_layers dimension!!"211 CALL abort212 CALL abort_physic("start2archive", 213 & "No subsurface_layers dimension",1) 212 214 ENDIF 213 215 ierr = NF_INQ_DIMLEN(nid1,varid,nsoilmx) 214 216 IF (ierr .NE. NF_NOERR) THEN 215 PRINT*, "start2archive: Failed reading subsurface_layers value!!"216 CALL abort217 CALL abort_physic("start2archive", 218 & "Failed reading subsurface_layers value",1) 217 219 ENDIF 218 220 ierr = NF_CLOSE(nid1) … … 265 267 IF (ierr.NE.NF_NOERR) THEN 266 268 write(6,*)' Pb d''ouverture du fichier'//trim(fichnom) 267 CALL ABORT 269 CALL abort_physic("start2archive", 270 & "Pb d'ouverture du ficher startfi.nc",1) 268 271 ENDIF 269 272 270 273 ierr = NF_INQ_VARID (nid1, "controle", varid) 271 274 IF (ierr .NE. NF_NOERR) THEN 272 PRINT*, "start2archive: Le champ <controle> est absent"273 CALL abort275 CALL abort_physic("start2archive", 276 & "Le champ <controle> est absent",1) 274 277 ENDIF 275 278 #ifdef NC_DOUBLE … … 279 282 #endif 280 283 IF (ierr .NE. NF_NOERR) THEN 281 PRINT*, "start2archive: Lecture echoue pour <controle>"282 CALL abort284 CALL abort_physic("start2archive", 285 & "Lecture echoue pour <controle>",1) 283 286 ENDIF 284 287 … … 290 293 c----------------------------------------------------------------------- 291 294 !mars a voir if ((day_ini_fi.ne.day_ini).or.(abs(timefi-timedyn).gt.1.e-10)) 292 if ((day_ini_fi.ne.day_ini)) 293 & stop ' Probleme de Synchro entre start et startfi !!!' 295 if ((day_ini_fi.ne.day_ini)) then 296 call abort_physic("start2archive", 297 & "Probleme de Synchro entre start et startfi",1) 298 endif 294 299 295 300 … … 447 452 if (ierr.ne.NF_NOERR) then 448 453 write(*,*) "time matter ",NF_STRERROR(ierr) 449 stop 454 call abort_physic("start2archive", 455 & "Wibbly wobbly timey wimey error",1) 450 456 endif 451 457 -
trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/write_archive.F
r3100 r3893 116 116 write(*,*) "***** PUT_VAR problem in write_archive" 117 117 write(*,*) "***** with ",trim(nom)," ",nf_STRERROR(ierr) 118 call abort 118 call abort_physic("write_archive", 119 & "PUT_VAR problem in write_archive", 1) 119 120 endif 120 121 … … 133 134 ierr=NF_INQ_DIMID(nid,"subsurface_layers",id(3)) 134 135 if (ierr.ne.NF_NOERR) then 135 write(*,*)"write_archive: dimension <subsurface_layers>", 136 & " is missing !!!" 137 call abort 136 call abort_physic("write_archive", 137 & "dimension <subsurface_layers> is missing",1) 138 138 endif 139 139 ierr=NF_INQ_DIMID(nid,"Time",id(4)) … … 175 175 ierr=NF_INQ_DIMID(nid,"ocean_layers",id(3)) 176 176 if (ierr.ne.NF_NOERR) then 177 write(*,*)"write_archive: dimension <ocean_layers>", 178 & " is missing !!!" 179 call abort 177 call abort_physic("write_archive", 178 & "dimension <ocean_layers> is missing",1) 180 179 endif 181 180 ierr=NF_INQ_DIMID(nid,"Time",id(4)) … … 246 245 write(*,*) "***** PUT_VAR problem in write_archive" 247 246 write(*,*) "***** with ",nom,nf_STRERROR(ierr) 248 call abort 247 call abort_physic("write_archive", 248 & "PUT_VAR problem in write_archive",1) 249 249 endif 250 250 … … 282 282 write(*,*) "***** PUT_VAR problem in write_archive" 283 283 write(*,*) "***** with ",nom,nf_STRERROR(ierr) 284 call abort 284 call abort_physic("write_archive", 285 & "PUT_VAR problem in write_archive",1) 285 286 endif 286 287 287 288 else 288 289 write(*,*) "write_archive: dim=",dim," ?!?" 289 call abort 290 call abort_physic("write_archive", 291 & "dim problem in write_archive",1) 290 292 endif ! of if (dim.eq.3) else if (dim.eq.-3) .... 291 293
Note: See TracChangeset
for help on using the changeset viewer.
