/[MITgcm]/MITgcm/pkg/aim/aim_dyn2aim.F
ViewVC logotype

Annotation of /MITgcm/pkg/aim/aim_dyn2aim.F

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

Revision 1.3 - (hide annotations) (download)
Wed Aug 15 13:50:43 2001 UTC (22 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: icebear2, checkpoint44h_pre, release1_p12, release1_p10, release1_p11, release1_p16, release1_p15, ecco_c44_e17, ecco_c44_e16, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, icebear5, icebear4, checkpoint44f_pre, icebear3, checkpoint46f_post, release1_p13_pre, checkpoint46d_pre, checkpoint46e_post, release1-branch_tutorials, release1_p14, checkpoint44g_post, checkpoint46h_pre, checkpoint44h_post, release1_p12_pre, checkpoint44e_post, checkpoint46e_pre, ecco-branch-mod4, checkpoint43a-release1mods, checkpoint45d_post, checkpoint45b_post, checkpoint46b_pre, chkpt44a_pre, release1-branch-end, release1_final_v1, ecco_c44_e19, checkpoint46, ecco_c44_e20, checkpoint44, release1_p13, ecco_c44_e18, checkpoint44f_post, checkpoint40pre8, release1_p17, release1_b1, checkpoint44b_post, chkpt44c_post, chkpt44d_post, checkpoint42, release1_p9, release1_p8, checkpoint43, checkpoint46g_pre, release1_p2, release1_p3, release1_p4, release1_p6, checkpoint46a_post, chkpt44a_post, checkpoint44b_pre, release1_p1, checkpoint46a_pre, ecco-branch-mod1, checkpoint45c_post, release1_p5, checkpoint44e_pre, chkpt44c_pre, checkpoint40pre9, release1_p7, ecco_ice2, ecco_ice1, checkpoint46b_post, checkpoint46d_post, ecco-branch-mod2, checkpoint46g_post, checkpoint45a_post, checkpoint46c_pre, ecco-branch-mod3, ecco-branch-mod5, ecco_c44_e22, release1_beta1, ecco_c44_e23, release1-branch_branchpoint, checkpoint46c_post, checkpoint40, checkpoint45, checkpoint46h_post, release1_chkpt44d_post, ecco_c44_e25, checkpoint41
Branch point for: c24_e25_ice, ecco-branch, release1_coupled, icebear, release1_final, release1-branch, release1, release1_50yr
Changes since 1.2: +7 -4 lines 
allow to run AIM with partial cell.
1 jmc 1.3 C $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/aim_dyn2aim.F,v 1.2 2001/06/18 17:39:58 cnh Exp $
2 cnh 1.1 C $Name: $
3    
4     #include "AIM_OPTIONS.h"
5    
6     CStartOfInterface
7     SUBROUTINE AIM_DYN2AIM(
8     I bi,bj, myCurrentTime,myThid)
9     C /==========================================================\
10     C | S/R AIM_DYN2AIM |
11     C | o Map dynamics conforming arrays to AIM internal arrays. |
12     C |==========================================================|
13     C | Currently this routine transfers grid information to |
14     C | AIM. Other information needs to be included from |
15     C | aim_do_atmos_physics(). |
16     C \==========================================================/
17     IMPLICIT rEAL*8 (A-H,O-Z)
18    
19     C == Global data ==
20     #include "SIZE.h"
21     #include "EEPARAMS.h"
22     #include "PARAMS.h"
23     #include "GRID.h"
24     #include "DYNVARS.h"
25    
26     #ifdef ALLOW_AIM
27     #include "atparam0.h"
28     #include "atparam1.h"
29     INTEGER NGP
30     INTEGER NLON
31     INTEGER NLAT
32     INTEGER NLEV
33     PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
34     #include "Lev_def.h"
35     #include "com_physvar.h"
36     #include "com_physcon.h"
37     #endif /* ALLOW_AIM */
38    
39     C == Routine arguments ==
40     INTEGER bi, bj
41     _RL myCurrentTime
42     INTEGER myThid
43     CEndOfInterface
44    
45     #ifdef ALLOW_AIM
46     C == Local variables ==
47     C Loop counters
48     INTEGER I, J, I2, K
49    
50     C Legendre polynomials
51     DO J=1,NLAT
52 cnh 1.2 DO I=1,NLON
53     I2=J
54     I2=NLON*(J-1)+I
55     FMU(I2,1,myThid)=SIN(yC(I,J,bi,bj)*deg2rad)
56     FMU(I2,2,myThid)=1.5 _d 0*FMU(I2,1,myThid)**2-0.5 _d 0
57     ENDDO
58 cnh 1.1 ENDDO
59 cnh 1.2 C DO J=1,NLAT
60     C FMU(J,1,myThid)=SIN(yC(1,J,bi,bj)*deg2rad)
61     C FMU(J,2,myThid)=1.5 _d 0*FMU(J,1,myThid)**2-0.5 _d 0
62     C ENDDO
63 cnh 1.1
64     C Depth in levels
65     I2=0
66     DO J=1,NLAT
67     DO I=1,NLON
68     I2=I2+1
69     NLEVxy(I2,myThid) =0
70     NLEVxyU(I2,myThid)=0
71     NLEVxyV(I2,myThid)=0
72     DO k=1,NLEV
73 jmc 1.3 NLEVxy(I2,myThid) = NLEVxy(I2,myThid)
74     & + INT( maskC(I,J,K,bi,bj) )
75     NLEVxyU(I2,myThid)= NLEVxyU(I2,myThid)
76     & + INT( maskW(I,J,K,bi,bj) )
77     NLEVxyV(I2,myThid)= NLEVxyV(I2,myThid)
78     & + INT( maskS(I,J,K,bi,bj) )
79 cnh 1.1 ENDDO
80     ENDDO
81     ENDDO
82    
83     #endif /* ALLOW_AIM */
84    
85     RETURN
86     END
  ViewVC Help
Powered by ViewVC 1.1.22