2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
5 |
#define OLD_PSTAR_SLOPE_TERM |
#undef OLD_PSTAR_SLOPE_TERM |
6 |
|
|
7 |
CBOP |
CBOP |
8 |
C !ROUTINE: CALC_GRAD_PHI_HYD |
C !ROUTINE: CALC_GRAD_PHI_HYD |
103 |
ELSEIF (Ro_surf(i,j,bi,bj).NE.0. _d 0) THEN |
ELSEIF (Ro_surf(i,j,bi,bj).NE.0. _d 0) THEN |
104 |
factPI = (rC(k)/Ro_surf(i,j,bi,bj))**atm_kappa |
factPI = (rC(k)/Ro_surf(i,j,bi,bj))**atm_kappa |
105 |
varLoc(i,j) = phiHydC(i,j) |
varLoc(i,j) = phiHydC(i,j) |
106 |
& *(rStarFacC(i,j,bi,bj)**atm_kappa - factPI) |
& *(pStarFacK(i,j,bi,bj) - factPI) |
107 |
& /(1. _d 0 -factPI) |
& /(1. _d 0 -factPI) |
108 |
& + phi0surf(i,j,bi,bj) |
& + phi0surf(i,j,bi,bj) |
109 |
ENDIF |
ENDIF |