/[MITgcm]/MITgcm/model/src/ini_linear_phisurf.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_linear_phisurf.F

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

revision 1.1 by jmc, Fri Jul 6 21:39:37 2001 UTC revision 1.2 by jmc, Mon Aug 27 18:43:46 2001 UTC
# Line 76  C--- Bo = d/d_p(Phi_surf) = tRef_surf*d/ Line 76  C--- Bo = d/d_p(Phi_surf) = tRef_surf*d/
76           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
77            DO J=1-Oly,sNy+Oly            DO J=1-Oly,sNy+Oly
78             DO I=1-Olx,sNx+Olx             DO I=1-Olx,sNx+Olx
79              IF (Ro_surf(I,J,bi,bj).GT.0. _d 0) THEN              IF ( Ro_surf(I,J,bi,bj).GT.0. _d 0
80         &          .AND. ksurfC(I,J,bi,bj).LE.Nr ) THEN
81                dPIdp = (atm_cp*atm_kappa/atm_po)*                dPIdp = (atm_cp*atm_kappa/atm_po)*
82       &         (Ro_surf(I,J,bi,bj)/atm_po)**(atm_kappa-1. _d 0)       &         (Ro_surf(I,J,bi,bj)/atm_po)**(atm_kappa-1. _d 0)
83                Bo_surf(I,J,bi,bj) = dPIdp*tRef(k_surf(I,J,bi,bj))                Bo_surf(I,J,bi,bj) = dPIdp*tRef(ksurfC(I,J,bi,bj))
84                recip_Bo(I,J,bi,bj) = 1. _d 0 / Bo_surf(I,J,bi,bj)                recip_Bo(I,J,bi,bj) = 1. _d 0 / Bo_surf(I,J,bi,bj)
85              ELSE              ELSE
86                Bo_surf(I,J,bi,bj) = 0.                Bo_surf(I,J,bi,bj) = 0.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22