- Timestamp:
- 12/15/14 17:04:28 (10 years ago)
- Location:
- codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf
- Files:
-
- 1 added
- 4 deleted
- 45 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3dpar/initialize_physics.F90
r253 r313 1 1 SUBROUTINE initialize_unstructured_physics(nbp, nlayer, communicator, nb_proc, distrib, & 2 punjours, pdayref, ptimestep, &2 punjours, pdayref,time0,ptimestep, & 3 3 nb_vertex, lat, lon, area, bounds_lon, bounds_lat, & 4 4 prad,pg,pr,pcpp, preff,ap,bp ) … … 37 37 REAL,INTENT(in) :: bounds_lat(nbp,nb_vertex) ! latitude boundaries of cell 38 38 INTEGER,INTENT(in) :: pdayref ! reference day of for the simulation 39 INTEGER,INTENT(in) :: time0 ! initialtime (s) 39 40 REAL,INTENT(in) :: ptimestep ! physics time step (s) 40 41 REAL,INTENT(in) :: preff ! reference surface pressure (Pa) … … 79 80 rlatd,rlond,airephy,prad,pg,pr,pcpp) 80 81 81 CALL initialize_xios_output 82 CALL initialize_xios_output(time0) 82 83 !$OMP END PARALLEL 83 84 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/aeropacity.F90
r298 r313 208 208 end do 209 209 210 call abort 210 call abort_physiq 211 211 endif 212 212 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/ave_stelspec.F90
r227 r313 88 88 file_id='/stellar_spectra/gj644.txt' 89 89 print*,'Find out tstellar before using this star!' 90 call abort 90 call abort_physiq 91 91 file_id_lam='/stellar_spectra/lam.txt' 92 92 Nfine=5000 … … 113 113 Case Default 114 114 print*,'Error: unknown star type chosen' 115 call abort 115 call abort_physiq 116 116 End Select 117 117 … … 128 128 write(*,*)' datadir = /absolute/path/to/datagcm' 129 129 write(*,*)' Also check that there is a ',trim(file_id_lam),' there.' 130 call abort 130 call abort_physiq 131 131 else 132 132 do ifine=1,Nfine … … 147 147 write(*,*)' datadir = /absolute/path/to/datagcm' 148 148 write(*,*)' Also check that there is a ',trim(file_id),' there.' 149 call abort 149 call abort_physiq 150 150 else 151 151 do ifine=1,Nfine -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/calc_cpp_mugaz.F90
r298 r313 79 79 else 80 80 print*,'Error in calc_cpp_mugaz: Gas species not recognised!' 81 call abort 81 call abort_physiq 82 82 endif 83 83 endif … … 102 102 print*,' Either adjust cpp / mugaz via newstart to calculated values,' 103 103 print*,' or set check_cpp_match to .false. in callphys.def.' 104 stop104 CALL abort_physiq 105 105 else 106 106 if (is_master) print*,'--> OK. Settings match composition.' -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callcorrk.F90
r298 r313 211 211 if(naerkind.gt.4)then 212 212 print*,'Code not general enough to deal with naerkind > 4 yet.' 213 call abort 213 call abort_physiq 214 214 endif 215 215 call su_aer_radii(ngrid,nlayer,reffrad,nueffrad) … … 235 235 if((igcm_h2o_vap.eq.0) .and. varactive)then 236 236 print*,'varactive in callcorrk but no h2o_vap tracer.' 237 stop237 CALL abort_physiq 238 238 endif 239 239 … … 250 250 if (global1d.and.diurnal) then 251 251 print*,'if global1d is true, diurnal must be set to false' 252 stop252 CALL abort_physiq 253 253 endif 254 254 … … 407 407 print*,'Serious problems with qsvaer values' 408 408 print*,'in callcorrk' 409 call abort 409 call abort_physiq 410 410 endif 411 411 if(qsvaer(k,nw,iaer).gt.qxvaer(k,nw,iaer))then … … 418 418 print*,'Serious problems with qsiaer values' 419 419 print*,'in callcorrk' 420 call abort 420 call abort_physiq 421 421 endif 422 422 if(qsiaer(k,nw,iaer).gt.qxiaer(k,nw,iaer))then … … 457 457 print*,'For open lower boundary in callcorrk must' 458 458 print*,'have surface albedo set to zero!' 459 call abort 459 call abort_physiq 460 460 endif 461 461 … … 556 556 if(nq.gt.1)then 557 557 print*,'Need 1 tracer only to run kcm1d.e' 558 stop558 CALL abort_physiq 559 559 endif 560 560 do l=1,nlayer … … 639 639 print*,'Minimum pressure is outside the radiative' 640 640 print*,'transfer kmatrix bounds, exiting.' 641 call abort 641 call abort_physiq 642 642 elseif(plevrad(L_LEVELS).gt.pgasmax)then 643 643 print*,'Maximum pressure is outside the radiative' 644 644 print*,'transfer kmatrix bounds, exiting.' 645 call abort 645 call abort_physiq 646 646 endif 647 647 … … 654 654 print*,"tgasmin=",tgasmin 655 655 if (strictboundcorrk) then 656 call abort 656 call abort_physiq 657 657 else 658 658 print*,'***********************************************' … … 669 669 print*,"tgasmax=",tgasmax 670 670 if (strictboundcorrk) then 671 call abort 671 call abort_physiq 672 672 else 673 673 print*,'***********************************************' … … 687 687 print*,"tgasmin=",tgasmin 688 688 if (strictboundcorrk) then 689 call abort 689 call abort_physiq 690 690 else 691 691 print*,'***********************************************' … … 702 702 print*,"tgasmax=",tgasmax 703 703 if (strictboundcorrk) then 704 call abort 704 call abort_physiq 705 705 else 706 706 print*,'***********************************************' … … 787 787 print*,'temp= ',pt(ig,:) 788 788 print*,'pplay= ',pplay(ig,:) 789 call abort 789 call abort_physiq 790 790 endif 791 791 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callsedim.F
r227 r313 81 81 write(*,*) "callsedim error: water=.true.", 82 82 & " but igcm_h2o_ice=0" 83 stop83 CALL abort_physiq 84 84 endif 85 85 if (iaero_h2o.eq.0) then 86 86 write(*,*) "callsedim error: water=.true.", 87 87 & " but iaero_ho2=0" 88 stop88 CALL abort_physiq 89 89 endif 90 90 endif -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/condense_cloud.F90
r227 r313 175 175 print*,'In condens_cloud but no CO2 ice tracer, exiting.' 176 176 print*,'Still need generalisation to arbitrary species!' 177 stop177 CALL abort_physiq 178 178 endif 179 179 … … 401 401 write(116,*) 0.0, pplev(1,1), 0.0, 0.0 402 402 close(116) 403 call abort 403 call abort_physiq 404 404 endif 405 405 … … 448 448 PRINT*,'Ps = ',pplev(ig,1) 449 449 PRINT*,'d Ps = ',pdpsrf(ig) 450 STOP450 CALL abort_physiq 451 451 ENDIF 452 452 END IF -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/convadj.F
r227 r313 97 97 if (noms(iq).eq."co2") then 98 98 print*,'dont go there' 99 stop99 CALL abort_physiq 100 100 ico2=iq 101 101 m_co2 = 44.01E-3 ! CO2 molecular mass (kg/mol) … … 297 297 zalpha=1. 298 298 ELSE 299 ! IF(zalpha.LT.0.) STOP299 ! IF(zalpha.LT.0.) CALL abort_physiq 300 300 IF(zalpha.LT.1.e-4) zalpha=1.e-4 301 301 ENDIF … … 381 381 print*,'jadrs=',jadrs 382 382 383 call abort 383 call abort_physiq 384 384 endif 385 385 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/cp_neutral.F90
r222 r313 22 22 else 23 23 print*,'Gas not recognised in cp_neutral!' 24 call abort 24 call abort_physiq 25 25 endif 26 26 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/datareadnc.F
r222 r313 118 118 write(*,*)' can be obtained online on:' 119 119 write(*,*)' http://www.lmd.jussieu.fr/~forget/datagcm/datafile' 120 STOP120 CALL abort_physiq 121 121 ENDIF 122 122 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/def_var.F90
r222 r313 39 39 write(*,*) "def_var: Failed defining variable "//trim(name) 40 40 write(*,*) NF_STRERROR(ierr) 41 stop ""41 CALL abort_physiq 42 42 endif 43 43 … … 48 48 write(*,*) "def_var: Failed writing title attribute for "//trim(name) 49 49 write(*,*) NF_STRERROR(ierr) 50 stop ""50 CALL abort_physiq 51 51 endif 52 52 … … 56 56 write(*,*) "def_var: Failed writing units attribute for "//trim(name) 57 57 write(*,*) NF_STRERROR(ierr) 58 stop ""58 CALL abort_physiq 59 59 endif 60 60 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/gfluxi.F
r265 r313 62 62 63 63 64 IF (NLL .GT. NLP) STOP 'PARAMETER NL TOO SMALL IN GFLUXI' 65 64 IF (NLL .GT. NLP) THEN 65 PRINT*, 'PARAMETER NL TOO SMALL IN GFLUXI' 66 CALL abort_physiq 67 ENDIF 68 66 69 NLAYER = L_NLAYRAD 67 70 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/gradients_kcm.F90
r222 r313 43 43 if(ngasmx.eq.1)then 44 44 print*,'Cannot have moist adiabat with one gas...' 45 stop45 CALL abort_physiq 46 46 endif 47 47 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/hydrol.F90
r227 r313 151 151 print*,'How are we supposed to average the ocean' 152 152 print*,'temperature, when there are no oceans?' 153 call abort 153 call abort_physiq 154 154 endif 155 155 … … 157 157 print*,'You have enabled runoff, but you have no oceans.' 158 158 print*,'Where did you think the water was going to go?' 159 call abort 159 call abort_physiq 160 160 endif 161 161 … … 315 315 print*,'Surface type not recognised in hydrol.F!' 316 316 print*,'Exiting...' 317 call abort 317 call abort_physiq 318 318 319 319 endif -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iniaerosol.F
r298 r313 71 71 print*, 'or change options in callphys.def' 72 72 print*, 'Abort in iniaerosol.F' 73 call abort 73 call abort_physiq 74 74 endif 75 75 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inifis.F
r310 r313 215 215 if (is_master) print*,'If diurnal=true, we should turn off', 216 216 & ' tlocked.' 217 stop217 CALL abort_physiq 218 218 endif 219 219 … … 306 306 if (is_master) print*,'We need a CO2 ice tracer to', 307 307 & ' condense CO2' 308 call abort 308 call abort_physiq 309 309 endif 310 310 … … 393 393 if (kastprof.and.(ngrid.gt.1)) then 394 394 if (is_master) print*,'kastprof can only be used in 1D!' 395 call abort 395 call abort_physiq 396 396 endif 397 397 … … 411 411 if (is_master) print*,'nosurf not compatible with soil', 412 412 & ' scheme! ... got to make a choice!' 413 call abort 413 call abort_physiq 414 414 endif 415 415 … … 430 430 if (is_master) print*,'newtonian not compatible with ', 431 431 & 'correlated-k!' 432 call abort 432 call abort_physiq 433 433 endif 434 434 if (newtonian.and.calladj) then 435 435 if (is_master) print*,'newtonian not compatible with ', 436 436 & 'adjustment!' 437 call abort 437 call abort_physiq 438 438 endif 439 439 if (newtonian.and.calldifv) then 440 440 if (is_master) print*,'newtonian not compatible with a ', 441 441 & 'boundary layer!' 442 call abort 442 call abort_physiq 443 443 endif 444 444 … … 452 452 if (testradtimes.and.(ngrid.gt.1)) then 453 453 if (is_master) print*,'testradtimes can only be used in 1D!' 454 call abort 454 call abort_physiq 455 455 endif 456 456 … … 624 624 if (varactive.and.varfixed) then 625 625 if (is_master) print*,'if varactive, varfixed must be OFF!' 626 stop626 CALL abort_physiq 627 627 endif 628 628 … … 642 642 if (is_master) print*,'if water is ON, tracer must be ', 643 643 & 'ON too!' 644 stop644 CALL abort_physiq 645 645 endif 646 646 … … 655 655 if (is_master) print*,'if watercond is used, water should ', 656 656 & 'be used too' 657 stop657 CALL abort_physiq 658 658 endif 659 659 … … 706 706 if (is_master) PRINT *, "mugaz must be set if ', 707 707 & 'force_cpp = T" 708 STOP708 CALL abort_physiq 709 709 ELSE 710 710 if (is_master) write(*,*) "inifis: mugaz=",mugaz … … 722 722 if (is_master) PRINT *, "cpp must be set if ', 723 723 & 'force_cpp = T" 724 STOP724 CALL abort_physiq 725 725 ELSE 726 726 if (is_master) write(*,*) "inifis: cpp=",cpp … … 739 739 if (is_master) write(*,*) 'Cannot read file callphys.def.", 740 740 & " Is it here ?' 741 stop741 CALL abort_physiq 742 742 ENDIF 743 743 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inistats.F
r222 r313 27 27 nsteppd=nint(daysec/dtphys) 28 28 write (*,*) 'nsteppd=',nsteppd 29 if (abs(float(nsteppd)-daysec/dtphys).gt.1.e-8*daysec) 30 & stop'Dans Instat: 1jour .ne. n pas physiques' 29 if (abs(float(nsteppd)-daysec/dtphys).gt.1.e-8*daysec) THEN 30 PRINT *,'Dans Instat: 1jour .ne. n pas physiques' 31 CALL abort_physiq 32 endif 31 33 32 if(mod(nsteppd,istime).ne.0) 33 & stop'Dans Instat: 1jour .ne. n*istime pas physiques' 34 if(mod(nsteppd,istime).ne.0) THEN 35 PRINT*,'Dans Instat: 1jour .ne. n*istime pas physiques' 36 CALL abort_physiq 37 endif 34 38 35 39 istats=nsteppd/istime … … 50 54 if (ierr.ne.NF_NOERR) then 51 55 write (*,*) NF_STRERROR(ierr) 52 stop ""56 CALL abort_physiq 53 57 endif 54 58 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/interpolateH2H2.F90
r298 r313 61 61 print*,'really want to run simulations with hydrogen at T > 400 K, contact' 62 62 print*,'Robin Wordsworth [rwordsworth@uchicago.edu].' 63 stop63 CALL abort_physiq 64 64 endif 65 65 … … 82 82 write(*,*) 'datadir = /absolute/path/to/datagcm' 83 83 write(*,*) 'Also check that the continuum data continuum_data/H2-H2_norm_2011.cia is there.' 84 call abort 84 call abort_physiq 85 85 else 86 86 … … 92 92 print*,'is ',nres,', which does not match nS.' 93 93 print*,'Please adjust nS value in interpolateH2H2.F90' 94 stop94 CALL abort_physiq 95 95 endif 96 96 temp_arr(iT)=Ttemp -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iostart.F90
r270 r313 57 57 write(*,*)'open_startphy: problem opening file '//trim(filename) 58 58 write(*,*)trim(nf90_strerror(ierr)) 59 !CALL ABORT59 !CALL abort_physiq 60 60 found=.false. 61 61 ELSE … … 125 125 //trim(field_name) 126 126 write(*,*)trim(nf90_strerror(ierr)) 127 CALL ABORT127 CALL abort_physiq 128 128 ENDIF 129 129 ENDIF … … 176 176 //trim(field_name) 177 177 write(*,*)trim(nf90_strerror(ierr)) 178 CALL ABORT178 CALL abort_physiq 179 179 ENDIF 180 180 ENDIF … … 315 315 IF (.NOT. tmp_found) THEN 316 316 PRINT*, 'get_field_rgen: Field <'//field_name//'> not found' 317 CALL abort 317 CALL abort_physiq 318 318 ENDIF 319 319 ENDIF … … 335 335 ! IF (ierr/=NF90_NOERR) THEN 336 336 ! PRINT*, 'phyetat0: Lecture echouee aussi en 2D pour <'//field_name//'>' 337 ! CALL abort 337 ! CALL abort_physiq 338 338 ! ELSE 339 339 ! PRINT*, 'phyetat0: La variable <'//field_name//'> lu sur surface seulement'!, selon ancien format, le reste mis a zero' 340 340 ! END IF 341 341 ! ELSE 342 CALL abort 342 CALL abort_physiq 343 343 ! ENDIF 344 344 ENDIF … … 435 435 IF (ierr/=NF90_NOERR) THEN 436 436 PRINT*, 'phyetat0: Failed loading <'//trim(var_name)//'>' 437 CALL abort 437 CALL abort_physiq 438 438 ENDIF 439 439 tmp_found=.TRUE. … … 455 455 IF (.NOT. tmp_found) THEN 456 456 PRINT*, 'phyetat0: Variable <'//trim(var_name)//'> not found' 457 CALL abort 457 CALL abort_physiq 458 458 ENDIF 459 459 ENDIF … … 490 490 write(*,*)'open_restartphy: problem creating file '//trim(filename) 491 491 write(*,*)trim(nf90_strerror(ierr)) 492 CALL ABORT492 CALL abort_physiq 493 493 ENDIF 494 494 already_created=.true. … … 499 499 write(*,*)'open_restartphy: problem opening file '//trim(filename) 500 500 write(*,*)trim(nf90_strerror(ierr)) 501 CALL ABORT501 CALL abort_physiq 502 502 ENDIF 503 503 return … … 515 515 write(*,*)'open_restartphy: problem defining index dimension ' 516 516 write(*,*)trim(nf90_strerror(ierr)) 517 CALL ABORT517 CALL abort_physiq 518 518 ENDIF 519 519 … … 522 522 write(*,*)'open_restartphy: problem defining physical_points dimension ' 523 523 write(*,*)trim(nf90_strerror(ierr)) 524 CALL ABORT524 CALL abort_physiq 525 525 ENDIF 526 526 … … 529 529 write(*,*)'open_restartphy: problem defining subsurface_layers dimension ' 530 530 write(*,*)trim(nf90_strerror(ierr)) 531 CALL ABORT531 CALL abort_physiq 532 532 ENDIF 533 533 … … 536 536 write(*,*)'open_restartphy: problem defining nlayer_plus_1 dimension ' 537 537 write(*,*)trim(nf90_strerror(ierr)) 538 CALL ABORT538 CALL abort_physiq 539 539 ENDIF 540 540 … … 545 545 write(*,*)'open_restartphy: problem defining number_of_advected_fields dimension ' 546 546 write(*,*)trim(nf90_strerror(ierr)) 547 CALL ABORT547 CALL abort_physiq 548 548 ENDIF 549 549 endif … … 553 553 write(*,*)'open_restartphy: problem defining nlayer dimension ' 554 554 write(*,*)trim(nf90_strerror(ierr)) 555 CALL ABORT555 CALL abort_physiq 556 556 ENDIF 557 557 … … 560 560 write(*,*)'open_restartphy: problem defining Time dimension ' 561 561 write(*,*)trim(nf90_strerror(ierr)) 562 CALL ABORT562 CALL abort_physiq 563 563 ENDIF 564 564 … … 567 567 write(*,*)'open_restartphy: problem defining oceanic layer dimension ' 568 568 write(*,*)trim(nf90_strerror(ierr)) 569 CALL ABORT569 CALL abort_physiq 570 570 ENDIF 571 571 … … 575 575 write(*,*)'open_restartphy: problem ending definition mode ' 576 576 write(*,*)trim(nf90_strerror(ierr)) 577 CALL ABORT577 CALL abort_physiq 578 578 ENDIF 579 579 ENDIF … … 884 884 PRINT *, "Error phyredem(put_field_rgen) : wrong dimension for ",trim(field_name) 885 885 write(*,*) " field_size =",field_size 886 CALL ABORT886 CALL abort_physiq 887 887 ENDIF 888 888 … … 891 891 write(*,*) " Error phyredem(put_field_rgen) : failed writing ",trim(field_name) 892 892 write(*,*)trim(nf90_strerror(ierr)) 893 call abort 893 call abort_physiq 894 894 endif 895 895 … … 996 996 write(*,*)'put_var_rgen: problem writing Time' 997 997 write(*,*)trim(nf90_strerror(ierr)) 998 CALL ABORT998 CALL abort_physiq 999 999 ENDIF 1000 1000 return ! nothing left to do … … 1011 1011 PRINT *, "put_var_rgen error : wrong dimension" 1012 1012 write(*,*) " var_size =",var_size 1013 CALL abort 1013 CALL abort_physiq 1014 1014 1015 1015 ENDIF ! of IF (var_size==length) THEN … … 1032 1032 write(*,*)'put_var_rgen: problem writing '//trim(var_name) 1033 1033 write(*,*)trim(nf90_strerror(ierr)) 1034 CALL ABORT1034 CALL abort_physiq 1035 1035 ENDIF 1036 1036 ENDIF ! of IF (is_master) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/kcm1d.F90
r227 r313 117 117 write(*,*) ' )' 118 118 write(*,*) ' ... might as well stop here ...' 119 stop119 CALL abort_physiq 120 120 else 121 121 close(90) … … 133 133 write(*,*) 'Please remove the file and restart the run.' 134 134 write(*,*) 'Runtime parameters are supposed to be in kcm1d.def' 135 stop135 CALL abort_physiq 136 136 else 137 137 call system('touch run.def') … … 144 144 if(.not.global1d)then 145 145 print*,'Error, kcm1d must have global1d=.true. in kcm1d.def!' 146 stop146 CALL abort_physiq 147 147 end if 148 148 … … 200 200 print*,"In 1D modeling, check_cpp_match is supposed to be F" 201 201 print*,"Please correct callphys.def" 202 stop202 CALL abort_physiq 203 203 endif 204 204 … … 221 221 write(*,*) "kcm1d: error reading number of tracers" 222 222 write(*,*) " (first line of traceur.def) " 223 stop223 CALL abort_physiq 224 224 endif 225 225 nqtot=nq … … 235 235 if (ierr.ne.0) then 236 236 write(*,*) 'kcm1d: error reading tracer names...' 237 stop237 CALL abort_physiq 238 238 endif 239 239 enddo !of do iq=1,nq -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/kcmprof_fn.F90
r227 r313 95 95 print*,'Must have Psat_max=0 if no variable species' 96 96 psat_max=0.0 97 ! stop97 !CALL abort_physiq 98 98 endif 99 99 print*, 'Assuming pure atmosphere' … … 109 109 m_v = 16.04/1000. 110 110 tcrit = 1.91d2 111 stop111 CALL abort_physiq 112 112 else 113 113 print*,'Variable gas not recognised!' 114 call abort 114 call abort_physiq 115 115 endif 116 116 … … 326 326 if(ilay.eq.1)then 327 327 print*,'Error in create_profils: Psurf here less than Psurf in RCM!' 328 call abort 328 call abort_physiq 329 329 endif 330 330 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/lect_start_archive.F
r227 r313 21 21 c 22 22 c Objet: Lecture des variables d'un fichier "start_archive" 23 c Plus besoin de r égler ancienne valeurs grace23 c Plus besoin de rï¿œgler ancienne valeurs grace 24 24 c a l'allocation dynamique de memoire (Yann Wanherdrick) 25 25 c … … 341 341 IF (ierr .NE. NF_NOERR) THEN 342 342 PRINT*, "Lect_start_archive: champ <controle> est absent" 343 CALL abort 343 CALL abort_physiq 344 344 ENDIF 345 345 #ifdef NC_DOUBLE … … 350 350 IF (ierr .NE. NF_NOERR) THEN 351 351 PRINT*, "lect_start_archive: Lecture echoue pour <controle>" 352 CALL abort 352 CALL abort_physiq 353 353 ENDIF 354 354 c … … 362 362 IF (ierr .NE. NF_NOERR) THEN 363 363 PRINT*, "lect_start_archive: Le champ <rlonv> est absent" 364 CALL abort 364 CALL abort_physiq 365 365 ENDIF 366 366 #ifdef NC_DOUBLE … … 371 371 IF (ierr .NE. NF_NOERR) THEN 372 372 PRINT*, "lect_start_archive: Lecture echouee pour <rlonv>" 373 CALL abort 373 CALL abort_physiq 374 374 ENDIF 375 375 c … … 377 377 IF (ierr .NE. NF_NOERR) THEN 378 378 PRINT*, "lect_start_archive: Le champ <rlatu> est absent" 379 CALL abort 379 CALL abort_physiq 380 380 ENDIF 381 381 #ifdef NC_DOUBLE … … 386 386 IF (ierr .NE. NF_NOERR) THEN 387 387 PRINT*, "lect_start_archive: Lecture echouee pour <rlatu>" 388 CALL abort 388 CALL abort_physiq 389 389 ENDIF 390 390 c … … 392 392 IF (ierr .NE. NF_NOERR) THEN 393 393 PRINT*, "lect_start_archive: Le champ <rlonu> est absent" 394 CALL abort 394 CALL abort_physiq 395 395 ENDIF 396 396 #ifdef NC_DOUBLE … … 401 401 IF (ierr .NE. NF_NOERR) THEN 402 402 PRINT*, "lect_start_archive: Lecture echouee pour <rlonu>" 403 CALL abort 403 CALL abort_physiq 404 404 ENDIF 405 405 c … … 407 407 IF (ierr .NE. NF_NOERR) THEN 408 408 PRINT*, "lect_start_archive: Le champ <rlatv> est absent" 409 CALL abort 409 CALL abort_physiq 410 410 ENDIF 411 411 #ifdef NC_DOUBLE … … 416 416 IF (ierr .NE. NF_NOERR) THEN 417 417 PRINT*, "lect_start_archive: Lecture echouee pour <rlatv>" 418 CALL abort 418 CALL abort_physiq 419 419 ENDIF 420 420 c … … 447 447 IF (ierr .NE. NF_NOERR) THEN 448 448 PRINT*, "Nothing to do..." 449 CALL abort 449 CALL abort_physiq 450 450 ENDIF 451 451 ENDIF … … 457 457 IF (ierr .NE. NF_NOERR) THEN 458 458 PRINT*, "lect_start_archive: Lecture echouee pour <bps>" 459 CALL abort 459 CALL abort_physiq 460 460 END IF 461 461 … … 482 482 if (ierr .NE. NF_NOERR) then 483 483 PRINT*, "lect_start_archive: Failed reading <soildepth>" 484 CALL abort 484 CALL abort_physiq 485 485 endif 486 486 … … 510 510 if (ierr .NE. NF_NOERR) then 511 511 PRINT*, "lect_start_archive: Failed reading <soildepth>" 512 CALL abort 512 CALL abort_physiq 513 513 endif 514 514 endif ! of if (depthinterpol) … … 535 535 if (ierr .NE. NF_NOERR) then 536 536 PRINT*, "lect_start_archive: Failed reading <inertiedat>" 537 CALL abort 537 CALL abort_physiq 538 538 endif 539 539 endif … … 546 546 IF (ierr .NE. NF_NOERR) THEN 547 547 PRINT*, "lect_start_archive: Le champ <phisinit> est absent" 548 CALL abort 548 CALL abort_physiq 549 549 ENDIF 550 550 #ifdef NC_DOUBLE … … 555 555 IF (ierr .NE. NF_NOERR) THEN 556 556 PRINT*, "lect_start_archive: Lecture echouee pour <phisinit>" 557 CALL abort 557 CALL abort_physiq 558 558 ENDIF 559 559 … … 575 575 IF (ierr .NE. NF_NOERR) THEN 576 576 PRINT*, "lect_start_archive: Le champ <Time> est absent" 577 CALL abort 577 CALL abort_physiq 578 578 endif 579 579 ENDIF … … 590 590 IF (ierr .NE. NF_NOERR) THEN 591 591 PRINT*, "lect_start_archive: Lecture echouee pour <Time>" 592 CALL abort 592 CALL abort_physiq 593 593 ENDIF 594 594 c … … 649 649 ! IF (ierr .NE. NF_NOERR) THEN 650 650 ! PRINT*, "lect_start_archive: Le champ <co2ice> est absent" 651 ! CALL abort 651 ! CALL abort_physiq 652 652 ! ENDIF 653 653 !#ifdef NC_DOUBLE … … 659 659 ! PRINT*, "lect_start_archive: Lecture echouee pour <co2ice>" 660 660 ! PRINT*, NF_STRERROR(ierr) 661 ! CALL abort 661 ! CALL abort_physiq 662 662 ! ENDIF 663 663 c … … 665 665 IF (ierr .NE. NF_NOERR) THEN 666 666 PRINT*, "lect_start_archive: Le champ <emis> est absent" 667 CALL abort 667 CALL abort_physiq 668 668 ENDIF 669 669 #ifdef NC_DOUBLE … … 674 674 IF (ierr .NE. NF_NOERR) THEN 675 675 PRINT*, "lect_start_archive: Lecture echouee pour <emis>" 676 CALL abort 676 CALL abort_physiq 677 677 ENDIF 678 678 c … … 680 680 IF (ierr .NE. NF_NOERR) THEN 681 681 PRINT*, "lect_start_archive: Le champ <ps> est absent" 682 CALL abort 682 CALL abort_physiq 683 683 ENDIF 684 684 #ifdef NC_DOUBLE … … 689 689 IF (ierr .NE. NF_NOERR) THEN 690 690 PRINT*, "lect_start_archive: Lecture echouee pour <ps>" 691 CALL abort 691 CALL abort_physiq 692 692 ENDIF 693 693 c … … 695 695 IF (ierr .NE. NF_NOERR) THEN 696 696 PRINT*, "lect_start_archive: Le champ <tsurf> est absent" 697 CALL abort 697 CALL abort_physiq 698 698 ENDIF 699 699 #ifdef NC_DOUBLE … … 704 704 IF (ierr .NE. NF_NOERR) THEN 705 705 PRINT*, "lect_start_archive: Lecture echouee pour <tsurf>" 706 CALL abort 706 CALL abort_physiq 707 707 ENDIF 708 708 c … … 710 710 IF (ierr .NE. NF_NOERR) THEN 711 711 PRINT*, "lect_start_archive: Le champ <q2surf> est absent" 712 CALL abort 712 CALL abort_physiq 713 713 ENDIF 714 714 #ifdef NC_DOUBLE … … 719 719 IF (ierr .NE. NF_NOERR) THEN 720 720 PRINT*, "lect_start_archive: Lecture echouee pour <q2surf>" 721 CALL abort 721 CALL abort_physiq 722 722 ENDIF 723 723 c … … 833 833 834 834 ! print*,'RDW has added hack to let me continue...' 835 ! CALL abort 835 ! CALL abort_physiq 836 836 ENDIF 837 837 #ifdef NC_DOUBLE … … 871 871 PRINT*, "lect_start_archive: ", 872 872 & "Field <","Tg"//str2,"> not found" 873 CALL abort 873 CALL abort_physiq 874 874 ENDIF 875 875 #ifdef NC_DOUBLE … … 883 883 PRINT*, "lect_start_archive: ", 884 884 & "Failed reading <","Tg"//str2,">" 885 CALL abort 885 CALL abort_physiq 886 886 ENDIF 887 887 c … … 897 897 if (ierr.ne.NF_NOERR) then 898 898 write(*,*)"lect_start_archive: Cannot find <tsoil>" 899 call abort 899 call abort_physiq 900 900 else 901 901 #ifdef NC_DOUBLE … … 915 915 ! if (ierr.ne.NF_NOERR) then 916 916 ! write(*,*)"lect_start_archive: Cannot find <inertiedat>" 917 ! call abort 917 ! call abort_physiq 918 918 ! else 919 919 !#ifdef NC_DOUBLE … … 936 936 IF (ierr .NE. NF_NOERR) THEN 937 937 PRINT*, "lect_start_archive: Le champ <temp> est absent" 938 CALL abort 938 CALL abort_physiq 939 939 ENDIF 940 940 #ifdef NC_DOUBLE … … 945 945 IF (ierr .NE. NF_NOERR) THEN 946 946 PRINT*, "lect_start_archive: Lecture echouee pour <temp>" 947 CALL abort 947 CALL abort_physiq 948 948 ENDIF 949 949 c … … 951 951 IF (ierr .NE. NF_NOERR) THEN 952 952 PRINT*, "lect_start_archive: Le champ <u> est absent" 953 CALL abort 953 CALL abort_physiq 954 954 ENDIF 955 955 #ifdef NC_DOUBLE … … 960 960 IF (ierr .NE. NF_NOERR) THEN 961 961 PRINT*, "lect_start_archive: Lecture echouee pour <u>" 962 CALL abort 962 CALL abort_physiq 963 963 ENDIF 964 964 c … … 966 966 IF (ierr .NE. NF_NOERR) THEN 967 967 PRINT*, "lect_start_archive: Le champ <v> est absent" 968 CALL abort 968 CALL abort_physiq 969 969 ENDIF 970 970 #ifdef NC_DOUBLE … … 975 975 IF (ierr .NE. NF_NOERR) THEN 976 976 PRINT*, "lect_start_archive: Lecture echouee pour <v>" 977 CALL abort 977 CALL abort_physiq 978 978 ENDIF 979 979 c … … 981 981 IF (ierr .NE. NF_NOERR) THEN 982 982 PRINT*, "lect_start_archive: Le champ <q2atm> est absent" 983 CALL abort 983 CALL abort_physiq 984 984 ENDIF 985 985 #ifdef NC_DOUBLE … … 990 990 IF (ierr .NE. NF_NOERR) THEN 991 991 PRINT*, "lect_start_archive: Lecture echouee pour <q2atm>" 992 CALL abort 992 CALL abort_physiq 993 993 ENDIF 994 994 c … … 1006 1006 PRINT*, "lect_start_archive: ", 1007 1007 & " Tracer <",trim(txt),"> not found" 1008 ! CALL abort 1008 ! CALL abort_physiq 1009 1009 ENDIF 1010 1010 #ifdef NC_DOUBLE … … 1226 1226 !print*,'Problem in lect_start_archive interpolating' 1227 1227 !print*,'to new resolution!!' 1228 !call abort 1228 !call abort_physiq 1229 1229 !endif 1230 1230 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/mass_redistribution.F90
r227 r313 176 176 PRINT*,'Ps = ',pplev(ig,1) 177 177 PRINT*,'d Ps = ',pdpsrfmr(ig)*ptimestep 178 STOP178 CALL abort_physiq 179 179 ENDIF 180 180 enddo ! of DO ig=1,ngrid … … 379 379 print*,q 380 380 print*,qm 381 stop381 CALL abort_physiq 382 382 end if 383 383 else ! if(w(l+1).lt.0) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/mod_phys_lmdz_mpi_transfert.F90
r245 r313 60 60 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 61 61 62 !! -- Les chaine de charact ère -- !!62 !! -- Les chaine de charactï¿œre -- !! 63 63 64 64 SUBROUTINE bcast_mpi_c(var1) … … 1321 1321 if (ierr.ne.MPI_SUCCESS) then 1322 1322 write(*,*) "bcast_mpi error: ierr=",ierr 1323 stop1323 CALL abort_physiq 1324 1324 endif 1325 1325 #endif -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/newstart.F
r227 r313 231 231 write(6,*)' Problem opening file:',fichnom 232 232 write(6,*)' ierr = ', ierr 233 CALL ABORT233 CALL abort_physiq 234 234 ENDIF 235 235 tab0 = 50 … … 247 247 write(6,*)' Problem opening file:',fichnom 248 248 write(6,*)' ierr = ', ierr 249 CALL ABORT249 CALL abort_physiq 250 250 ENDIF 251 251 … … 255 255 write(6,*)' Problem opening file:',fichnom 256 256 write(6,*)' ierr = ', ierr 257 CALL ABORT257 CALL abort_physiq 258 258 ENDIF 259 259 … … 1024 1024 if (igcm_h2o_vap.eq.0) then 1025 1025 write(*,*) "No water vapour tracer! Can't use this option" 1026 stop1026 CALL abort_physiq 1027 1027 endif 1028 1028 DO l=1,llm … … 1044 1044 if (igcm_h2o_ice.eq.0) then 1045 1045 write(*,*) "No water ice tracer! Can't use this option" 1046 stop1046 CALL abort_physiq 1047 1047 endif 1048 1048 do ig=1,ngridmx … … 1061 1061 if (igcm_h2o_ice.eq.0) then 1062 1062 write(*,*) "No water ice tracer! Can't use this option" 1063 stop1063 CALL abort_physiq 1064 1064 endif 1065 1065 … … 1103 1103 if (igcm_h2o_ice.eq.0) then 1104 1104 write(*,*) "No water ice tracer! Can't use this option" 1105 stop1105 CALL abort_physiq 1106 1106 endif 1107 1107 … … 1144 1144 if (igcm_h2o_ice.eq.0) then 1145 1145 write(*,*) "No water ice tracer! Can't use this option" 1146 stop1146 CALL abort_physiq 1147 1147 endif 1148 1148 DO j=1,jjp1 … … 1173 1173 if (igcm_h2o_ice.eq.0) then 1174 1174 write(*,*) "No water ice tracer! Can't use this option" 1175 stop1175 CALL abort_physiq 1176 1176 endif 1177 1177 DO j=1,jjp1 … … 1205 1205 if (igcm_h2o_ice.eq.0) then 1206 1206 write(*,*) "No water ice tracer! Can't use this option" 1207 stop1207 CALL abort_physiq 1208 1208 endif 1209 1209 DO j=1,jjp1 … … 1229 1229 if (igcm_h2o_ice.eq.0) then 1230 1230 write(*,*) "No water ice tracer! Can't use this option" 1231 stop1231 CALL abort_physiq 1232 1232 endif 1233 1233 DO j=1,jjp1 … … 1258 1258 if (igcm_h2o_ice.eq.0) then 1259 1259 write(*,*) "No water ice tracer! Can't use this option" 1260 stop1260 CALL abort_physiq 1261 1261 endif 1262 1262 DO j=1,jjp1 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/phyetat0.F90
r227 r313 16 16 !====================================================================== 17 17 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818 18 ! Adaptation àMars : Yann Wanherdrick18 ! Adaptation ï¿œ Mars : Yann Wanherdrick 19 19 ! Objet: Lecture de l etat initial pour la physique 20 20 !====================================================================== … … 113 113 ! IF (ierr.NE.NF_NOERR) THEN 114 114 ! PRINT*, 'phyetat0: Le champ <latitude> est absent' 115 ! CALL abort 115 ! CALL abort_physiq 116 116 ! ENDIF 117 117 !#ifdef NC_DOUBLE … … 122 122 ! IF (ierr.NE.NF_NOERR) THEN 123 123 ! PRINT*, 'phyetat0: Lecture echouee pour <latitude>' 124 ! CALL abort 124 ! CALL abort_physiq 125 125 ! ENDIF 126 126 !c … … 130 130 ! IF (ierr.NE.NF_NOERR) THEN 131 131 ! PRINT*, 'phyetat0: Le champ <longitude> est absent' 132 ! CALL abort 132 ! CALL abort_physiq 133 133 ! ENDIF 134 134 !#ifdef NC_DOUBLE … … 139 139 ! IF (ierr.NE.NF_NOERR) THEN 140 140 ! PRINT*, 'phyetat0: Lecture echouee pour <longitude>' 141 ! CALL abort 141 ! CALL abort_physiq 142 142 ! ENDIF 143 143 !c … … 147 147 ! IF (ierr.NE.NF_NOERR) THEN 148 148 ! PRINT*, 'phyetat0: Le champ <area> est absent' 149 ! CALL abort 149 ! CALL abort_physiq 150 150 ! ENDIF 151 151 !#ifdef NC_DOUBLE … … 156 156 ! IF (ierr.NE.NF_NOERR) THEN 157 157 ! PRINT*, 'phyetat0: Lecture echouee pour <area>' 158 ! CALL abort 158 ! CALL abort_physiq 159 159 ! ENDIF 160 160 ! xmin = 1.0E+20 … … 168 168 if (.not.found) then 169 169 write(*,*) "phyetat0: Failed loading <phisfi>" 170 call abort 170 call abort_physiq 171 171 else 172 172 write(*,*) "phyetat0: surface geopotential <phisfi> range:", & … … 178 178 if (.not.found) then 179 179 write(*,*) "phyetat0: Failed loading <albedodat>" 180 call abort 180 call abort_physiq 181 181 else 182 182 write(*,*) "phyetat0: Bare ground albedo <albedodat> range:", & … … 188 188 if (.not.found) then 189 189 write(*,*) "phyetat0: Failed loading <ZMEA>" 190 call abort 190 call abort_physiq 191 191 else 192 192 write(*,*) "phyetat0: <ZMEA> range:", & … … 198 198 if (.not.found) then 199 199 write(*,*) "phyetat0: Failed loading <ZSTD>" 200 call abort 200 call abort_physiq 201 201 else 202 202 write(*,*) "phyetat0: <ZSTD> range:", & … … 208 208 if (.not.found) then 209 209 write(*,*) "phyetat0: Failed loading <ZSIG>" 210 call abort 210 call abort_physiq 211 211 else 212 212 write(*,*) "phyetat0: <ZSIG> range:", & … … 218 218 if (.not.found) then 219 219 write(*,*) "phyetat0: Failed loading <ZGAM>" 220 call abort 220 call abort_physiq 221 221 else 222 222 write(*,*) "phyetat0: <ZGAM> range:", & … … 228 228 if (.not.found) then 229 229 write(*,*) "phyetat0: Failed loading <ZTHE>" 230 call abort 230 call abort_physiq 231 231 else 232 232 write(*,*) "phyetat0: <ZTHE> range:", & … … 238 238 if (.not.found) then 239 239 write(*,*) "phyetat0: Failed loading <tsurf>" 240 call abort 240 call abort_physiq 241 241 else 242 242 write(*,*) "phyetat0: Surface temperature <tsurf> range:", & … … 248 248 if (.not.found) then 249 249 write(*,*) "phyetat0: Failed loading <emis>" 250 call abort 250 call abort_physiq 251 251 else 252 252 write(*,*) "phyetat0: Surface emissivity <emis> range:", & … … 258 258 if (.not.found) then 259 259 write(*,*) "phyetat0: Failed loading <cloudfrac>" 260 call abort 260 call abort_physiq 261 261 else 262 262 write(*,*) "phyetat0: Cloud fraction <cloudfrac> range:", & … … 268 268 if (.not.found) then 269 269 write(*,*) "phyetat0: Failed loading <totcloudfrac>" 270 call abort 270 call abort_physiq 271 271 else 272 272 write(*,*) "phyetat0: Total cloud fraction <totcloudfrac> range:", & … … 278 278 if (.not.found) then 279 279 write(*,*) "phyetat0: Failed loading <hice>" 280 ! call abort 280 ! call abort_physiq 281 281 do ig=1,ngrid 282 282 hice(ig)=0. … … 361 361 if (.not.found) then 362 362 write(*,*) "phyetat0: Failed loading <q2>" 363 call abort 363 call abort_physiq 364 364 else 365 365 write(*,*) "phyetat0: PBL wind variance <q2> range:", & -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90
r310 r313 587 587 if (is_master) write(*,*) "taking dynamics day for physics: ",day_ini 588 588 day_ini=pday 589 !ym stop589 !ym CALL abort_physiq 590 590 endif 591 591 … … 670 670 if (is_master) write(*,*) "physiq: Error! No num_run file!" 671 671 if (is_master) write(*,*) " (which is needed for sourceevol option)" 672 stop672 CALL abort_physiq 673 673 endif 674 674 read(128,*) num_run … … 797 797 print*,'I need values for flatten, J2, Rmean and MassPlanet to compute glat (else set oblate=.false.)' 798 798 799 call abort 799 call abort_physiq 800 800 else 801 801 … … 940 940 if(kastprof)then 941 941 print*,'kastprof should not = true here' 942 call abort 942 call abort_physiq 943 943 endif 944 944 if(water) then … … 1287 1287 if (.not.tracer) then 1288 1288 print*,'We need a CO2 ice tracer to condense CO2' 1289 call abort 1289 call abort_physiq 1290 1290 endif 1291 1291 call condense_cloud(ngrid,nlayer,nq,ptimestep, & … … 1392 1392 call planetwide_sumval(cpp*massarea(:,:)*dtlscale(:,:)/totarea_planet,dEtot) 1393 1393 ! if(isnan(dEtot)) then ! NB: isnan() is not a standard function... 1394 ! print*,'Nan in largescale, abort '1395 ! STOP1394 ! print*,'Nan in largescale, abort_physiq' 1395 ! CALL abort_physiq 1396 1396 ! endif 1397 1397 if (is_master) print*,'In largescale atmospheric energy change =',dEtot,' W m-2' … … 1794 1794 print*,'temp= ',pt(ig,:) 1795 1795 print*,'pplay= ',pplay(ig,:) 1796 call abort 1796 call abort_physiq 1797 1797 endif 1798 1798 end do … … 1844 1844 print*,'As testradtimes enabled,' 1845 1845 print*,'exiting physics on first call' 1846 call abort 1846 call abort_physiq 1847 1847 endif 1848 1848 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/profile.F
r222 r313 184 184 185 185 GOTO 201 186 101 STOP'fichier profile inexistant' 186 101 PRINT*, 'fichier profile inexistant' 187 CALL abort_physiq 187 188 201 CONTINUE 188 189 CLOSE(10) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/radii_mod.F90
r298 r313 93 93 print*,'The code still needs generalisation to arbitrary' 94 94 print*,'aerosol kinds and number.' 95 call abort 95 call abort_physiq 96 96 endif 97 97 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/rain.F90
r227 r313 370 370 if(zrfl(i).lt.0.0)then 371 371 print*,'Droplets of negative rain are falling...' 372 call abort 372 call abort_physiq 373 373 endif 374 374 IF (t(i,1) .LT. T_h2O_ice_liq) THEN -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/setspi.F90
r298 r313 90 90 write(*,*)' datadir = /absolute/path/to/datagcm' 91 91 write(*,*)'Also check that the corrkdir you chose in callphys.def exists.' 92 call abort 92 call abort_physiq 93 93 endif 94 94 … … 110 110 if(nb.ne.L_NSPECTI) then 111 111 write(*,*) 'MISMATCH !! I stop here' 112 call abort 112 call abort_physiq 113 113 endif 114 114 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/setspv.F90
r298 r313 68 68 write(*,*)' datadir = /absolute/path/to/datagcm' 69 69 write(*,*)'Also check that the corrkdir you chose in callphys.def exists.' 70 call abort 70 call abort_physiq 71 71 endif 72 72 … … 87 87 if(nb.ne.L_NSPECTV) then 88 88 write(*,*) 'MISMATCH !! I stop here' 89 call abort 89 call abort_physiq 90 90 endif 91 91 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/soil_settings.F
r222 r313 102 102 if (ierr.ne.0) then 103 103 write(*,*) 'soil_settings: failed allocation of oldmlayer!' 104 stop104 CALL abort_physiq 105 105 endif 106 106 endif … … 171 171 if (.not.found) then 172 172 write(*,*) "soil_settings: Failed loading <inertiedat>" 173 call abort 173 call abort_physiq 174 174 endif 175 175 … … 188 188 write(*,*) 'soil_settings: failed allocation of ', 189 189 & 'oldinertiedat!' 190 stop190 CALL abort_physiq 191 191 endif 192 192 endif ! of if (.not.allocated(oldinertiedat)) … … 194 194 if (.not.found) then 195 195 write(*,*) "soil_settings: Failed loading <inertiedat>" 196 call abort 196 call abort_physiq 197 197 endif 198 198 else ! put values in therm_i … … 200 200 if (.not.found) then 201 201 write(*,*) "soil_settings: Failed loading <inertiedat>" 202 call abort 202 call abort_physiq 203 203 endif 204 204 ! endif … … 224 224 write(*,*) 'soil_settings: failed allocation of ', 225 225 & 'oldtsoil!' 226 stop226 CALL abort_physiq 227 227 endif 228 228 endif … … 230 230 if (.not.found) then 231 231 write(*,*) "soil_settings: Failed loading <tsoil>" 232 call abort 232 call abort_physiq 233 233 endif 234 234 else ! put values in tsoil … … 236 236 if (.not.found) then 237 237 write(*,*) "soil_settings: Failed loading <tsoil>" 238 call abort 238 call abort_physiq 239 239 endif 240 240 endif ! of if (interpol) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/start2archive.F
r227 r313 176 176 IF (ierr.NE.NF_NOERR) THEN 177 177 write(6,*)' Pb d''ouverture du fichier'//trim(fichnom) 178 CALL ABORT178 CALL abort_physiq 179 179 ENDIF 180 180 … … 182 182 IF (ierr .NE. NF_NOERR) THEN 183 183 PRINT*, "start2archive: Le champ <controle> est absent" 184 CALL abort 184 CALL abort_physiq 185 185 ENDIF 186 186 #ifdef NC_DOUBLE … … 191 191 IF (ierr .NE. NF_NOERR) THEN 192 192 PRINT*, "start2archive: Lecture echoue pour <controle>" 193 CALL abort 193 CALL abort_physiq 194 194 ENDIF 195 195 … … 217 217 IF (ierr.NE.NF_NOERR) THEN 218 218 write(6,*)' Pb d''ouverture du fichier'//trim(fichnom) 219 CALL ABORT219 CALL abort_physiq 220 220 ENDIF 221 221 … … 223 223 IF (ierr .NE. NF_NOERR) THEN 224 224 PRINT*, "start2archive: Le champ <controle> est absent" 225 CALL abort 225 CALL abort_physiq 226 226 ENDIF 227 227 #ifdef NC_DOUBLE … … 232 232 IF (ierr .NE. NF_NOERR) THEN 233 233 PRINT*, "start2archive: Lecture echoue pour <controle>" 234 CALL abort 234 CALL abort_physiq 235 235 ENDIF 236 236 … … 242 242 c----------------------------------------------------------------------- 243 243 !mars a voir if ((day_ini_fi.ne.day_ini).or.(abs(timefi-timedyn).gt.1.e-10)) 244 if ((day_ini_fi.ne.day_ini)) 245 & stop ' Probleme de Synchro entre start et startfi !!!' 244 if ((day_ini_fi.ne.day_ini)) THEN 245 PRINT*, 'Probleme de Synchro entre start et startfi !!!' 246 CALL abort_physiq 247 endif 246 248 247 249 … … 401 403 if (ierr.ne.NF_NOERR) then 402 404 write(*,*) "time matter ",NF_STRERROR(ierr) 403 stop405 CALL abort_physiq 404 406 endif 405 407 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/su_gases.F90
r227 r313 30 30 write(*,*) "sugases.F90: error reading number of gases" 31 31 write(*,*) " (first line of gases.def) " 32 call abort 32 call abort_physiq 33 33 endif 34 34 … … 40 40 if (ierr.ne.0) then 41 41 write(*,*) 'sugases.F90: error reading gas names in gases.def...' 42 call abort 42 call abort_physiq 43 43 endif 44 44 enddo !of do igas=1,ngasmx … … 50 50 if (ierr.ne.0) then 51 51 write(*,*) 'sugases.F90: error reading gas molar fractions in gases.def...' 52 call abort 52 call abort_physiq 53 53 endif 54 54 … … 60 60 print*,'You seem to be choosing two variable gases' 61 61 print*,'Check that gases.def is correct' 62 call abort 62 call abort_physiq 63 63 endif 64 64 endif … … 120 120 else 121 121 write(*,*) 'Cannot find required file "gases.def"' 122 call abort 122 call abort_physiq 123 123 endif 124 124 close(90) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/suaer_corrk.F90
r298 r313 222 222 write(*,*)' http://www.lmd.jussieu.fr/',& 223 223 '~forget/datagcm/datafile' 224 CALL ABORT224 CALL abort_physiq 225 225 ENDIF 226 226 OPEN(UNIT=file_unit,& … … 248 248 WRITE(*,*) 'readoptprop: ',& 249 249 'Error while loading optical properties.' 250 CALL ABORT250 CALL abort_physiq 251 251 END SELECT reading1_seq ! ============================== 252 252 ENDIF … … 323 323 WRITE(*,*) 'readoptprop: ',& 324 324 'Error while loading optical properties.' 325 CALL ABORT325 CALL abort_physiq 326 326 END SELECT reading2_seq ! ============================== 327 327 ENDIF -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/tabfi.F
r298 r313 164 164 if (.not.found) then 165 165 write(*,*)"tabfi: Failed reading <controle> array" 166 call abort 166 call abort_physiq 167 167 else 168 168 if (is_master) write(*,*)'tabfi: tab_cntrl',tab_cntrl … … 276 276 write(*,*) "tabfi: Error modifying tab_control should", 277 277 & " only happen in serial mode (eg: by newstart)" 278 stop278 CALL abort_physiq 279 279 endif 280 280 c----------------------------------------------------------------------- -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/vlz_fi.F
r227 r313 169 169 c wq(ij,l+1)= (MQtot + (-w(ij,l+1)-Mtot)*qm(ij,1)) 170 170 write(*,*) 'a rather weird situation in vlz_fi !' 171 stop171 CALL abort_physiq 172 172 end if 173 173 endif -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/write_archive.F
r222 r313 115 115 write(*,*) "***** PUT_VAR matter in write_archive" 116 116 write(*,*) "***** with ",nom," ",nf_STRERROR(ierr) 117 call abort 117 call abort_physiq 118 118 endif 119 119 … … 134 134 write(*,*)"write_archive: dimension <subsurface_layers>", 135 135 & " is missing !!!" 136 call abort 136 call abort_physiq 137 137 endif 138 138 ierr=NF_INQ_DIMID(nid,"Time",id(4)) … … 202 202 write(*,*) "***** PUT_VAR matter in write_archive" 203 203 write(*,*) "***** with ",nom,nf_STRERROR(ierr) 204 call abort 204 call abort_physiq 205 205 endif 206 206 … … 238 238 write(*,*) "***** PUT_VAR matter in write_archive" 239 239 write(*,*) "***** with ",nom,nf_STRERROR(ierr) 240 call abort 240 call abort_physiq 241 241 endif 242 242 243 243 else 244 244 write(*,*) "write_archive: dim=",dim," ?!?" 245 call abort 245 call abort_physiq 246 246 endif ! of if (dim.eq.3) else if (dim.eq.-3) .... 247 247 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/writediagfi.F
r253 r313 150 150 if (n.ge.n_nom_def_max) then 151 151 write(*,*)"n_nom_def_max too small in writediagfi.F:",n 152 stop152 CALL abort_physiq 153 153 end if 154 154 n_nom_def=n-1 … … 183 183 write(*,*) " firstnom string not long enough!!" 184 184 write(*,*) " increase its size to at least ",len_trim(nom) 185 stop185 CALL abort_physiq 186 186 endif 187 187 … … 457 457 write(*,*) "writediagfi error: dim=1 not implemented ", 458 458 & "in parallel mode" 459 stop459 CALL abort_physiq 460 460 endif 461 461 ! Passage variable physique --> physique dynamique -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/writediagsoil.F90
r245 r313 77 77 write(*,*) " firstname string not long enough!!" 78 78 write(*,*) " increase its size to at least ",len_trim(name) 79 stop79 CALL abort_physiq 80 80 endif 81 81 … … 89 89 if (ierr.ne.NF_NOERR) then 90 90 write(*,*)'writediagsoil: Error, failed creating file '//trim(filename) 91 stop91 CALL abort_physiq 92 92 endif 93 93 endif ! of if (is_master) … … 133 133 if (ierr.ne.NF_NOERR) then 134 134 write(*,*)"writediagsoil: Failed writing date to time variable" 135 stop135 CALL abort_physiq 136 136 endif 137 137 endif ! of if (is_master) … … 285 285 #ifdef CPP_PARA 286 286 write(*,*) "writediagsoil: dimps==0 case not implemented in // mode!!" 287 stop287 CALL abort_physiq 288 288 #endif 289 289 ! A. Copy data value -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/writediagspecIR.F
r245 r313 138 138 write(*,*) " firstnom string not long enough!!" 139 139 write(*,*) " increase its size to at least ",len_trim(nom) 140 stop140 CALL abort_physiq 141 141 endif 142 142 … … 214 214 write(*,*) "***** with time" 215 215 write(*,*) 'ierr=', ierr 216 c call abort 216 c call abort_physiq 217 217 endif 218 218 … … 295 295 write(*,*) "***** with ",nom 296 296 write(*,*) 'ierr=', ierr 297 call abort 297 call abort_physiq 298 298 endif 299 299 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/writediagspecVI.F
r245 r313 137 137 write(*,*) " firstnom string not long enough!!" 138 138 write(*,*) " increase its size to at least ",len_trim(nom) 139 stop139 CALL abort_physiq 140 140 endif 141 141 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/wstats.F90
r245 r313 222 222 if (ierr.ne.NF_NOERR) then 223 223 write (*,*) NF_STRERROR(ierr) 224 stop ""224 CALL abort_physiq 225 225 endif 226 226 … … 237 237 if (ierr.ne.NF_NOERR) then 238 238 write (*,*) NF_STRERROR(ierr) 239 stop ""239 CALL abort_physiq 240 240 endif 241 241 endif … … 405 405 write(*,*) "def_var_stats: Failed defining variable "//trim(name) 406 406 write(*,*) NF_STRERROR(ierr) 407 stop ""407 CALL abort_physiq 408 408 endif 409 409 … … 414 414 write(*,*) "def_var_stats: Failed writing title attribute for "//trim(name) 415 415 write(*,*) NF_STRERROR(ierr) 416 stop ""416 CALL abort_physiq 417 417 endif 418 418 … … 422 422 write(*,*) "def_var_stats: Failed writing units attribute for "//trim(name) 423 423 write(*,*) NF_STRERROR(ierr) 424 stop ""424 CALL abort_physiq 425 425 endif 426 426 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/xios_output.F90
r270 r313 14 14 CONTAINS 15 15 16 SUBROUTINE initialize_xios_output 16 SUBROUTINE initialize_xios_output(time) 17 17 USE comgeomphy 18 18 USE xios … … 24 24 IMPLICIT NONE 25 25 INCLUDE "comcstfi.h" 26 26 INTEGER,INTENT(IN) :: time 27 27 TYPE(xios_context) :: ctx_hdl 28 28 TYPE(xios_time) :: dtime … … 35 35 INTEGER :: l 36 36 REAL :: presnivs(klev) 37 INTEGER :: time0 38 CHARACTER(LEN=256) :: start_date 39 INTEGER :: year,month,day,hour,minute,second 37 40 38 41 CALL gather_omp(rlond,lon_mpi) … … 43 46 !$OMP BARRIER 44 47 !$OMP MASTER 48 time0=time 49 year=time0/(86400*360) 50 time0=time0-year*(86400*360) 51 month=time0/(86400*30) 52 time0=time0-month*(86400*30) 53 month=month+1 54 day=time0/86400 55 time0=time0-day*86400 56 day=day+1 57 hour=time0/3600 58 time0=time0-hour*3600 59 minute=time0/60 60 time0=time0-minute*60 61 second=time0 62 write(start_date,'(i0.4,"-",i0.2,"-",i0.2," ",i0.2,":",i0.2,":",i0.2)'),year,month,day,hour,minute,second 63 64 PRINT *,"time0 ",time," Start Date ",TRIM(start_date) 45 65 CALL xios_context_initialize(context_id,comm_lmdz_phy) 46 66 CALL xios_get_handle(context_id, ctx_hdl) 47 67 CALL xios_set_current_context(ctx_hdl) 68 CALL xios_set_attr(ctx_hdl,start_date=TRIM(start_date)) 48 69 49 70 ! lev_value(:) = (/ (l,l=1,klev) /)
Note: See TracChangeset
for help on using the changeset viewer.