/[MITgcm]/MITgcm/pkg/fizhi/slprs.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/slprs.F

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

revision 1.2 by molod, Thu May 18 22:14:23 2006 UTC revision 1.3 by molod, Mon Jun 19 23:37:34 2006 UTC
# Line 3  C $Name$ Line 3  C $Name$
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                      
# Line 63  C*                   COMPUTE SEA LEVEL P Line 63  C*                   COMPUTE SEA LEVEL P
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)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22