/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.18 by jahn, Fri Jul 24 21:47:26 2015 UTC revision 1.20 by jahn, Fri Apr 1 00:02:17 2016 UTC
# Line 599  c sum wavebands for total PAR at the mid Line 599  c sum wavebands for total PAR at the mid
599  #else /* DAR_RADTRANS */  #else /* DAR_RADTRANS */
600  c ------------ FULL RADIATIVE TRANSFER CODE ----------------------------  c ------------ FULL RADIATIVE TRANSFER CODE ----------------------------
601           do ilam = 1,tlam           do ilam = 1,tlam
602            Edwsf(ilam) = oasim_ed(i,j,ilam,bi,bj)            Edwsf(ilam)=oasim_ed(i,j,ilam,bi,bj)*(1 _d 0-fice(i,j,bi,bj))
603            Eswsf(ilam) = oasim_es(i,j,ilam,bi,bj)            Eswsf(ilam)=oasim_es(i,j,ilam,bi,bj)*(1 _d 0-fice(i,j,bi,bj))
604           enddo           enddo
605    
606  #ifdef DAR_RADTRANS_USE_MODEL_CALENDAR  #ifdef DAR_RADTRANS_USE_MODEL_CALENDAR
# Line 610  C cos(solz) is average over daylight per Line 610  C cos(solz) is average over daylight per
610       O                       solz)       O                       solz)
611    
612  #else /* not DAR_RADTRANS_USE_MODEL_CALENDAR */  #else /* not DAR_RADTRANS_USE_MODEL_CALENDAR */
613  C use calendar date for full solar zenith angle computation  C Use calendar date for full solar zenith angle computation.
614  C oj: average light effective at noon?  C Use local noon zenith angle to avoid problems with zero cosine and
615           solz = 0.0 _d 0  C non-zero light.  One should really use a zenith angle compatible with
616           isec = 12*3600        C the light fields, in particular averaged over the same time period.
617             isec = MOD(36.*3600. - 240.*XC(i,j,bi,bj), 86400.)
618           call radtrans_sfcsolz(rad,iyr,imon,iday,isec,           call radtrans_sfcsolz(rad,iyr,imon,iday,isec,
619       I                         XC(i,j,bi,bj),YC(i,j,bi,bj),       I                         XC(i,j,bi,bj),YC(i,j,bi,bj),
620       O                         solz)       O                         solz)
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22