1 | MODULE istate |
---|
2 | !!====================================================================== |
---|
3 | !! *** MODULE istate *** |
---|
4 | !! Ocean state : initial state setting |
---|
5 | !!===================================================================== |
---|
6 | |
---|
7 | !!---------------------------------------------------------------------- |
---|
8 | !! istate_init : initial state setting |
---|
9 | !!---------------------------------------------------------------------- |
---|
10 | !! * Modules used |
---|
11 | USE oce ! ocean dynamics and active tracers |
---|
12 | USE dom_oce ! ocean space and time domain |
---|
13 | USE ldftra_oce ! ocean active tracers: lateral physics |
---|
14 | USE zdf_oce ! ocean vertical physics |
---|
15 | USE in_out_manager ! I/O manager |
---|
16 | USE phycst ! physical constants |
---|
17 | |
---|
18 | IMPLICIT NONE |
---|
19 | PRIVATE |
---|
20 | |
---|
21 | !! * Routine accessibility |
---|
22 | PUBLIC istate_init ! routine called by step.F90 |
---|
23 | |
---|
24 | !! * Substitutions |
---|
25 | # include "domzgr_substitute.h90" |
---|
26 | # include "vectopt_loop_substitute.h90" |
---|
27 | !!---------------------------------------------------------------------- |
---|
28 | !! OPA 9.0 , LOCEAN-IPSL (2005) |
---|
29 | !! $Id$ |
---|
30 | !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
31 | !!---------------------------------------------------------------------- |
---|
32 | |
---|
33 | CONTAINS |
---|
34 | |
---|
35 | SUBROUTINE istate_init |
---|
36 | !!---------------------------------------------------------------------- |
---|
37 | !! *** ROUTINE istate_init *** |
---|
38 | !! |
---|
39 | !! ** Purpose : Initialization of the dynamics and tracers. |
---|
40 | !! |
---|
41 | !! ** Method : |
---|
42 | !! |
---|
43 | !! History : |
---|
44 | !! 4.0 ! 91-03 () Original code |
---|
45 | !! ! 91-11 (G. Madec) |
---|
46 | !! 9.0 ! 03-09 (G. Madec) F90: Free form, modules, orthogonality |
---|
47 | !!---------------------------------------------------------------------- |
---|
48 | !! * Local declarations |
---|
49 | !!---------------------------------------------------------------------- |
---|
50 | |
---|
51 | |
---|
52 | ! Initialization to zero |
---|
53 | ! ---------------------- |
---|
54 | |
---|
55 | ! before fields ! now fields ! after fields ! |
---|
56 | ; un (:,:,:) = 0.e0 ; ua (:,:,:) = 0.e0 |
---|
57 | ; vn (:,:,:) = 0.e0 ; va (:,:,:) = 0.e0 |
---|
58 | ; ; wn (:,:,:) = 0.e0 |
---|
59 | ; hdivn(:,:,:) = 0.e0 ; |
---|
60 | |
---|
61 | ; tn (:,:,:) = 0.e0 ; ta (:,:,:) = 0.e0 |
---|
62 | ; sn (:,:,:) = 0.e0 ; sa (:,:,:) = 0.e0 |
---|
63 | |
---|
64 | rhd (:,:,:) = 0.e0 |
---|
65 | rhop (:,:,:) = 0.e0 |
---|
66 | rn2 (:,:,:) = 0.e0 |
---|
67 | |
---|
68 | |
---|
69 | END SUBROUTINE istate_init |
---|
70 | |
---|
71 | !!===================================================================== |
---|
72 | END MODULE istate |
---|