c file: 'dbgdmp.for' subroutine dbgdmp(day,sr) C ****************************************************************** wjr C The dumps variables that have gone out of range C C EDIT HISTORY C 01-Mar-99 wjr original coding C integer day integer sr C c + + + GLOBAL COMMON BLOCKS + + + *$noereference include 'p1werm.inc' include 'p1unconv.inc' include 'wpath.inc' include 'm1subr.inc' include 'm1sim.inc' include 'm1geo.inc' include 'm1flag.inc' include 'm1dbug.inc' include 's1layr.inc' include 's1surf.inc' include 's1phys.inc' include 's1agg.inc' include 's1dbh.inc' include 's1dbc.inc' include 's1layd.inc' include 's1sgeo.inc' include 'c1info.inc' include 'c1gen.inc' include 'c1db1.inc' include 'c1db2.inc' include 'c1db3.inc' include 'c1db4.inc' include 'c1db5.inc' include 'c1glob.inc' include 'd1glob.inc' include 'b1glob.inc' include 'w1clig.inc' include 'w1wind.inc' include 'w1pavg.inc' include 'h1hydro.inc' include 'h1scs.inc' include 'h1db1.inc' include 'h1temp.inc' c + + + LOCAL COMMON BLOCKS + + + include 'main/main.inc' include 'main/output.inc' include 'manage/man.inc' include 'erosion/m2geo.inc' include 'erosion/e2erod.inc' *$INCLUDE fsublib.fi *$reference C integer idx,jdx C real tstmin parameter (tstmin=1e-10) C real tstmax parameter (tstmax=1e10) C logical dmpflg data dmpflg /.true./ C prototype, remove before compiling C if (#(sr).lt.tstmin.or.#(sr).gt.tstmax) C * write(*,*) 'day ',day,' # ', #(sr) C C C s1surf C if (dmpflg) write(*,*) 's1surf' C if (aszcr(sr).lt.0.0.or.aszcr(sr).gt.23.0) * write(*,*) 'day ',day,' aszcr ', aszcr(sr) C if (asfcr(sr).lt.0.0.or.asfcr(sr).gt.1.0) * write(*,*) 'day ',day,' asfcr ', asfcr(sr) C if (asmlos(sr).lt.0.0.or.asmlos(sr).gt.2.0) * write(*,*) 'day ',day,' asmlos ', asmlos(sr) C if (asflos(sr).lt.0.0.or.asflos(sr).gt.1.0) * write(*,*) 'day ',day,' asflos ', asflos(sr) C C wjr, test values based on definition if (asdcr(sr).lt.0.6.or.asdcr(sr).gt.2.0) * write(*,*) 'day ',day,' asdcr ', asdcr(sr) C if (asecr(sr).lt.0.1.or.asecr(sr).gt.7.0) * write(*,*) 'day ',day,' asecr ', asecr(sr) C if (asfald(sr).lt.0.05.or.asfald(sr).gt.0.25) * write(*,*) 'day ',day,' asfald ', asfald(sr) C if (asfalw(sr).lt.0.05.or.asfalw(sr).gt.0.2) * write(*,*) 'day ',day,' asfalw ', asfalw(sr) C C w1info C C *** if (aw0cln.lt.70.0.or.aw0cln.gt.170.0) C *** * write(*,*) 'day ',day,' aw0cln ', aw0cln C C *** if (aw0clt.lt.15.0.or.aw0clt.gt.75.0) C *** * write(*,*) 'day ',day,' aw0clt ', aw0clt C C *** if (aw0wln.lt.70.0.or.aw0wln.gt.170.0) C *** * write(*,*) 'day ',day,' aw0wln ', aw0wln C C *** if (aw0wlt.lt.15.0.or.aw0wlt.gt.75.0) C *** * write(*,*) 'day ',day,' aw0wlt ', aw0wlt C C s1sgeo C if (dmpflg) write(*,*) 's1sgeo' C if (aszrgh(sr).lt.0.0.or.aszrgh(sr).gt.500.0) * write(*,*) 'day ',day,' aszrgh ', aszrgh(sr) C if (asxrgw(sr).lt.10.0.or.asxrgw(sr).gt.4000.0) * write(*,*) 'day ',day,' asxrgw ', asxrgw(sr) C if (asxrgs(sr).lt.10.0.or.asxrgs(sr).gt.2000.0) * write(*,*) 'day ',day,' asxrgs ', asxrgs(sr) C if (asargo(sr).lt.0.0.or.asargo(sr).gt.179.0) * write(*,*) 'day ',day,' asargo ', asargo(sr) C C wjr, test values based on definition if (asxdks(sr).lt.0.0.or.asxdks(sr).gt.1000.0) * write(*,*) 'day ',day,' asxdks ', asxdks(sr) C C wjr, test values based on definition if (asxdkh(sr).lt.0.0.or.asxdkh(sr).gt.1000.0) * write(*,*) 'day ',day,' asxdkh ', asxdkh(sr) C if (as0rrk(sr).lt.tstmin.or.as0rrk(sr).gt.tstmax) * write(*,*) 'day ',day,' as0rrk ', as0rrk(sr) C if (aslrrc(sr).lt.tstmin.or.aslrrc(sr).gt.tstmax) * write(*,*) 'day ',day,' aslrrc ', aslrrc(sr) C if (aslrr(sr).lt.1.0.or.aslrr(sr).gt.30.0) * write(*,*) 'day ',day,' aslrr ', aslrr(sr) C C w1wind C if (dmpflg) write(*,*) 'w1wind' C C wjr, test values based on definition if (awadir.lt.0.0.or.awadir.gt.360.0) * write(*,*) 'day ',day,' awadir ', awadir C if (awhrmx.lt.1.0.or.awhrmx.gt.24.0) * write(*,*) 'day ',day,' awhrmx ', awhrmx C if (awrmxn.lt.0.0.or.awrmxn.gt.tstmax) * write(*,*) 'day ',day,' awrmxn ', awrmxn C C wjr, test values based on definition if (awudmx.lt.0.0.or.awudmx.gt.50.0) * write(*,*) 'day ',day,' awudmx ', awudmx C C wjr, test values based on definition if (awudmn.lt.0.0.or.awudmn.gt.25.0) * write(*,*) 'day ',day,' awudmn ', awudmn C C wjr, test values based on definition if (awudav.lt.0.0.or.awudav.gt.35.0) * write(*,*) 'day ',day,' awudav ', awudav C do 10 idx=1,mntime C wjr, test values based on definition if (awu(idx).lt.0.0.or.awu(idx).gt.35.0) * write(*,*) 'day ',day,' awu(',idx,') ', awu(idx) 10 continue C C w1pagv C C wjr, test values based on definition if (awtmmx.lt.-10.0.or.awtmmx.gt.40.0) * write(*,*) 'day ',day,' awtmmx ', awtmmx C C wjr, test values based on definition if (awtmmn.lt.-20.0.or.awtmmn.gt.30.0) * write(*,*) 'day ',day,' awtmmn ', awtmmn C C wjr, test values based on definition if (awdair.lt.0.0.or.awdair.gt.tstmax) * write(*,*) 'day ',day,' awdair ', awdair C C wjr, test values based on definition if (awztpt.lt.0.0.or.awztpt.gt.500.0) * write(*,*) 'day ',day,' awztpt ', awztpt C C wjr, test values based on definition if (awtpav.lt.-20.0.or.awtpav.gt.40) * write(*,*) 'day ',day,' awtpav ', awtpav C C *** if (awepir.lt.tstmin.or.awepir.gt.tstmax) C *** * write(*,*) 'day ',day,' awepir ', awepir C C wjr, test values based on definition if (awupav.lt.0.0.or.awupav.gt.30.0) * write(*,*) 'day ',day,' awupav ', awupav C C wjr, test values based on definition if (awnuet.lt.0.or.awnuet.gt.31.0) * write(*,*) 'day ',day,' awnuet ', awnuet C C wjr, test values based on definition if (aweuet.lt.0.0.or.aweuet.gt.tstmax) * write(*,*) 'day ',day,' aweuet ', aweuet C C b1geom C if (dmpflg) write(*,*) 'b1geom' C C wjr, test values based on definition if (abrsai(sr).lt.0.0.or.abrsai(sr).gt.1.0) * write(*,*) 'day ',day,' abrsai ', abrsai(sr) C C wjr, test values based on definition if (abrlai(sr).lt.0.0.or.abrlai(sr).gt.1.0) * write(*,*) 'day ',day,' abrlai ', abrlai(sr) C do 20 idx=1,mncz C C wjr, test values based on definition if (abrsaz(idx, sr).lt.0.0.or.abrsaz(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' abrsaz(',idx,') ', abrsaz(idx, sr) C C wjr, test values based on definition if (abrlaz(idx, sr).lt.0.0.or.abrlaz(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' abrlaz(',idx,') ', abrlaz(idx, sr) C 20 continue C C wjr, test values based on definition if (abffcv(sr).lt.0.0.or.abffcv(sr).gt.1.0) * write(*,*) 'day ',day,' abffcv ', abffcv(sr) C C wjr, test values based on definition if (abfscv(sr).lt.0.0.or.abfscv(sr).gt.1.0) * write(*,*) 'day ',day,' abfscv ', abfscv(sr) C C wjr, test values based on definition if (abftcv(sr).lt.0.0.or.abftcv(sr).gt.1.0) * write(*,*) 'day ',day,' abftcv ', abftcv(sr) C C w1clig C if (dmpflg) write(*,*) 'w1clig' C C wjr, test values based on definition if (awrrh.lt.0.0.or.awrrh.gt.100.0) * write(*,*) 'day ',day,' awrrh ', awrrh C C wjr, test values based on definition if (awtdav.lt.-20.0.or.awtdav.gt.50.0) * write(*,*) 'day ',day,' awtdav ', awtdav C C wjr, test values based on definition if (awtyav.lt.0.0.or.awtyav.gt.30.0) * write(*,*) 'day ',day,' awtyav ', awtyav C do 30 idx=1,12 C wjr, test values based on definition if (awtmav(idx).lt.-10.0.or.awtmav(idx).gt.40.0) * write(*,*) 'day ',day,' awtmav(',idx,') ', awtmav(idx) 30 continue C C wjr, test values based on definition if (awtdmx.lt.0.0.or.awtdmx.gt.50.0) * write(*,*) 'day ',day,' awtdmx ', awtdmx C C wjr, test values based on definition if (awtdmn.lt.-20.0.or.awtdmn.gt.40.0) * write(*,*) 'day ',day,' awtdmn ', awtdmn C C wjr, test values based on definition if (awtdpt.lt.0.0.or.awtdpt.gt.40.0) * write(*,*) 'day ',day,' awtdpt ', awtdpt C C wjr, test values based on definition if (awzdpt.lt.0.0.or.awzdpt.gt.1000.0) * write(*,*) 'day ',day,' awzdpt ', awzdpt C C wjr, test values based on definition if (aweirr.lt.0.0.or.aweirr.gt.tstmax) * write(*,*) 'day ',day,' aweirr ', aweirr C C s1psd C C *** do 40 idx=1,mnsz C *** if (aslsgm(idx, sr).lt.0.0.or.aslsgm(idx, sr).gt.10.0) C *** * write(*,*) 'day ',day,' aslsgm ', aslsgm(idx, sr) C C *** if (as0sgs(idx, sr).lt.0.0.or.as0sgs(idx, sr).gt.10.0) C *** * write(*,*) 'day ',day,' as0sgs ', as0sgs(idx, sr) C *** 40 continue C C s1layd C if (dmpflg) write(*,*) 's1layd' C do 50 idx=1,nslay(sr) if (asdsbk(idx, sr).lt.tstmin.or.asdsbk(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asdsbk(',idx,') ', asdsbk(idx, sr) C if (aszlyd(idx, sr).lt.tstmin.or.aszlyd(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' aszlyd(',idx,') ', aszlyd(idx, sr) C C *** if (aszlym(idx, sr).lt.tstmin.or.aszlym(idx, sr).gt.tstmax) C *** * write(*,*) 'day ',day,' aszlym(',idx,') ', aszlym(idx, sr) C C *** if (aszmpt(idx, sr).lt.tstmin.or.aszmpt(idx, sr).gt.tstmax) C *** * write(*,*) 'day ',day,' aszmpt(',idx,') ', aszmpt(idx, sr) 50 continue C C s1layr C if (dmpflg) write(*,*) 's1layr' C if (nslay(sr).lt.1.or.nslay(sr).gt.10) * write(*,*) 'day ',day,' nslay ', nslay(sr) C if (aszlyt(1, sr).lt.10.0.or.aszlyt(1, sr).gt.10.0) * write(*,*) 'day ',day,' aszlyt(1) ', aszlyt(1, sr) C if (nslay(sr).gt.1.and. * (aszlyt(2, sr).lt.40.0.or.aszlyt(2, sr).gt.40.0)) * write(*,*) 'day ',day,' aszlyt(2) ', aszlyt(2, sr) C if (nslay(sr).gt.2.and. * (aszlyt(3, sr).lt.50.0.or.aszlyt(3, sr).gt.100.0)) * write(*,*) 'day ',day,' aszlyt(3) ', aszlyt(3, sr) C if (nslay(sr).gt.3.and. * (aszlyt(4, sr).lt.50.0.or.aszlyt(4, sr).gt.100.0)) * write(*,*) 'day ',day,' aszlyt(4) ', aszlyt(4, sr) C do 60 idx=5,mnsz+1 if (nslay(sr).ge.idx.and. * (aszlyt(idx, sr).lt.1.0.or.aszlyt(idx, sr).gt.1000.0)) * write(*,*) 'day ',day,' aszlyt(',idx,') ', aszlyt(idx, sr) 60 continue C C s1phys C if (dmpflg) write(*,*) 's1phys' C do 70 idx=0,mnsz if (asdblk(idx, sr).lt.0.50.or.asdblk(idx, sr).gt.2.5) * write(*,*) 'day ',day,' asdblk(',idx,') ', asdblk(idx, sr) 70 continue C C s1dbh C if (dmpflg) write(*,*) 's1dbh' C do 80 idx=0,mnsz if (asfsan(idx, sr).lt.0.0.or.asfsan(idx, sr).gt.1.0) * write(*,*) 'day ',day,' asfsan(',idx,') ', asfsan(idx, sr) C if (asfsil(idx, sr).lt.0.0.or.asfsil(idx, sr).gt.1.0) * write(*,*) 'day ',day,' asfsil(',idx,') ', asfsil(idx, sr) C if (asfcla(idx, sr).lt.0.0.or.asfcla(idx, sr).gt.1.0) * write(*,*) 'day ',day,' asfcla(',idx,') ', asfcla(idx, sr) C if (asvroc(idx, sr).lt.0.0.or.asvroc(idx, sr).gt.1.0) * write(*,*) 'day ',day,' asvroc(',idx,') ', asvroc(idx, sr) 80 continue C C s1agg C if (dmpflg) write(*,*) 's1agg' C do 90 idx=0,mnsz if (asdagd(idx, sr).lt.0.6.or.asdagd(idx, sr).gt.2.5) * write(*,*) 'day ',day,' asdagd(',idx,') ', asdagd(idx, sr) C if (aseags(idx, sr).lt.0.1.or.aseags(idx, sr).gt.7.0) * write(*,*) 'day ',day,' aseags(',idx,') ', aseags(idx, sr) C if (aslagm(idx, sr).lt.0.03.or.aslagm(idx, sr).gt.30.0) * write(*,*) 'day ',day,' aslagm(',idx,') ', aslagm(idx, sr) C if (aslagn(idx, sr).lt.0.001.or.aslagn(idx, sr).gt.5.0) * write(*,*) 'day ',day,' aslagn(',idx,') ', aslagn(idx, sr) C if (aslagx(idx, sr).lt.1.0.or.aslagx(idx, sr).gt.1000.0) * write(*,*) 'day ',day,' aslagx(',idx,') ', aslagx(idx, sr) C if (as0ags(idx, sr).lt.1.0.or.as0ags(idx, sr).gt.20.0) * write(*,*) 'day ',day,' as0ags(',idx,') ', as0ags(idx, sr) 90 continue C C s1dbc C if (dmpflg) write(*,*) 's1dbc' C if (asfom(0, sr).lt.tstmin.or.asfom(0, sr).gt.tstmax) * write(*,*) 'day ',day,' asfom(0) ', asfom(0, sr) C do 100 idx=1,mnsz if (as0ph(idx, sr).lt.0.0.or.as0ph(idx, sr).gt.14.0) * write(*,*) 'day ',day,' as0ph(',idx,') ', as0ph(idx, sr) C if (ascmg(idx, sr).lt.0.0.or.ascmg(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ascmg(',idx,') ', ascmg(idx, sr) C if (ascna(idx, sr).lt.0.0.or.ascna(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ascna(',idx,') ', ascna(idx, sr) C if (asfcce(idx, sr).lt.0.0.or.asfcce(idx, sr).gt.100.0) * write(*,*) 'day ',day,' asfcce(',idx,') ', asfcce(idx, sr) C if (asfcec(idx, sr).lt.0.0.or.asfcec(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asfcec(',idx,') ', asfcec(idx, sr) C if (asfesp(idx, sr).lt.0.0.or.asfesp(idx, sr).gt.100.0) * write(*,*) 'day ',day,' asfesp(',idx,') ', asfesp(idx, sr) C if (asfom(idx, sr).lt.0.0.or.asfom(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asfom(',idx,') ', asfom(idx, sr) C if (asfnoh(idx, sr).lt.0.0.or.asfnoh(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asfnoh(',idx,') ', asfnoh(idx, sr) C if (asfpoh(idx, sr).lt.0.0.or.asfpoh(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asfpoh(',idx,') ', asfpoh(idx, sr) C if (asfpsp(idx, sr).lt.0.0.or.asfpsp(idx, sr).gt.1.0) * write(*,*) 'day ',day,' asfpsp(',idx,') ', asfpsp(idx, sr) C if (asfsmb(idx, sr).lt.0.0.or.asfsmb(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asfsmb(',idx,') ', asfsmb(idx, sr) C if (asftap(idx, sr).lt.0.0.or.asftap(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asftap(',idx,') ', asftap(idx, sr) C if (asftan(idx, sr).lt.0.0.or.asftan(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' asftan(',idx,') ', asftan(idx, sr) 100 continue C if (asmno3(sr).lt.0.0.or.asmno3(sr).gt.tstmax) * write(*,*) 'day ',day,' asmno3 ', asmno3(sr) C C m1sim C if (dmpflg) write(*,*) 'm1sim' C if (ntstep.lt.1.or.ntstep.gt.96) * write(*,*) 'day ',day,' ntstep ', ntstep C if (am0jd.ne.day) * write(*,*) 'day ',day,' am0jd ', am0jd C if (amalat.lt.15.0.or.amalat.gt.75.0) * write(*,*) 'day ',day,' amalat ', amalat C if (amalon.lt.70.0.or.amalon.gt.170.0) * write(*,*) 'day ',day,' amalon ', amalon C if (amzele.lt.0.0.or.amzele.gt.2500.0) * write(*,*) 'day ',day,' amzele ', amzele C C m1subr C if (dmpflg) write(*,*) 'm1subr' C if (nsubr.lt.1.or.nsubr.gt.4) * write(*,*) 'day ',day,' nsubr ', nsubr C if (am0csr.lt.1.or.am0csr.gt.4) * write(*,*) 'day ',day,' am0csr ', am0csr C if (am0irr(sr).lt.0.or.am0irr(sr).gt.2) * write(*,*) 'day ',day,' am0irr ', am0irr(sr) C if (amnryr(sr).lt.1.or.amnryr(sr).gt.10) * write(*,*) 'day ',day,' amnryr ', amnryr(sr) C if (amrslp(sr).lt.0.0.or.amrslp(sr).gt.1.0) * write(*,*) 'day ',day,' amrslp ', amrslp(sr) C C h1temp C if (dmpflg) write(*,*) 'h1temp' C do 110 idx=1,mnsz if (ahtsav(idx, sr).lt.-20.0.or.ahtsav(idx, sr).gt.50.0) * write(*,*) 'day ',day,' ahtsav(',idx,') ', ahtsav(idx, sr) C if (ahtsmx(idx, sr).lt.-20.0.or.ahtsmx(idx, sr).gt.50.0) * write(*,*) 'day ',day,' ahtsmx(',idx,') ', ahtsmx(idx, sr) C if (ahtsmn(idx, sr).lt.-20.0.or.ahtsmn(idx, sr).gt.50.0) * write(*,*) 'day ',day,' ahtsmn(',idx,') ', ahtsmn(idx, sr) 110 continue C C h1et C C *** if (ahzea.lt.0.0.or.ahzea.gt.50.0) C *** * write(*,*) 'day ',day,' ahzea ', ahzea C C *** if (ahzep.lt.0.0.or.ahzep.gt.50.0) C *** * write(*,*) 'day ',day,' ahzep ', ahzep C C *** if (ahzeta.lt.0.0.or.ahzeta.gt.50.0) C *** * write(*,*) 'day ',day,' ahzeta ', ahzeta C C *** if (ahzetp.lt.0.0.or.ahzetp.gt.50.0) C *** * write(*,*) 'day ',day,' ahzetp ', ahzetp C C *** if (ahzpta.lt.0.0.or.ahzpta.gt.50.0) C *** * write(*,*) 'day ',day,' ahzpta ', ahzpta C C *** if (ahzptp.lt.0.0.or.ahzptp.gt.50.0) C *** * write(*,*) 'day ',day,' ahzptp ', ahzptp C C *** if (ah0drat.lt.0.0.or.ah0drat.gt.1.0) C *** * write(*,*) 'day ',day,' ah0drat ', ah0drat C C h1hydro C if (dmpflg) write(*,*) 'h1hydro' C do 120 idx=1,mnsz if (ahrwc(idx, sr).lt.0.011.or.ahrwc(idx, sr).gt.0.379) * write(*,*) 'day ',day,' ahrwc(',idx,') ', ahrwc(idx, sr) C if (aheaep(idx, sr).lt.-17.91.or.aheaep(idx, sr).gt.0.0) * write(*,*) 'day ',day,' aheaep(',idx,') ', aheaep(idx, sr) C if (ahrsk(idx, sr).lt.0.0.or.ahrsk(idx, sr).gt.0.001) * write(*,*) 'day ',day,' ahrsk(',idx,') ', ahrsk(idx, sr) C if (ah0cb(idx, sr).lt.0.917.or.ah0cb(idx, sr).gt.27.927) * write(*,*) 'day ',day,' ah0cb(',idx,') ', ah0cb(idx, sr) 120 continue C if (ahfwsf(sr).lt.tstmin.or.ahfwsf(sr).gt.tstmax) * write(*,*) 'day ',day,' ahfwsf ', ahfwsf(sr) C if (ahzsno(sr).lt.0.0.or.ahzsno(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzsno ', ahzsno(sr) C if (ahzirr(sr).lt.0.0.or.ahzirr(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzirr ', ahzirr(sr) C if (ahzper(sr).lt.0.0.or.ahzper(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzper ', ahzper(sr) C if (ahzrun(sr).lt.0.0.or.ahzrun(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzrun ', ahzrun(sr) C if (ahzsmt(sr).lt.0.0.or.ahzsmt(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzsmt ', ahzsmt(sr) C C h1scs C if (dmpflg) write(*,*) 'h1scs' C if (ah0cng(sr).lt.6.0.or.ah0cng(sr).gt.91.0) * write(*,*) 'day ',day,' ah0cng ', ah0cng(sr) C if (ah0cnp(sr).lt.45.0.or.ah0cnp(sr).gt.94.0) * write(*,*) 'day ',day,' ah0cnp ', ah0cnp(sr) C C h1db1 C if (dmpflg) write(*,*) 'h1db1' C do 130 idx=1,mnsz if (ahrwcw(idx, sr).lt.0.005.or.ahrwcw(idx, sr).gt.0.242) * write(*,*) 'day ',day,' ahrwcw(',idx,') ', ahrwcw(idx, sr) C if (ahrwcf(idx, sr).lt.0.012.or.ahrwcf(idx, sr).gt.0.335) * write(*,*) 'day ',day,' ahrwcf(',idx,') ', ahrwcf(idx, sr) C if (ahrwcs(idx, sr).lt.0.208.or.ahrwcs(idx, sr).gt.0.440) * write(*,*) 'day ',day,' ahrwcs(',idx,') ', ahrwcs(idx, sr) C if (ahrwca(idx, sr).lt.0.0.or.ahrwca(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ahrwca(',idx,') ', ahrwca(idx, sr) 130 continue C do 140 idx=1,mnhhrs if (ahrwc0(idx, sr).lt.0.0.or.ahrwc0(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ahrwc0(',idx,') ', ahrwc0(idx, sr) 140 continue C if (ahzsnd(sr).lt.0.0.or.ahzsnd(sr).gt.tstmax) * write(*,*) 'day ',day,' ahzsnd ', ahzsnd(sr) C C p1werm C do 145 idx=1,mnbpls if (admf(idx,sr).lt.0.0.or.admf(idx,sr).gt.tstmax) * write(*,*) 'day ',day,' admf(',idx,')', admf(idx,sr) 145 continue C C Not used anymore c do 150 idx=0,mncz c if (adma(idx, sr).lt.0.0.or.adma(idx, sr).gt.tstmax) c * write(*,*) 'day ',day,' adma(',idx,') ', adma(idx, sr) c 150 continue C C if (admbt(sr).lt.0.0.or.admbt(sr).gt.tstmax) C * write(*,*) 'day ',day,' admbt ', admbt(sr) C do 160 idx=1,mnsz do 170 jdx=1,mnbpls if (admbgz(idx, jdx, sr).lt.0.0 .or. * admbgz(idx, jdx, sr).gt.tstmax) * write(*,*) 'day ',day,' admbgz(',idx,jdx,') ', * admbgz(idx, jdx, sr) C if (admrtz(idx, jdx, sr).lt.0.0 .or. * admrtz(idx, jdx, sr).gt.tstmax) * write(*,*) 'day ',day,' admrtz(',idx,jdx,') ', * admrtz(idx, jdx, sr) 170 continue 160 continue C C d1glob C if (dmpflg) write(*,*) 'd1glob' C do 180 idx=1,mnbpls if (admf(idx, sr).lt.tstmin.or.admf(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' admf(',idx,') ', admf(idx, sr) C if (admst(idx, sr).lt.tstmin.or.admst(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' admst(',idx,') ', admst(idx, sr) C if (adzht(idx, sr).lt.tstmin.or.adzht(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' adzht(',idx,') ', adzht(idx, sr) C if (addstm(idx, sr).lt.tstmin.or.addstm(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' addstm(',idx,') ', addstm(idx, sr) 180 continue C C d1gen C C this one I really don't understand C *** if (addpop(sr).lt.tstmin.or.addpop(sr).gt.tstmax) C *** * write(*,*) 'day ',day,' addpop ', addpop(sr) C C c1db3 C if (dmpflg) write(*,*) 'c1db3' C if (ac0bn1(sr).lt.0.0.or.ac0bn1(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bn1 ', ac0bn1(sr) C if (ac0bn2(sr).lt.0.0.or.ac0bn2(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bn2 ', ac0bn2(sr) C if (ac0bn3(sr).lt.0.0.or.ac0bn3(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bn3 ', ac0bn3(sr) C if (ac0bp1(sr).lt.0.0.or.ac0bp1(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bp1 ', ac0bp1(sr) C if (ac0bp2(sr).lt.0.0.or.ac0bp2(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bp2 ', ac0bp2(sr) C if (ac0bp3(sr).lt.0.0.or.ac0bp3(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bp3 ', ac0bp3(sr) C if (acfny(sr).lt.0.0.or.acfny(sr).gt.1.0) * write(*,*) 'day ',day,' acfny ', acfny(sr) C if (acfpy(sr).lt.0.0.or.acfpy(sr).gt.1.0) * write(*,*) 'day ',day,' acfpy ', acfpy(sr) C if (acfwy(sr).lt.0.0.or.acfwy(sr).gt.1.0) * write(*,*) 'day ',day,' acfwy ', acfwy(sr) C C c1db5 C if (dmpflg) write(*,*) 'c1db5' C if (ac0bw1(sr).lt.0.0.or.ac0bw1(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bw1 ', ac0bw1(sr) C if (ac0bw2(sr).lt.0.0.or.ac0bw2(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bw2 ', ac0bw2(sr) C if (ac0bw3(sr).lt.0.0.or.ac0bw3(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bw3 ', ac0bw3(sr) C if (ac0cmn(sr).lt.0.0.or.ac0cmn(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0cmn ', ac0cmn(sr) C C cldb4 C if (dmpflg) write(*,*) 'c1db4' C if (ac0hta(sr).lt.0.0.or.ac0hta(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0hta ', ac0hta(sr) C if (ac0htb(sr).lt.0.0.or.ac0htb(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0htb ', ac0htb(sr) C if (ac0lza(sr).lt.0.0.or.ac0lza(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lza ', ac0lza(sr) C if (ac0lzb(sr).lt.0.0.or.ac0lzb(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lzb ', ac0lzb(sr) C if (ac0saa(sr).lt.0.0.or.ac0saa(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0saa ', ac0saa(sr) C if (ac0sab(sr).lt.0.0.or.ac0sab(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0sab ', ac0sab(sr) C if (ac0lia(sr).lt.0.0.or.ac0lia(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lia ', ac0lia(sr) C if (ac0lib(sr).lt.0.0.or.ac0lib(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lib ', ac0lib(sr) C if (ac0cca(sr).lt.0.0.or.ac0cca(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0cca ', ac0cca(sr) C if (ac0ccb(sr).lt.0.0.or.ac0ccb(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0ccb ', ac0ccb(sr) C C c1gen C if (dmpflg) write(*,*) 'c1gen' C if (ac0rg(sr).lt.0.or.ac0rg(sr).gt.1) * write(*,*) 'day ',day,' ac0rg ', ac0rg(sr) C if (acdpop(sr).lt.0.0.or.acdpop(sr).gt.tstmax) * write(*,*) 'day ',day,' acdpop ', acdpop(sr) C if (acxrow(sr).lt.0.0.or.acxrow(sr).gt.1.0) * write(*,*) 'day ',day,' acxrow ', acxrow(sr) C C c1geom C if (dmpflg) write(*,*) 'c1geom' C if (acrsai(sr).lt.0.0.or.acrsai(sr).gt.tstmax) * write(*,*) 'day ',day,' acrsai ', acrsai(sr) C if (acrlai(sr).lt.0.0.or.acrlai(sr).gt.tstmax) * write(*,*) 'day ',day,' acrlai ', acrlai(sr) C do 191 idx=1,mncz if (acrsaz(idx,sr).lt.0.0.or.acrsaz(idx,sr).gt.tstmax) * write(*,*) 'day ',day,' acrsaz(',idx,') ', acrsaz(idx,sr) C if (acrlaz(idx,sr).lt.0.0.or.acrlaz(idx,sr).gt.tstmax) * write(*,*) 'day ',day,' acrlaz(',idx,') ', acrlaz(idx,sr) 191 continue C if (acffcv(sr).lt.0.0.or.acffcv(sr).gt.tstmax) * write(*,*) 'day ',day,' acffcv ', acffcv(sr) C if (acfscv(sr).lt.0.0.or.acfscv(sr).gt.tstmax) * write(*,*) 'day ',day,' acfscv ', acfscv(sr) C if (acftcv(sr).lt.0.0.or.acftcv(sr).gt.tstmax) * write(*,*) 'day ',day,' acftcv ', acftcv(sr) C C c1glob C if (dmpflg) write(*,*) 'c1glob' C if (aczht(sr).lt.0.0.or.aczht(sr).gt.3.0) * write(*,*) 'day ',day,' aczht ', aczht(sr) C if (aczrtd(sr).lt.0.0.or.aczrtd(sr).gt.3.0) * write(*,*) 'day ',day,' aczrtd ', aczrtd(sr) C if (acmst(sr).lt.0.0.or.acmst(sr).gt.tstmax) * write(*,*) 'day ',day,' acmst ', acmst(sr) C if (acmrt(sr).lt.0.0.or.acmrt(sr).gt.tstmax) * write(*,*) 'day ',day,' acmrt ', acmrt(sr) C if (acdstm(sr).lt.0.0.or.acdstm(sr).gt.tstmax) * write(*,*) 'day ',day,' acdstm ', acdstm(sr) C C b1mass C if (dmpflg) write(*,*) 'b1mass' C C no such variable anymore c if (abmres(sr).lt.0.0.or.abmres(sr).gt.tstmax) c * write(*,*) 'day ',day,' abmres ', abmres(sr) C c No such variable anymore c do 190 idx=1,mncz c if (abma(idx, sr).lt.0.0.or.abma(idx, sr).gt.tstmax) c * write(*,*) 'day ',day,' abma(',idx,') ', abma(idx, sr) C c if (abmb(idx, sr).lt.0.0.or.abmb(idx, sr).gt.tstmax) c * write(*,*) 'day ',day,' abmb(',idx,') ', abmb(idx, sr) c 190 continue C C b1glob C if (dmpflg) write(*,*) 'b1glob' C if (abzht(sr).lt.0.0.or.abzht(sr).gt.3.0) * write(*,*) 'day ',day,' abzht ', abzht(sr) C C no such variable anymore c if (abz(sr).lt.0.0.or.abz(sr).gt.3.0) c * write(*,*) 'day ',day,' abz ', abz(sr) C C No such variable anymore c if (abzrtd(sr).lt.0.0.or.abzrtd(sr).gt.3.0) c * write(*,*) 'day ',day,' abzrtd ', abzrtd(sr) C if (abmst(sr).lt.0.0.or.abmst(sr).gt.tstmax) * write(*,*) 'day ',day,' abmst ', abmst(sr) C if (abmrt(sr).lt.0.0.or.abmrt(sr).gt.tstmax) * write(*,*) 'day ',day,' abmrt ', abmrt(sr) C c no such variable anymore c if (abrahi(sr).lt.0.0.or.abrahi(sr).gt.tstmax) c * write(*,*) 'day ',day,' abrahi ', abrahi(sr) C C cldb2 C if (dmpflg) write(*,*) 'c1db2' C if (actopt(sr).lt.0.0.or.actopt(sr).gt.40.0) * write(*,*) 'day ',day,' actopt ', actopt(sr) C if (actmin(sr).lt.0.0.or.actmin(sr).gt.20.0) * write(*,*) 'day ',day,' actmin ', actmin(sr) C if (acfdla(sr).lt.0.0.or.acfdla(sr).gt.1.0) * write(*,*) 'day ',day,' acfdla ', acfdla(sr) C if (acrdla(sr).lt.0.0.or.acrdla(sr).gt.tstmax) * write(*,*) 'day ',day,' acrdla ', acrdla(sr) C if (ac0caf(sr).lt.0.0.or.ac0caf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0caf ', ac0caf(sr) C if (ac0psf(sr).lt.0.0.or.ac0psf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0psf ', ac0psf(sr) C do 200 idx=1,2 if (ac0pt1(idx, sr).lt.0.0.or.ac0pt1(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0pt1(',idx,') ', ac0pt1(idx, sr) C if (ac0pt2(idx, sr).lt.0.0.or.ac0pt2(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0pt2(',idx,') ', ac0pt2(idx, sr) C if (ac0fd1(idx, sr).lt.0.0.or.ac0fd1(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0fd1(',idx,') ', ac0fd1(idx, sr) C if (ac0fd2(idx, sr).lt.0.0.or.ac0fd2(idx, sr).gt.1.0) * write(*,*) 'day ',day,' ac0fd2(',idx,') ', ac0fd2(idx, sr) 200 continue C if (ac0ck(sr).lt.0.0.or.ac0ck(sr).gt.1.0) * write(*,*) 'day ',day,' ac0ck ', ac0ck(sr) C C c1mass C C *** if (acmyld(sr).lt.tstmin.or.acmyld(sr).gt.tstmax) C *** * write(*,*) 'day ',day,' acmyld ', acmyld(sr) C C *** do 210 idx=1,mnsz C *** if (acmbgr(idx, sr).lt.tstmin.or.acmbgr(idx, sr).gt.tstmax) C *** * write(*,*) 'day ',day,' acmbgr ', acmbgr(idx, sr) C *** 210 continue C C c1db1 C if (dmpflg) write(*,*) 'c1db1' C if (acrcn(sr).lt.0.0.or.acrcn(sr).gt.tstmax) * write(*,*) 'day ',day,' acrcn ', acrcn(sr) C if (ac0id(sr).lt.0.or.ac0id(sr).gt.100) * write(*,*) 'day ',day,' ac0id ', ac0id(sr) C if (acrhi(sr).lt.0.0.or.acrhi(sr).gt.tstmax) * write(*,*) 'day ',day,' acrhi ', acrhi(sr) C if (actdtm(sr).lt.0.0.or.actdtm(sr).gt.tstmax) * write(*,*) 'day ',day,' actdtm ', actdtm(sr) C if (aczmrt(sr).lt.0.0.or.aczmrt(sr).gt.tstmax) * write(*,*) 'day ',day,' aczmrt ', aczmrt(sr) C if (aczmxc(sr).lt.0.0.or.aczmxc(sr).gt.tstmax) * write(*,*) 'day ',day,' aczmxc ', aczmxc(sr) C if (acrmla(sr).lt.0.0.or.acrmla(sr).gt.tstmax) * write(*,*) 'day ',day,' acrmla ', acrmla(sr) C if (acrbe(sr).lt.0.0.or.acrbe(sr).gt.tstmax) * write(*,*) 'day ',day,' acrbe ', acrbe(sr) C if (acrbed(sr).lt.0.0.or.acrbed(sr).gt.tstmax) * write(*,*) 'day ',day,' acrbed ', acrbed(sr) C do 220 idx=1,mncz if (ac0lad(idx, sr).lt.0.0.or.ac0lad(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lad(',idx,') ', ac0lad(idx, sr) C if (ac0sad(idx, sr).lt.0.0.or.ac0sad(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0sad(',idx,') ', ac0sad(idx, sr) 220 continue C if (acehu0(sr).lt.0.0.or.acehu0(sr).gt.tstmax) * write(*,*) 'day ',day,' acehu0 ', acehu0(sr) C if (ac0bev(sr).lt.0.0.or.ac0bev(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bev ', ac0bev(sr) C do 230 idx =1,2 if (ac0be1(idx, sr).lt.0.0.or.ac0be1(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0be1(',idx,') ', ac0be1(idx, sr) C if (ac0be2(idx, sr).lt.0.0.or.ac0be2(idx, sr).gt.tstmax) * write(*,*) 'day ',day,' ac0be2(',idx,') ', ac0be2(idx, sr) 230 continue C if (acrmhi(sr).lt.0.0.or.acrmhi(sr).gt.tstmax) * write(*,*) 'day ',day,' acrmhi ', acrmhi(sr) C if (ac0alf(sr).lt.0.0.or.ac0alf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0alf ', ac0alf(sr) C if (ac0blf(sr).lt.0.0.or.ac0blf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0blf ', ac0blf(sr) C if (ac0clf(sr).lt.0.0.or.ac0clf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0clf ', ac0clf(sr) C if (ac0dlf(sr).lt.0.0.or.ac0dlf(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0dlf ', ac0dlf(sr) C if (ac0arp(sr).lt.0.0.or.ac0arp(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0arp ', ac0arp(sr) C if (ac0brp(sr).lt.0.0.or.ac0brp(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0brp ', ac0brp(sr) C if (ac0crp(sr).lt.0.0.or.ac0crp(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0crp ', ac0crp(sr) C if (ac0drp(sr).lt.0.0.or.ac0drp(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0drp ', ac0drp(sr) C if (ac0aht(sr).lt.0.0.or.ac0aht(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0aht ', ac0aht(sr) C if (ac0bht(sr).lt.0.0.or.ac0bht(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0bht ', ac0bht(sr) C if (ac0ssa(sr).lt.0.0.or.ac0ssa(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0ssa ', ac0ssa(sr) C if (ac0ssb(sr).lt.0.0.or.ac0ssb(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0ssb ', ac0ssb(sr) C if (ac0sla(sr).lt.0.0.or.ac0sla(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0sla ', ac0sla(sr) C if (ac0hue(sr).lt.0.0.or.ac0hue(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0hue ', ac0hue(sr) C if (ac0lfe(sr).lt.0.0.or.ac0lfe(sr).gt.tstmax) * write(*,*) 'day ',day,' ac0lfe ', ac0lfe(sr) C if (dmpflg) write(*,*) 'end dbgdmp' C end