WETA_BB90132 21 1 631 412 68 440 384 488 317 510 147 532 132 595 157 650 692 715 285 880 258 0fdb WETA_BB90132 21 1 632 412 72 440 379 488 355 510 141 532 129 595 158 650 643 715 283 880 256 0fd6 WETA_BB90132 21 1 633 412 72 440 377 488 334 510 143 532 126 595 161 650 675 715 284 880 263 0fcf WETA_BB90132 21 1 634 412 73 440 382 488 324 510 145 532 126 595 161 650 678 715 285 880 258 0fd6 WETA_BB90132 21 1 635 412 68 440 379 488 334 510 147 532 127 595 167 650 675 715 286 880 256 0fe7 WETA_BB90132 21 1 636 412 74 440 383 488 321 510 146 532 126 595 169 650 796 715 288 880 258 0fe4 WETA_BB90132 21 1 637 412 71 440 385 488 345 510 157 532 127 595 170 650 727 715 283 880 259 0fdb WETA_BB90132 21 1 638 412 68 440 386 488 373 510 153 532 127 595 173 650 841 715 287 880 257 0fe2 WETA_BB90132 21 1 639 412 71 440 389 488 311 510 153 532 129 595 167 650 732 715 285 880 257 0fda WETA_BB90132 21 1 640 412 73 440 391 488 331 510 147 532 125 595 174 650 668 715 283 880 259 0fd4 WETA_BB90132 21 1 641 412 72 440 395 488 346 510 151 532 130 595 175 650 677 715 286 880 259 0fd9 WETA_BB90132 21 1 642 412 72 440 395 488 338 510 150 532 138 595 166 650 663 715 281 880 255 0fd4 WETA_BB90132 21 1 643 412 71 440 396 488 327 510 156 532 145 595 166 650 690 715 287 880 258 0fe0 WETA_BB90132 21 1 644 412 78 440 401 488 339 510 156 532 147 595 164 650 782 715 285 880 257 0fdd WETA_BB90132 21 1 645 412 66 440 399 488 322 510 155 532 153 595 165 650 685 715 286 880 257 0fe3 WETA_BB90132 21 1 646 412 72 440 404 488 329 510 153 532 156 595 175 650 737 715 285 880 257 0fda WETA_BB90132 21 1 647 412 71 440 406 488 355 510 151 532 145 595 167 650 720 715 282 880 258 0fce WETA_BB90132 21 1 648 412 77 440 405 488 368 510 154 532 144 595 164 650 690 715 285 880 254 0fdc WETA_BB90132 21 1 649 412 66 440 412 488 326 510 155 532 151 595 166 650 679 715 283 880 259 0fde WETA_BB90132 21 1 650 412 71 440 409 488 328 510 156 532 152 595 166 650 685 715 282 880 259 0fd8 WETA_BB90132 21 1 651 412 72 440 410 488 325 510 154 532 148 595 166 650 722 715 287 880 263 0fca WETA_BB90132 21 1 652 412 70 440 413 488 322 510 150 532 148 595 163 650 673 715 284 880 258 0fc8 WETA_BB90132 21 1 653 412 74 440 413 488 329 510 155 532 151 595 167 650 678 715 285 880 262 0fd8 WETA_BB90132 21 1 654 412 73 440 417 488 334 510 161 532 149 595 173 650 716 715 285 880 259 0fd8 WETA_BB90132 21 1 655 412 71 440 418 488 321 510 156 532 147 595 176 650 963 715 284 880 253 0fd6 WETA_BB90132 21 1 656 412 66 440 419 488 320 510 163 532 146 595 172 650 687 715 285 880 257 0fdc WETA_BB90132 21 1 657 412 73 440 424 488 322 510 148 532 148 595 164 650 684 715 283 880 256 0fd9 WETA_BB90132 21 1 658 412 67 440 422 488 324 510 165 532 150 595 160 650 680 715 284 880 255 0fcd WETA_BB90132 21 1 659 412 71 440 428 488 320 510 152 532 147 595 161 650 681 715 284 880 257 0fd1 WETA_BB90132 21 1 660 412 66 440 424 488 322 510 159 532 148 595 159 650 686 715 285 880 256 0fdf WETA_BB90132 21 1 661 412 73 440 426 488 321 510 146 532 147 595 152 650 685 715 286 880 258 0fd5 WETA_BB90132 21 1 662 412 66 440 432 488 324 510 156 532 150 595 145 650 687 715 286 880 258 0fd7 WETA_BB90132 21 1 663 412 76 440 427 488 324 510 149 532 153 595 148 650 677 715 283 880 257 0fe0 WETA_BB90132 21 1 664 412 73 440 429 488 324 510 151 532 157 595 155 650 691 715 284 880 256 0fd7 WETA_BB90132 21 1 665 412 76 440 433 488 336 510 149 532 157 595 162 650 700 715 284 880 264 0fd4 WETA_BB90132 21 1 666 412 69 440 438 488 328 510 150 532 144 595 164 650 703 715 281 880 266 0fd5 WETA_BB90132 21 1 667 412 74 440 462 488 324 510 157 532 140 595 161 650 1406 715 359 880 362 0fff WETA_BB90132 21 1 668 412 71 440 478 488 366 510 152 532 143 595 162 650 782 715 315 880 273 0fd8 WETA_BB90132 21 1 669 412 77 440 456 488 240 510 155 532 137 595 151 650 681 715 284 880 244 0fd7 WETA_BB90132 21 1 670 412 73 440 459 488 224 510 154 532 130 595 144 650 676 715 291 880 253 0fcc WETA_BB90132 21 1 671 412 70 440 443 488 225 510 151 532 152 595 162 650 627 715 285 880 256 0fc7 WETA_BB90132 21 1 672 412 70 440 443 488 253 510 149 532 149 595 167 650 633 715 282 880 254 0fd3 WETA_BB90132 21 1 673 412 80 440 445 488 214 510 143 532 161 595 169 650 670 715 287 880 257 0fd3 WETA_BB90132 21 1 674 412 74 440 446 488 215 510 151 532 161 595 169 650 625 715 280 880 253 0fcd WETA_BB90132 21 1 675 412 75 440 444 488 212 510 154 532 156 595 164 650 653 715 288 880 257 0fd9 WETA_BB90132 21 1 676 412 70 440 450 488 216 510 155 532 161 595 165 650 761 715 282 880 255 0fcc WETA_BB90132 21 1 677 412 71 440 448 488 212 510 155 532 156 595 169 650 644 715 287 880 258 0fe1 WETA_BB90132 21 1 678 412 68 440 451 488 208 510 157 532 151 595 168 650 633 715 286 880 257 0fdf WETA_BB90132 21 1 679 412 72 440 455 488 216 510 150 532 160 595 181 650 623 715 281 880 255 0fca WETA_BB90132 21 1 680 412 76 440 455 488 211 510 152 532 156 595 181 650 924 715 299 880 268 0fd9 WETA_BB90132 21 1 681 412 67 440 461 488 219 510 173 532 153 595 171 650 1007 715 278 880 274 1001 WETA_BB90132 21 1 682 412 69 440 479 488 247 510 170 532 154 595 163 650 983 715 310 880 321 0fd5 WETA_BB90132 21 1 683 412 72 440 471 488 222 510 155 532 164 595 178 650 937 715 295 880 234 0fd9 WETA_BB90132 21 1 684 412 77 440 477 488 251 510 154 532 163 595 167 650 735 715 298 880 256 0fe6 WETA_BB90132 21 1 685 412 70 440 465 488 209 510 159 532 159 595 170 650 604 715 288 880 251 0fd9 WETA_BB90132 21 1 686 412 71 440 471 488 210 510 159 532 152 595 166 650 766 715 322 880 251 0fcc WETA_BB90132 21 1 687 412 66 440 467 488 213 510 162 532 153 595 167 650 807 715 313 880 246 0fd5 WETA_BB90132 21 1 688 412 68 440 469 488 227 510 149 532 154 595 170 650 1082 715 300 880 238 1008 WETA_BB90132 21 1 689 412 73 440 473 488 209 510 144 532 160 595 177 650 725 715 297 880 493 0fe4 WETA_BB90132 21 1 690 412 67 440 489 488 235 510 152 532 157 595 174 650 902 715 300 880 431 0fcd WETA_BB90132 21 1 691 412 73 440 470 488 222 510 157 532 180 595 171 650 1225 715 362 880 421 0ff1 WETA_BB90132 21 1 692 412 74 440 507 488 260 510 179 532 156 595 178 650 1326 715 329 880 492 1011 WETA_BB90132 21 1 693 412 71 440 475 488 224 510 163 532 156 595 176 650 1375 715 326 880 255 1008 WETA_BB90132 21 1 694 412 77 440 520 488 253 510 158 532 167 595 185 650 1356 715 354 880 493 1012 WETA_BB90132 21 1 695 412 67 440 496 488 213 510 176 532 160 595 172 650 1135 715 273 880 260 1002 WETA_BB90132 21 1 696 412 72 440 480 488 189 510 161 532 163 595 178 650 879 715 296 880 319 0fef vgTstep, avgZ, avgUb, avgVb & , avgU, avgV, avgT, avgR & , avgO, avgW, avgVisc, avgAkv, avgAkt, avgAks & , avgHbl, avgHbbl & , avgHel & , avgBostr & , wrthis & , wrtavg character*128 date_str, title character*128 ininame, grdname, hisname & , rstname, frcname, bulkname, usrname & , avgname & , clmname character*52 vname(4, 150) common /cncscrum/ date_str, title & , ininame, grdname, hisname & , rstname, frcname, bulkname, usrname & , avgname & , clmname & , vname real cff integer i,ierr, lstr,lvar,lenstr, nf_fread, advance_cycle if (may_day_flag.ne.0) return if (itsss.eq.0 .or. iic.eq.0) then lstr=lenstr(frcname) if (ncidfrc.eq.-1) then ierr=nf_open(frcname(1:lstr), nf_nowrite, ncidfrc) if (ierr. ne. nf_noerr) goto 4 endif ierr=nf_inq_varid (ncidfrc, 'sss_time', sss_tid) if (ierr .ne. nf_noerr) then write(stdout,3) 'sss_time', frcname(1:lstr) goto 99 endif lvar=lenstr(vname(1,indxSSS)) ierr=nf_inq_varid (ncidfrc, vname(1,indxSSS)(1:lvar), sss_id) if (ierr .eq. nf_noerr) then ierr=nf_inq_varndims (ncidfrc, sss_id, i) if (ierr. eq. nf_noerr) then if (i.gt.1) then lsssgrd=1 else lsssgrd=0 endif endif endif if (ierr .ne. nf_noerr) then write(stdout,3) vname(1,indxSSS)(1:lvar), frcname(1:lstr) goto 99 endif call set_cycle (ncidfrc, sss_tid, ntsss, & sss_cycle, sss_ncycle, sss_rec) if (may_day_flag.ne.0) return itsss=2 sss_time(1)=-1.E+20 sss_time(2)=-1.E+20 endif 1 i=3-itsss cff=time+0.5*dt if (sss_time(i).le.cff .and. cff.lt.sss_time(itsss)) return ierr=advance_cycle (sss_cycle, ntsss, sss_ncycle, sss_rec) if (ierr .ne. 0) goto 7 ierr=nf_get_var1_double (ncidfrc, sss_tid, sss_rec, cff) if (ierr .ne. nf_noerr) then write(stdout,6) 'sss_time', sss_rec goto 99 endif sss_time(i)=cff*day2sec+sss_cycle*sss_ncycle if (sss_time(itsss).eq.-1.E+20) sss_time(itsss)=sss_time(i) if (lsssgrd.eq.1) then ierr=nf_fread (sssg(-1,-1,i), ncidfrc, sss_id, & sss_rec, r2dvar) else ierr=nf_get_var1_double (ncidfrc, sss_id, sss_rec, sssp(i)) endif if (ierr .ne. nf_noerr) then write(stdout,6) 'SSS', sss_rec goto 99 endif itsss=i write(stdout,'(6x,A,16x,A,1x,g12.4,1x,I4)') & 'GET_SSS -- Read SSS fields', & 'for time =', cff & , mynode if (ntsss.gt.1) goto 1 if (ntsss.eq.1) return 3 format(/,' GET_SSS - ERROR: unable to find forcing variable', & ': ',a,/,11x,'in forcing NetCDF file: ',a) 4 write(stdout,5) frcname(1:lstr) 5 format(/,' GET_SSS - ERROR: unable to open forcing NetCDF ', & 'file: ',a) goto 99 6 format(/,' GET_SSS - ERROR while reading variable: ',a,2x, & ' at TIME index = ',i4) 7 write(stdout,8) sss_rec, ntsss, frcname(1:lstr), tdays, & sss_time(itsss)*sec2day 8 format(/,' GET_SSS - ERROR: requested time record ',I4, & 1x,'exeeds the last available', /, 11x,'record ',I4, & 1x,'in forcing NetCDF file: ', a, /, 11x,'TDAYS = ', & g12.4,2x,'last available SSS_TIME = ',g12.4) 99 may_day_flag=2 return end subroutine set_sss_tile (Istr,Iend,Jstr,Jend) implicit none integer Istr,Iend,Jstr,Jend, i,j, it1,it2 real cff, cff1,cff2, cff3,cff4, val1 integer LLm,Lm,MMm,Mm,N, LLm0,MMm0 parameter (LLm0=336, MMm0=544, N=40) integer Lmmpi,Mmmpi,iminmpi,imaxmpi,jminmpi,jmaxmpi common /comm_setup_mpi/ Lmmpi,Mmmpi, & iminmpi,imaxmpi,jminmpi,jmaxmpi parameter (LLm=LLm0, MMm=MMm0) integer NSUB_X, NSUB_E, NPP integer NP_XI, NP_ETA, NNODES parameter (NP_XI=6, NP_ETA=8, NNODES=NP_XI*NP_ETA) parameter (NSUB_X=1, NSUB_E=1, NPP=1) integer NT, itemp, & ntrc_salt, ntrc_pas, ntrc_bio, ntrc_sed & , ntrc_diats, ntrc_diauv, ntrc_diabio & ,itrc_bio & , isalt & , iNO3, iSIO4, iNH4, iS1, iS2,iZZ1, iZZ2, iDD, iDDSi & , iPO4, iOX, iTCO2, iTALK parameter (itemp=1) parameter (ntrc_salt=1) parameter (isalt=itemp+1) parameter (ntrc_pas=0) parameter (ntrc_bio=13,itrc_bio=itemp+ntrc_salt+ntrc_pas+1) parameter (iNO3=itrc_bio, iSIO4=iNO3+1, & iNH4=iNO3+2, iS1=iNO3+3, & iS2=iNO3+4, iZZ1=iNO3+5, & iZZ2=iNO3+6, iDD=iNO3+7, & iDDSi=iNO3+8, & iPO4=iNO3+9, iOX=iNO3+10, & iTCO2=iNO3+11, iTALK=iNO3+12) parameter (ntrc_sed=0) parameter (NT=itemp+ntrc_salt+ntrc_pas+ntrc_bio+ntrc_sed) parameter (ntrc_diats=0, ntrc_diabio=0) parameter (ntrc_diauv=0) integer Ntides parameter (Ntides=8) integer stdout, Np, padd_X,padd_E parameter (stdout=6, Np=N+1) parameter (Lm=(LLm+NP_XI-1)/NP_XI, Mm=(MMm+NP_ETA-1)/NP_ETA) parameter (padd_X=(Lm+2)/2-(Lm+1)/2) parameter (padd_E=(Mm+2)/2-(Mm+1)/2) integer NSA, N2d,N3d, size_XI,size_ETA integer se,sse, sz,ssz parameter (NSA=28) parameter (size_XI=7+(Lm+NSUB_X-1)/NSUB_X) parameter (size_ETA=7+(Mm+NSUB_E-1)/NSUB_E) parameter (sse=size_ETA/Np, ssz=Np/size_ETA) parameter (se=sse/(sse+ssz), sz=1-se) parameter (N2d=size_XI*(se*size_ETA+sz*Np)) parameter (N3d=size_XI*size_ETA*Np) integer NWEIGHT parameter (NWEIGHT=137) real sustr(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real svstr(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /forces_sustr/sustr /forces_svstr/svstr real sustrg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real svstrg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /smsdat_sustrg/sustrg /smsdat_svstrg/svstrg real sustrp(2), svstrp(2), sms_time(2) real sms_cycle, sms_scale integer itsms, sms_ncycle, sms_rec, lsusgrd, & lsvsgrd,sms_tid, susid, svsid common /smsdat1/ sustrp, svstrp, sms_time common /smsdat2/ sms_cycle, sms_scale common /smsdat3/ itsms, sms_ncycle, sms_rec, lsusgrd, & lsvsgrd,sms_tid, susid, svsid real bustr(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real bvstr(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /forces_bustr/bustr /forces_bvstr/bvstr real bustrg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real bvstrg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /bmsdat_bustrg/bustrg /bmsdat_bvstrg/bvstrg real bms_tintrp(2), bustrp(2), bvstrp(2), tbms(2) real bmsclen, bms_tstart, bms_tend, tsbms, sclbms integer itbms, bmstid,busid, bvsid, tbmsindx logical bmscycle, bms_onerec, lbusgrd, lbvsgrd common /bmsdat/ & bms_tintrp, bustrp, bvstrp, tbms, & bmsclen, bms_tstart, bms_tend, tsbms, sclbms, & itbms, bmstid,busid, bvsid, tbmsindx, & bmscycle, bms_onerec, lbusgrd, lbvsgrd real stflx(-1:Lm+2+padd_X,-1:Mm+2+padd_E,NT) common /forces_stflx/stflx real stflxg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2,NT) common /stfdat_stflxg/stflxg real stflxp(2,NT), stf_time(2,NT) real stf_cycle(NT), stf_scale(NT) integer itstf(NT), stf_ncycle(NT), stf_rec(NT), & lstfgrd(NT), stf_tid(NT), stf_id(NT) common /stfdat/ stflxp, stf_time, stf_cycle, stf_scale, & itstf, stf_ncycle, stf_rec, lstfgrd, & stf_tid, stf_id real btflx(-1:Lm+2+padd_X,-1:Mm+2+padd_E,NT) common /forces_btflx/btflx real dqdt(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real sst(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /forces_dqdt/dqdt /forces_sst/sst real dqdtg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real sstg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /sstdat_dqdtg/dqdtg /sstdat_sstg/sstg real sstp(2), dqdtp(2), sst_time(2) real sst_cycle, scldqdt integer itsst, sst_ncycle, sst_rec, sst_tid, sst_id, & dqdt_id, lsstgrd, sstunused common /sstdat1/ sstp, dqdtp, sst_time common /sstdat2/ sst_cycle, scldqdt common /sstdat3/ itsst, sst_ncycle, sst_rec, sst_tid, sst_id, & dqdt_id, lsstgrd, sstunused real sss(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /forces_sss/sss real sssg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /sssdat_sssg/sssg real sssp(2), sss_time(2) real sss_cycle integer itsss, sss_ncycle, sss_rec, sss_tid, sss_id, & lsssgrd, sssunused common /sssdat/ & sssp, sss_time, sss_cycle, & itsss, sss_ncycle, sss_rec, sss_tid, sss_id, & lsssgrd, sssunused real tair(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real rhum(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real prate(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real radlw(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real radsw(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real uwnd(-1:Lm+2+padd_X,-1:Mm+2+padd_E) real vwnd(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /bulk_tair/tair & /bulk_rhum/rhum & /bulk_prate/prate & /bulk_radlw/radlw & /bulk_radsw/radsw & /bulk_uwnd/uwnd & /bulk_vwnd/vwnd real tairg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real rhumg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real prateg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real radlwg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real radswg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real uwndg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) real vwndg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /bulkdat_tairg/tairg & /bulkdat_rhumg/rhumg & /bulkdat_prateg/prateg & /bulkdat_radlwg/radlwg & /bulkdat_radswg/radswg & /bulk_uwndg/uwndg & /bulk_vwndg/vwndg real tairp(2),rhump(2),pratep(2),radlwp(2),radswp(2) & ,uwndp(2),vwndp(2) real bulk_time(2), bulk_cycle integer tair_id,rhum_id,prate_id,radlw_id,radsw_id, & ltairgrd,lrhumgrd,lprategrd,lradlwgrd,lradswgrd & ,uwnd_id,vwnd_id,luwndgrd,lvwndgrd integer itbulk,bulk_ncycle,bulk_rec,bulk_tid, & bulkunused common /bulkdat1/ & tair_id,rhum_id,prate_id,radlw_id,radsw_id, & ltairgrd,lrhumgrd,lprategrd,lradlwgrd,lradswgrd, & itbulk, bulk_ncycle, bulk_rec, bulk_tid, & bulkunused & ,uwnd_id,vwnd_id,luwndgrd,lvwndgrd common /bulkdat2/ & tairp,rhump,pratep,radlwp,radswp, & bulk_time, bulk_cycle & ,uwndp,vwndp real srflx(-1:Lm+2+padd_X,-1:Mm+2+padd_E) common /forces_srflx/srflx real srflxg(-1:Lm+2+padd_X,-1:Mm+2+padd_E,2) common /srfdat_srflxg/srflxg real srflxp(2),srf_time(2) real srf_cycle, srf_scale integer itsrf, srf_ncycle, srf_rec, & lsrfgrd, srf_tid, srf_id common /srfdat/ & srflxp,srf_time, srf_cycle, srf_scale, & itsrf, srf_ncycle, srf_rec, lsrfgrd, srf_tid, srf_id real dt, dtfast, time, time_start, tdays integer iic, kstp, krhs, knew, next_kstp & , iif, nstp, nrhs, nnew logical PREDICTOR_2D_STEP common /time_indices/ dt,dtfast, time,time_start, tdays, & iic, kstp, krhs, knew, next_kstp, & iif, nstp, nrhs, nnew, & PREDICTOR_2D_STEP real time_avg, rho0 & , rdrg, rdrg2, Cdb_min, Cdb_max, Zob & , xl, el, visc2, visc4, gamma2 real theta_s, theta_b, Tcline, hc real sc_w(0:N), Cs_w(0:N), sc_r(N), Cs_r(N) real rx0, rx1 real tnu2(NT),tnu4(NT) real weight(2,0:NWEIGHT) real x_sponge, v_sponge real tauT_in, tauT_out, tauM_in, tauM_out integer numthreads, ntstart, ntimes, ninfo & , ndtfast,nfast, nrrec, nrst, nwrt & , ntsavg, navg logical ldefhis logical got_tini(NT) common /scalars_main/ & time_avg, rho0, rdrg, rdrg2 & , Zob, Cdb_min, Cdb_max & , xl, el, visc2, visc4, gamma2 & , theta_s, theta_b, Tcline, hc & , sc_w, Cs_w, sc_r, Cs_r & , rx0, rx1, tnu2, tnu4 & , weight & , x_sponge, v_sponge & , tauT_in, tauT_out, tauM_in, tauM_out & , numthreads, ntstart, ntimes, ninfo & , ndtfast,nfast, nrrec, nrst, nwrt & , ntsavg, navg & , ldefhis & , got_tini logical synchro_flag common /sync_flag/ synchro_flag integer may_day_flag integer tile_count, first_time, bc_count & , bio_count common /communicators_i/ & may_day_flag, tile_count, first_time, bc_count & , bio_count real hmin, hmax, grdmin, grdmax, Cu_min, Cu_max common /communicators_r/ & hmin, hmax, grdmin, grdmax, Cu_min, Cu_max real*16 volume, avgke, avgpe, avgkp, bc_crss & , global_sum(0:2*NT+1) common /communicators_rq/ & volume, avgke, avgpe, avgkp, bc_crss & , global_sum real*4 CPU_time(0:31,0:NPP) integer proc(0:31,0:NPP),trd_count common /timers/CPU_time,proc,trd_count logical EAST_INTER, WEST_INTER, NORTH_INTER, SOUTH_INTER integer mynode, ii,jj, p_W,p_E,p_S,p_N, p_SW,p_SE, p_NW,p_NE common /comm_setup/ mynode, ii,jj, p_W,p_E,p_S,p_N, p_SW,p_SE, & p_NW,p_NE, EAST_INTER, WEST_INTER, NORTH_INTER, SOUTH_INTER real pi, deg2rad, rad2deg parameter (pi=3.14159265358979323846, deg2rad=pi/180., & rad2deg=180./pi) real Eradius, g, day2sec,sec2day, jul_off, & year2day,day2year parameter (Eradius=6371315.0, day2sec=86400., & sec2day=1./86400., jul_off=2440000., & year2day=365.25, day2year=1./365.25) parameter (g=9.81) real Cp parameter (Cp=3985.0) real vonKar parameter (vonKar=0.41) integer IstrR,IendR,JstrR,JendR if (Istr.eq.1) then if (WEST_INTER) then IstrR=Istr-2 else IstrR=Istr-1 endif else IstrR=Istr endif if (Iend.eq.Lmmpi) then if (EAST_INTER) then IendR=Iend+2 else IendR=Iend+1 endif else IendR=Iend endif if (Jstr.eq.1) then if (SOUTH_INTER) then JstrR=Jstr-2 else JstrR=Jstr-1 endif else JstrR=Jstr endif if (Jend.eq.Mmmpi) then if (NORTH_INTER) then JendR=Jend+2 else JendR=Jend+1 endif else JendR=Jend endif it1=3-itsss it2=itsss cff=time+0.5*dt cff1=sss_time(it2)-cff cff2=cff-sss_time(it1) if (sss_cycle.lt.0.) then if (iic.eq.ntstart) then if (lsssgrd.eq.1) then do j=JstrR,JendR do i=IstrR,IendR sss(i,j)=sssg(i,j,itsss) enddo enddo else val1=sssp(itsss) do j=JstrR,JendR do i=IstrR,IendR sss(i,j)=val1 enddo enddo endif endif elseif (cff1.ge.0. .and. cff2.ge.0.) then if (Istr+Jstr.eq.2 .and. cff1.lt.dt) synchro_flag=.TRUE. cff=scldqdt/(cff1+cff2) cff3=cff1*cff cff4=cff2*cff cff=1./(cff1+cff2) cff1=cff1*cff cff2=cff2*cff if (lsssgrd.eq.1) then do j=JstrR,JendR do i=IstrR,IendR sss(i,j)=cff1*sssg(i,j,it1)+cff2*sssg(i,j,it2) enddo enddo else val1=cff1*sssp(it1)+cff2*sssp(it2) do j=JstrR,JendR do i=IstrR,IendR sss(i,j)=val1 enddo enddo endif else if (Istr+Jstr.eq.2) then write(stdout,1) 'sss_time', tdays, sss_time(it2)*sec2day 1 format(/,' SET_SSS - current model time exceeds ending', & 1x,'value for variable: ',a,/,11x,'TDAYS = ',g12.4, & 2x,'TEND = ',g12.4) may_day_flag=2 endif endif return end ELF>P@@  "'83 I0_@uP% ;(Q0g8}@HPX`hpx-CYo 5Kaw (`p'0=Si@@@@@@ @@ @6@K @`@ @v @ @@@@@@@@H&@;@Q @j`@@@ @`@ @!@ #@'`$@>%@U&@l (@`)@*@+@ -@`.@/@ ) @? PU Ht h x p     7 N i        ! 8 O f } (   0   @ P H5 XL `c 8~ p h x     6 M d {      3Jax0(8P @H0hGX^up`x-D[r*A8XHo@X(0` 'P>hUplx >Ul2I` w(HPx8@0X/hF]pt`3Jax( 0G@^uH08PX`p-hD[xr*AXo0'>(U l8@PHXh` $p;xRi(?V5 E O  Z  - (8 (C Q (_ (g @p hIy @ @q P &; &" 'pk 0(p  h 8e z $5`-38<=>@ Gm} ( ( ( ( ( ( ( ( ( ( ( ( ( (! (/ (< (J (X (f (v ( ( ( ( ( ( ( ( ( :A :A hI @ @  @+ @7 @B @MXi0z P < ( ( P P0  ( ( P  P/ 7 X'E NT 8\ X'j (u ( ( ( ( P P P P P P P  L    P-  HHLIHHc D0OE^2LcMiI;H"LHIHLd$ fEE(HG$ID$D)$E,HC($ID$E3A(,HC(4ID$X)$A<HGID$X$PE( HCID)$)$AHC I)$$`AHG<I$h$pE,HG$I$xD$EHGID$D$A$HG I$D$EHC<I$D$A4HC,I$$xA$HCI$$pAHCI$$hA HCI$$E<HC ID$0$AHG<ID$@$8A H$(C I$`AH$HG<ID$AH$G<ID$ AH$G<ID$A$G<ID$D$D$CI$D)t$pLl$(G<ID$Lt$0L|$8G<ID$H\$@Hl$HG(<ID)<$D$D$G(<D)|$$XHc.HcH;HHIIH$H$H$L$LT$PHL$XH|$`Ht$hH$f.E?H=LH$HIIHcwHIHi@HHHOH$HIIL$PK IO, L$H3H$HO<0L$@H$ML$8I($H$@H$HML$0IL$0ML$(IL$(ML$ IL$ ML$IML$IML$IML$IL$ML$IML$L$L$8L$PCL4IHHfB28J20p~+H ,fB18J100lfXDmf_fC08fC00ff_fB38fB30.fعf_fB08fB00pi0f0f_fC18fC10D0qfDX/rfD_fG28fG20D޲fDjfD_fG78fG70DfDؤfD_fG38fG30DpT5fD5fD_H$fF58fF50D0vfDXwfD_H$fF58fF50DɷfDUfD_H$fF58fF50f؏f_fC48fC40p?:f:f_ؐH$fB58fB500z{fX|HPf_fC58H(fC50cL-IHIHE3H$HKl5MD=HL OdLM\5H5H$H>5M,9A5I:fHfAL5LB)TAfOtO|fALLB)dhAhf@}OtO| fA@}LhL@}B)t ADfAfLHffDfC)C)\C)l Ot(O|0C)|0Ol8KT@F)D0II(L$0H5IIL$A(A I(H^EBIL$MA6A~(F.IE \AYAYDY^UEaAAhI^]D^EBjA[>HIE\\EY^mD^eBlIFD>fEnfHnIIHFd.fvfBo,fBoLfBo6N,>H$L>H$D(D$D( $H$H$L$0fApfofvDPB(TAt{G(tfDofD^fDbfEoETGUlfEoEVG)tF(lF(tfEYfEYfD^fD^GU|ETETEVGUdEVG)lG)tfDoffDffDofDfEEPAG$G(lfGdfE\fD^fEbfEoETGUdEVG)lF,F(dfFlfE\fEYfD^fEoETGUlEVG)dF$F(lfFdfE\fEYfD^ETGUtEVG)lfDfDofAfDvfADPAG$fGdfEfE\fD^fbfDoDTGUlfDoFd$@YQD$H$\$d$(l$ht$8|$0L H$HIHIHA H^EC$DYYEXl$xDXA<9C4\YA1XYY|$PD\D$pDXA(L HHI$HIH@XHA^CYA\ A,9C$\Yd$xA 1XYL$pYl$PXXL HHHIIHA, fYl$XA9YT$PA\,1C\Y\$HC$XYd$hHX\AD$pH=HHHIfYT$P,7\Y\$XB$XYd$hHX@\D$pHHHHIHIIBIHYD$`Ht$8Id$PD(B D8HYX0FY (DYYDX FXDYHBB\Y,0\Yl$hY\$HDDXD\<8YD\EXDYL$XEXAXLH\$XHHHIHA 2IYI@ AYA\ :C$ YL$PAX YXCX$CYT$HYd$8\XLH\$XHHHIHA 2IYI@8AYA\ :C$ YL$PAX YXCX$CYT$HYd$8\XL HT$XHIHHHA 1IY@@ACYA\ 9YYL$PA\ CXXY\$8XHHH YL$P@HD$ \LHHHl$XHIHA<IA IA3IYXYXA\;A C$YT$HYY\$@\CX$ XYd$C4XYD$`CX4+Yt$(X(LHIHd$XHIIA HC<YL$PYHA2l$H@PAC4\YC\< YXXA:YY|$8\XHIHIDT$(Hl$@A(|$XIB4HBHYIYY 0B\4DD$HIB$(HDHAYDYXHF$\YDXd$`DXBD 0AYEYD\E\DXEY8EXYT$XYAXHHH YL$0^@(D$ XH HHHYT$X 1YL$(\XT$`XX(LHHId$0HIIAA;fC,YYYA^ C^C^,+H$HH5I\ @0HD$pAYhH$fD$D$p\H|=H|$D$ Ht$ HD$(HT$(HHH(YHQHXYVb(XY^L$oXL$A@`\$0d$8l$@d$8l$@D(|$(D^Yf\t$0YY\\$A(^H$D\^*,$^D^H$H$H$AU!DHĐA]Ht&ATAUAVH`IIH4$IHT$@LHt$HT$@H$fAAMf/$v\$PL$HHD$\$HL$PH$HT$(H HIX\H|$YHt$BT HT$@T$8D$L$ L$(D$0D$0AL$(T$t$ ((|$@f\d$H\\$PYY\\Yf/s2H H(HYXHHTTf/v-T$P(t$H(\Yt$8Xf.u$z"((^(\f.zHH(A$H Tf/wwD$H|$t$ Ht$L$(HT$@D$0T$fD$0L$(t$ HD$f/v D$PH$ D$HHD$AAdH`A^A]A\ÐHt&H|'UHHHH IIIIDQH ?LAi8L AD((D(DO(D(A^DMD(D^AHAb@IEXIDYXYXEYD^C4DXYEYEYDYEXuEXAYAKHIDpHYX(MAYXDXDXmEXC,E(YMEYAXmAk(DUD(DMDXYEYA^DHDmIEYA\s0DmLDED(DUD(Ei0D^@ D^QPDmDXD^UDXDXAR Ye(DxXYA(A^DaXEYD^H A(^D_XE^A(A^EpXEYD^YPA(A^E\EY^E(^U\E^(\YME^E^DpXDYuE^(^E^XxY^U\A^(XUYYUD\DYDY}}D^mE\^6E\D]E\D^DeEXEYE\A^DuD^E^u(YAYYE\DYEDXY\mYDY]^^D^DXD\E\D\D:H]HHATAUAVAWSUVMMIHIHHHHHL ^QY HY Y  $AYE $HD$@A\YA^Y][A_A^A]A\LHHII(YI4 IYKO  \YYA)\XÐHt&H|'UHHHHH \X^mHH(YH$(IYHHH8(FD(D,IIHYX$0DYYDYDYX HBHF$IDYDXEYDYXFX$EXDXDYDYD(EY $EXEXEXA(HH $^YYH]H&ATAUAVSHxHIT$8MQMBHL$@IT$0XD$(LHHIL$8HHHA0^\$(A$YA\8CYAX$YXD$ XfA,\Y,$H5HHHL$8H,$f^T$(YX>D$XYX\Y(HHt$0IHIDD$8ID0HDYHB HAYDX8BHYD\BIEY\d$(DXD^D$HYFX<0YXHD 8EYYD,D\AYE\DXDYT$ EXAXAD$*H=HHHIHY^ YL$(^BYT$@7$\YHHHL$Y YT$@L$\Y$HHHL$Y YT$@L$\Y$L$YL$HHD$((L$^($LH$HHHHIA ^A0YAYAX8AXYXCXY\$8XHHH$ ^L$(\HHHD$H ^L$(\EHIIDL$(Il$8HHHF0HB4IE^YIB^48AYYt$0$H \AYY\YD$@XBDA\DY\$ B<A\DYd$pt$`|$XDD$PH5HHHHD$h^T$( YL$ \D$0X>YY\$8\H=HL$HT$Hd$hH3HDI\$,7<Y YDY^Yd$pd$h$(AE\A $IL$0T$ F F'D$(L$@T$HD(D(Dd$`D(D$0D(D^D^HHDXEXE^D(A\DX\$XDd$PE^D\$hA\D^D(A\D^d$pD(AXEYE(EYEXEYEXEXDYE^E(D(Dl$ A\^D$(D^DYD\A\DY(DXEQEXEYA^\T fD/sA(dD$(t$A(((A^^DL$A(E((AYXE^AYX^H$AXX^^DYh^,$E^(@DHHx[A^A]A\HAUHHHI'($$^HI$$HHIN IA H ^YH<0HHXH<0H YNX)IA\YY\X(AYmXHHYHA]ÐHt&H?h㈵>h㈵>??@@Y@Y@??.@.@9@9@''''((((@9@33333ƚ@Zd;O?Q?\(\?@@Tư>333333?@@?Y@?????ffffff?Q?p= ף?p= ף?{Gzt??333333333333??@!@@333333@333333Dl?>@(\?@@Qw@{Gz?RQ?333333?.@?@{Gz?8@@"@@@&S @! ?Hzg_333332@O?ư>fffffq@Y@C3@@@Gz?k ? -ZdSt$V@}bIk@Ϡ?9ѮBO?1*ԙ?tlt?\W@ͪN@[7@'?jIG9?݃Ds?$@fffff@/$O@ 0*?ǘ#@IӠh?S@HPsג?ʕ@mO?C6&a@0?^)ˀb@HV?a4o8&䃞ͪ?(\9@?5^I|SS?3333sV}b?(\@\@n|ys?s2@}"?vZZӂe@R[7ۊe?h˹W?Zd;O;@3?{G djt$2@ WF'>?; ~@Z*oG8E1@ezFwP?q= ףX]@d/$U3@Hzg@@ǘ??ޫV&?1*T(~ @q= ף|-Cb@aTR'?{G 7@{G]@ r@n?@@"~a@+7@@R@Q@ʡE\@QIt@|?5G@:MھW?ؘ@x&1H)@ffffff?IӠh.?Zd;%@Q?X@t>?>6e]:?A@|=Ac]b@ @7@= ףpũ@2U0*S@ZӼ(@$I??@:0yE>@@/tG??Y@/neL F%a@\`@S8?JXƆ"?S'ц>A{GzTF@zRx  & f8    6  (08~>80( dp4|p _ T` C $  Ƅ4 (08@80( , d@  (& ( ,l biology_.f.strtab.symtab.data.bss.textbulk_uwnd_bulk_vwnd_scalars_main_forces_srflx_mask_r_grid_hz_grid_zw_grid_zr_ocean_t_time_indices_biology_tile__GLOBAL_OFFSET_TABLE_.rodata.L_2il0floatpacket.11.L_2il0floatpacket.12.L_2il0floatpacket.13.L_2il0floatpacket.66.L_2il0floatpacket.14.L_2il0floatpacket.15.L_2il0floatpacket.16.L_2il0floatpacket.17.L_2il0floatpacket.18.L_2il0floatpacket.19.L_2il0floatpacket.20.L_2il0floatpacket.21.L_2il0floatpacket.22.L_2il0floatpacket.23.L_2il0floatpacket.24.L_2il0floatpacket.25.L_2il0floatpacket.26.L_2il0floatpacket.27.L_2il0floatpacket.28.L_2il0floatpacket.29.L_2il0floatpacket.30.L_2il0floatpacket.31.L_2il0floatpacket.32.L_2il0floatpacket.33.L_2il0floatpacket.34.L_2il0floatpacket.35.L_2il0floatpacket.36.L_2il0floatpacket.37.L_2il0floatpacket.38.L_2il0floatpacket.39.L_2il0floatpacket.40.L_2il0floatpacket.41.L_2il0floatpacket.42.L_2il0floatpacket.43.L_2il0floatpacket.44.L_2il0floatpacket.45.L_2il0floatpacket.46.L_2il0floatpacket.47.L_2il0floatpacket.48.L_2il0floatpacket.49.L_2il0floatpacket.50.L_2il0floatpacket.51.L_2il0floatpacket.52.L_2il0floatpacket.53.L_2il0floatpacket.54.L_2il0floatpacket.55.L_2il0floatpacket.57.L_2il0floatpacket.58.L_2il0floatpacket.59.L_2il0floatpacket.60.L_2il0floatpacket.61.L_2il0floatpacket.62.L_2il0floatpacket.63.L_2il0floatpacket.65.L_2il0floatpacket.64expo2flux_co2flux_biology_tile_$OXFLX.0.1biology_tile_$CO2FLX.0.1biology_tile_$TEMP.0.1biology_tile_$SALT.0.1biology_tile_$NO3.0.1biology_tile_$SIO4.0.1biology_tile_$NH4.0.1biology_tile_$S1.0.1biology_tile_$S2.0.1biology_tile_$ZZ1.0.1biology_tile_$ZZ2.0.1biology_tile_$DD.0.1biology_tile_$DDSI.0.1biology_tile_$PO4.0.1biology_tile_$OX.0.1biology_tile_$TCO2.0.1biology_tile_$TALK.0.1biology_tile_$ZW.0.1biology_tile_$ZR.0.1biology_tile_$DZ0.0.1biology_tile_$SINKS2.0.1biology_tile_$SINKDD.0.1biology_tile_$SKDDSI.0.1biology_tile_$ADPT.0.1biology_tile_$PAR.0.1biology_tile_$REMVZ2.0.1biology_tile_$SMS1.0.1biology_tile_$SMS3.0.1biology_tile_$SMS4.0.1biology_tile_$SMS5.0.1biology_tile_$SMS6.0.1biology_tile_$SMS7.0.1biology_tile_$SMS8.0.1biology_tile_$SMS2.0.1biology_tile_$SMS9.0.1biology_tile_$SMS10.0.1biology_tile_$SMS11.0.1biology_tile_$SMS12.0.1biology_tile_$SMS13.0.1__NLITPACK_0.0.1scco2_.L_2il0floatpacket.78.L_2il0floatpacket.80.L_2il0floatpacket.79co2calc_.L_2il0floatpacket.83.L_2il0floatpacket.86.L_2il0floatpacket.85.L_2il0floatpacket.87.L_2il0floatpacket.88const_species_.L_2il0floatpacket.104.L_2il0floatpacket.105.L_2il0floatpacket.106.L_2il0floatpacket.206.L_2il0floatpacket.107log.L_2il0floatpacket.110.L_2il0floatpacket.108.L_2il0floatpacket.109.L_2il0floatpacket.111.L_2il0floatpacket.114.L_2il0floatpacket.113.L_2il0floatpacket.117.L_2il0floatpacket.116.L_2il0floatpacket.118.L_2il0floatpacket.112.L_2il0floatpacket.115.L_2il0floatpacket.119.L_2il0floatpacket.123.L_2il0floatpacket.122.L_2il0floatpacket.124.L_2il0floatpacket.121.L_2il0floatpacket.120.L_2il0floatpacket.126.L_2il0floatpacket.128.L_2il0floatpacket.127.L_2il0floatpacket.129.L_2il0floatpacket.130.L_2il0floatpacket.125pow.L_2il0floatpacket.132.L_2il0floatpacket.131.L_2il0floatpacket.133.L_2il0floatpacket.134.L_2il0floatpacket.136.L_2il0floatpacket.135.L_2il0floatpacket.140.L_2il0floatpacket.138.L_2il0floatpacket.143.L_2il0floatpacket.139.L_2il0floatpacket.141.L_2il0floatpacket.144.L_2il0floatpacket.142.L_2il0floatpacket.137.L_2il0floatpacket.146.L_2il0floatpacket.145.L_2il0floatpacket.149.L_2il0floatpacket.148.L_2il0floatpacket.150.L_2il0floatpacket.153.L_2il0floatpacket.147.L_2il0floatpacket.151.L_2il0floatpacket.152.L_2il0floatpacket.156.L_2il0floatpacket.155.L_2il0floatpacket.157.L_2il0floatpacket.160.L_2il0floatpacket.154.L_2il0floatpacket.158.L_2il0floatpacket.159.L_2il0floatpacket.163.L_2il0floatpacket.161.L_2il0floatpacket.166.L_2il0floatpacket.164.L_2il0floatpacket.162.L_2il0floatpacket.165.L_2il0floatpacket.167.L_2il0floatpacket.168.L_2il0floatpacket.169.L_2il0floatpacket.171.L_2il0floatpacket.170.L_2il0floatpacket.174.L_2il0floatpacket.172.L_2il0floatpacket.176.L_2il0floatpacket.173.L_2il0floatpacket.175.L_2il0floatpacket.178.L_2il0floatpacket.181.L_2il0floatpacket.177.L_2il0floatpacket.179.L_2il0floatpacket.183.L_2il0floatpacket.182.L_2il0floatpacket.180.L_2il0floatpacket.189.L_2il0floatpacket.191.L_2il0floatpacket.185.L_2il0floatpacket.190.L_2il0floatpacket.184.L_2il0floatpacket.193.L_2il0floatpacket.187.L_2il0floatpacket.188.L_2il0floatpacket.194.L_2il0floatpacket.186.L_2il0floatpacket.192.L_2il0floatpacket.195.L_2il0floatpacket.196.L_2il0floatpacket.197.L_2il0floatpacket.198.L_2il