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. |