/[MITgcm]/MITgcm/pkg/gridalt/gridalt_update.F
ViewVC logotype

Diff of /MITgcm/pkg/gridalt/gridalt_update.F

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

revision 1.7 by molod, Mon May 23 20:51:09 2005 UTC revision 1.11 by jmc, Sat Jul 7 00:08:09 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "GRIDALT_OPTIONS.h"
5    
6         subroutine gridalt_update(myThid)         subroutine gridalt_update(myThid)
7  c***********************************************************************  c***********************************************************************
8  c subroutine gridalt_update  c subroutine gridalt_update
9  c  c
10  c Purpose: Update the pressure thicknesses of the layers of the  c Purpose: Update the pressure thicknesses of the layers of the
11  c          alternative vertical grid (used now for atmospheric physics).  c          alternative vertical grid (used now for atmospheric physics).
12  c  c
13  c Calculate: dpphys    - new delta r (p*) edge-to-edge of physics grid  c Calculate: dpphys    - new delta r (p*) edge-to-edge of physics grid
# Line 13  c                        using dpphys0 ( Line 15  c                        using dpphys0 (
15  c  c
16  c***********************************************************************  c***********************************************************************
17         implicit none         implicit none
18  c  
 #include "PACKAGES_CONFIG.h"  
 #include "CPP_OPTIONS.h"  
19  #include "SIZE.h"  #include "SIZE.h"
20  #include "EEPARAMS.h"  #include "EEPARAMS.h"
21  #include "PARAMS.h"  #include "PARAMS.h"
# Line 31  c Line 31  c
31         integer i,j,L         integer i,j,L
32         logical  diagnostics_is_on         logical  diagnostics_is_on
33         external diagnostics_is_on         external diagnostics_is_on
34           _RL temparr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy)
35    
36         im1 = 1-OLx         im1 = 1-OLx
37         im2 = sNx+OLx         im2 = sNx+OLx
# Line 59  c Line 60  c
60  c  c
61          enddo          enddo
62          enddo          enddo
63    
64         endif         endif
65    
66  #endif  #endif
67    
68  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
69         if(useDiagnostics) then         if(useDiagnostics) then
70          IF ( DIAGNOSTICS_IS_ON('DPPHYS  ',myThid) ) THEN         IF ( DIAGNOSTICS_IS_ON('DPPHYS  ',myThid) ) THEN
71           CALL DIAGNOSTICS_FILL(dpphys, 'DPPHYS  ',0,Nrphys,0,1,1,myThid)          do bj = myByLo(myThid), myByHi(myThid)
72          ENDIF          do bi = myBxLo(myThid), myBxHi(myThid)
73             do L = 1,Nrphys
74              do j = j1,j2
75              do i = i1,i2
76               temparr(i,j,L,bi,bj) = dpphys(i,j,L,bi,bj) * 1. _d -2
77              enddo
78              enddo
79             enddo
80            enddo
81            enddo
82            CALL DIAGNOSTICS_FILL(temparr,'DPPHYS  ',-1,Nrphys,0,1,1,myThid)
83           ENDIF
84         endif         endif
85  #endif  #endif
86    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22