[325] | 1 | MODULE istate |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE istate *** |
---|
[2444] | 4 | !! Ocean state : initial state setting, off-line case |
---|
[325] | 5 | !!===================================================================== |
---|
[2444] | 6 | !! History : 3.3 ! 2010-10 (C. Ethe) original code |
---|
| 7 | !!---------------------------------------------------------------------- |
---|
[325] | 8 | |
---|
| 9 | !!---------------------------------------------------------------------- |
---|
[2444] | 10 | !! istate_init : initial state set to zero |
---|
[325] | 11 | !!---------------------------------------------------------------------- |
---|
| 12 | USE oce ! ocean dynamics and active tracers |
---|
| 13 | USE dom_oce ! ocean space and time domain |
---|
| 14 | |
---|
| 15 | IMPLICIT NONE |
---|
| 16 | PRIVATE |
---|
| 17 | |
---|
[2444] | 18 | PUBLIC istate_init ! routine called by step.F90 |
---|
[325] | 19 | |
---|
| 20 | !! * Substitutions |
---|
| 21 | # include "domzgr_substitute.h90" |
---|
| 22 | # include "vectopt_loop_substitute.h90" |
---|
| 23 | !!---------------------------------------------------------------------- |
---|
[2287] | 24 | !! NEMO/OFF 3.3 , NEMO Consortium (2010) |
---|
| 25 | !! $Id$ |
---|
[2444] | 26 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
| 27 | !!--------------------------------------------------------------------- |
---|
[325] | 28 | CONTAINS |
---|
| 29 | |
---|
| 30 | SUBROUTINE istate_init |
---|
| 31 | !!---------------------------------------------------------------------- |
---|
| 32 | !! *** ROUTINE istate_init *** |
---|
| 33 | !! |
---|
[2444] | 34 | !! ** Purpose : Initialization to zero of the dynamics and tracers. |
---|
[325] | 35 | !!---------------------------------------------------------------------- |
---|
[2444] | 36 | ! |
---|
| 37 | ! now fields ! after fields ! |
---|
| 38 | un (:,:,:) = 0._wp ; ua(:,:,:) = 0._wp ! |
---|
| 39 | vn (:,:,:) = 0._wp ; va(:,:,:) = 0._wp ! |
---|
| 40 | wn (:,:,:) = 0._wp ! ! |
---|
| 41 | hdivn(:,:,:) = 0._wp ! ! |
---|
| 42 | tsn (:,:,:,:) = 0._wp ! ! |
---|
| 43 | ! |
---|
[325] | 44 | rhd (:,:,:) = 0.e0 |
---|
| 45 | rhop (:,:,:) = 0.e0 |
---|
| 46 | rn2 (:,:,:) = 0.e0 |
---|
[2444] | 47 | ! |
---|
[325] | 48 | END SUBROUTINE istate_init |
---|
| 49 | |
---|
| 50 | !!===================================================================== |
---|
| 51 | END MODULE istate |
---|