pkg/gridalt/gridalt_update.F

       subroutine gridalt_update(myThid)
c***********************************************************************
c subroutine gridalt_update
c 
c Purpose: Update the pressure thicknesses of the layers of the 
c          alternative vertical grid (used now for atmospheric physics).
c
c Calculate: dpphys    - new delta r (p*) edge-to-edge of physics grid
c                        using dpphys0 (initial value) and rstarfacC
c
c***********************************************************************
       implicit none
c
#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "fizhi_SIZE.h"
#include "gridalt_mapping.h"
#include "SURFACE.h"

       integer myThid

       integer im1,im2,jm1,jm2
       integer i1,i2,j1,j2,bi,bj
       integer i,j,L

       im1 = 1-OLx
       im2 = sNx+OLx
       jm1 = 1-OLy
       jm2 = sNy+OLy
       i1 = 1
       i2 = sNx
       j1 = 1
       j2 = sNy

#ifdef ALLOW_FIZHI
                                                                                         
       if(usefizhi)then
        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
c
         do L = 1,Nrphys
          do j = j1,j2
          do i = i1,i2
           dpphys(i,j,L,bi,bj) = dpphys0(i,j,L,bi,bj) *
     .                     rstarfacC(i,j,bi,bj)
          enddo
          enddo
         enddo
c
        enddo
        enddo
       endif

#endif

       return
       end
1	molod	1.1	subroutine gridalt_update(myThid)
2			c***********************************************************************
3			c subroutine gridalt_update
4			c
5			c Purpose: Update the pressure thicknesses of the layers of the
6			c alternative vertical grid (used now for atmospheric physics).
7			c
8			c Calculate: dpphys - new delta r (p*) edge-to-edge of physics grid
9			c using dpphys0 (initial value) and rstarfacC
10			c
11			c***********************************************************************
12			implicit none
13			c
14	molod	1.3	#include "PACKAGES_CONFIG.h"
15	molod	1.1	#include "CPP_OPTIONS.h"
16			#include "SIZE.h"
17	molod	1.4	#include "EEPARAMS.h"
18	molod	1.3	#include "PARAMS.h"
19	molod	1.1	#include "GRID.h"
20			#include "fizhi_SIZE.h"
21			#include "gridalt_mapping.h"
22			#include "SURFACE.h"
23
24			integer myThid
25
26			integer im1,im2,jm1,jm2
27			integer i1,i2,j1,j2,bi,bj
28			integer i,j,L
29
30			im1 = 1-OLx
31			im2 = sNx+OLx
32			jm1 = 1-OLy
33			jm2 = sNy+OLy
34			i1 = 1
35			i2 = sNx
36			j1 = 1
37			j2 = sNy
38	molod	1.3
39			#ifdef ALLOW_FIZHI
40	molod	1.1
41	molod	1.3	if(usefizhi)then
42			do bj = myByLo(myThid), myByHi(myThid)
43			do bi = myBxLo(myThid), myBxHi(myThid)
44	molod	1.1	c
45	molod	1.3	do L = 1,Nrphys
46			do j = j1,j2
47			do i = i1,i2
48			dpphys(i,j,L,bi,bj) = dpphys0(i,j,L,bi,bj) *
49	molod	1.2	. rstarfacC(i,j,bi,bj)
50	molod	1.3	enddo
51			enddo
52	molod	1.1	enddo
53	molod	1.3	c
54			enddo
55	molod	1.1	enddo
56	molod	1.3	endif
57
58			#endif
59	molod	1.1
60			return
61			end