3 |
|
|
4 |
#include "FIZHI_OPTIONS.h" |
#include "FIZHI_OPTIONS.h" |
5 |
subroutine slprs (PHIS,PLE,THZ,lwmask,im,jm,lm,SLP) |
subroutine slprs (PHIS,PLE,THZ,lwmask,im,jm,lm,SLP) |
6 |
C*********************************************************************** |
C********************************************************************** |
7 |
C INPUT |
C INPUT |
8 |
C PHIS .... SURFACE GEOPOTENTIAL (M2/S2) |
C PHIS .... SURFACE GEOPOTENTIAL (M2/S2) |
9 |
C THZ ..... POTENTIAL TEMPERATURE (K) ON Model LEVELS |
C THZ ..... POTENTIAL TEMPERATURE (K) ON Model LEVELS |
63 |
C*********************************************************************** |
C*********************************************************************** |
64 |
|
|
65 |
do i=1,im*jm |
do i=1,im*jm |
66 |
if( lwmask(i).eq.0.0 ) then |
if( lwmask(i).ne.0.0 ) then |
67 |
TM(I) = TM(I) * (PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |
TM(I) = TM(I) * (PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |
68 |
else |
else |
69 |
TM(I) = THZ(I,LM)*(PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |
TM(I) = THZ(I,LM)*(PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |