/[MITgcm]/MITgcm/pkg/seaice/groatb.F
ViewVC logotype

Annotation of /MITgcm/pkg/seaice/groatb.F

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

Revision 1.5 - (hide annotations) (download)
Tue Feb 18 05:33:55 2003 UTC (22 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint48i_post, checkpoint48h_post, checkpoint50b_pre, checkpoint49, checkpoint48g_post, checkpoint50, checkpoint50b_post, checkpoint50a_post
Changes since 1.4: +2 -2 lines 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 heimbach 1.2 C $Header:
2    
3     #include "SEAICE_OPTIONS.h"
4    
5     CStartOfInterface
6     SUBROUTINE groatb( A22, myThid )
7     C /==========================================================\
8     C | SUBROUTINE groatb |
9     C | o Calculate ice growth |
10     C |==========================================================|
11     C \==========================================================/
12     IMPLICIT NONE
13    
14     C === Global variables ===
15     #include "SIZE.h"
16     #include "EEPARAMS.h"
17     #include "PARAMS.h"
18     #include "DYNVARS.h"
19     #include "FFIELDS.h"
20     #include "SEAICE.h"
21     #include "SEAICE_FFIELDS.h"
22    
23     C === Routine arguments ===
24     C myThid - Thread no. that called this routine.
25     _RL A22
26     INTEGER myThid
27     CEndOfInterface
28    
29     #ifdef ALLOW_SEAICE
30    
31     C === Local variables ===
32 dimitri 1.4 C i,j,bi,bj - Loop counters
33 heimbach 1.2
34 dimitri 1.4 INTEGER i, j, bi, bj
35 heimbach 1.2 INTEGER KOPEN
36     _RL U1, V1
37    
38     _RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
39     _RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
40    
41     C if this is not done before then do it here
42 dimitri 1.4 _EXCH_XY_R8(UWIND, myThid)
43     _EXCH_XY_R8(VWIND, myThid)
44 heimbach 1.2
45 dimitri 1.3 c$taf loop = parallel
46 heimbach 1.2 DO bj=myByLo(myThid),myByHi(myThid)
47 dimitri 1.3 c$taf loop = parallel
48 heimbach 1.2 DO bi=myBxLo(myThid),myBxHi(myThid)
49    
50     C DETERMINE AMOUNT OF OPEN WATER AND ICE THICKNESS
51     DO J=1,sNy
52     DO I=1,sNx
53     AREA(I,J,2,bi,bj)=MAX(A22,AREA(I,J,2,bi,bj))
54 dimitri 1.3 FHEFF(I,J,bi,bj)=0.0 _d 0
55 heimbach 1.2 HICE(I,J)=HEFF(I,J,2,bi,bj)/AREA(I,J,2,bi,bj)
56     ENDDO
57     ENDDO
58    
59     C NOW DETERMINE MIXED LAYER TEMPERATURE
60     DO J=1,sNy
61     DO I=1,sNx
62 dimitri 1.3 TMIX(I,J,bi,bj)=theta(I,J,1,bi,bj)+273.16 _d +00
63 heimbach 1.2 #ifdef SEAICE_DEBUG
64 dimitri 1.3 TMIX(I,J,bi,bj)=MAX(TMIX(I,J,bi,bj),271.2 _d +00)
65 dimitri 1.5 #endif
66 heimbach 1.2 ENDDO
67     ENDDO
68    
69     DO J=1,sNy
70     DO I=1,sNx
71 dimitri 1.4 U1=0.25 _d 0*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj)
72     1 +UWIND(I,J-1,bi,bj)+UWIND(I,J,bi,bj))
73     V1=0.25 _d 0*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj)
74     1 +VWIND(I,J-1,bi,bj)+VWIND(I,J,bi,bj))
75 heimbach 1.2 UG(I,J)=SQRT(U1**2+V1**2)
76     ENDDO
77     ENDDO
78    
79     C NOW DETERMINE GROWTH RATES
80     C FIRST DO OPEN WATER
81     KOPEN=-1
82     CALL BUDGET(UG, TMIX, HICE, FO, KOPEN, bi, bj)
83     C NOW DO ICE
84     KOPEN=1
85     CALL BUDGET(UG, TICE, HICE, FICE, KOPEN, bi, bj)
86    
87     ENDDO
88     ENDDO
89    
90 dimitri 1.5 #endif /* ALLOW_SEAICE */
91 heimbach 1.2
92     RETURN
93     END
  ViewVC Help
Powered by ViewVC 1.1.22