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

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

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

Revision 1.10 - (show annotations) (download)
Thu Aug 3 19:03:15 2006 UTC (17 years, 9 months ago) by molod
Branch: MAIN
CVS Tags: checkpoint58u_post, checkpoint58w_post, checkpoint63p, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint58r_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58q_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint59, checkpoint58o_post, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint58v_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.9: +3 -2 lines 
use _d notation
1 C $Header: /u/gcmpack/MITgcm/pkg/gridalt/gridalt_update.F,v 1.9 2005/07/07 23:18:51 molod Exp $
2 C $Name: $
3
4 subroutine gridalt_update(myThid)
5 c***********************************************************************
6 c subroutine gridalt_update
7 c
8 c Purpose: Update the pressure thicknesses of the layers of the
9 c alternative vertical grid (used now for atmospheric physics).
10 c
11 c Calculate: dpphys - new delta r (p*) edge-to-edge of physics grid
12 c using dpphys0 (initial value) and rstarfacC
13 c
14 c***********************************************************************
15 implicit none
16 c
17 #include "PACKAGES_CONFIG.h"
18 #include "CPP_OPTIONS.h"
19 #include "SIZE.h"
20 #include "EEPARAMS.h"
21 #include "PARAMS.h"
22 #include "GRID.h"
23 #include "fizhi_SIZE.h"
24 #include "gridalt_mapping.h"
25 #include "SURFACE.h"
26
27 integer myThid
28
29 integer im1,im2,jm1,jm2
30 integer i1,i2,j1,j2,bi,bj
31 integer i,j,L
32 logical diagnostics_is_on
33 external diagnostics_is_on
34 _RL temparr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy)
35
36 im1 = 1-OLx
37 im2 = sNx+OLx
38 jm1 = 1-OLy
39 jm2 = sNy+OLy
40 i1 = 1
41 i2 = sNx
42 j1 = 1
43 j2 = sNy
44
45 #ifdef ALLOW_FIZHI
46 #ifdef NONLIN_FRSURF
47
48 if(usefizhi .and. (select_rStar.GT.0) )then
49 do bj = myByLo(myThid), myByHi(myThid)
50 do bi = myBxLo(myThid), myBxHi(myThid)
51 c
52 do L = 1,Nrphys
53 do j = j1,j2
54 do i = i1,i2
55 dpphys(i,j,L,bi,bj) = dpphys0(i,j,L,bi,bj) *
56 . rstarfacC(i,j,bi,bj)
57 enddo
58 enddo
59 enddo
60 c
61 enddo
62 enddo
63
64 endif
65
66 #endif
67
68 #ifdef ALLOW_DIAGNOSTICS
69 if(useDiagnostics) then
70 IF ( DIAGNOSTICS_IS_ON('DPPHYS ',myThid) ) THEN
71 do bj = myByLo(myThid), myByHi(myThid)
72 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
85 #endif
86
87 #endif
88 return
89 end
  ViewVC Help
Powered by ViewVC 1.1.22