| [325] | 1 | MODULE istate |
|---|
| 2 | !!====================================================================== |
|---|
| 3 | !! *** MODULE istate *** |
|---|
| [2528] | 4 | !! Ocean state : initial state setting, off-line case |
|---|
| [325] | 5 | !!===================================================================== |
|---|
| [2528] | 6 | !! History : 3.3 ! 2010-10 (C. Ethe) original code |
|---|
| 7 | !!---------------------------------------------------------------------- |
|---|
| [325] | 8 | |
|---|
| 9 | !!---------------------------------------------------------------------- |
|---|
| [2528] | 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 | |
|---|
| [2528] | 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 | !!---------------------------------------------------------------------- |
|---|
| [2528] | 24 | !! NEMO/OFF 3.3 , NEMO Consortium (2010) |
|---|
| 25 | !! $Id$ |
|---|
| 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 | !! |
|---|
| [2528] | 34 | !! ** Purpose : Initialization to zero of the dynamics and tracers. |
|---|
| [325] | 35 | !!---------------------------------------------------------------------- |
|---|
| [2528] | 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 |
|---|
| [2528] | 47 | ! |
|---|
| [325] | 48 | END SUBROUTINE istate_init |
|---|
| 49 | |
|---|
| 50 | !!===================================================================== |
|---|
| 51 | END MODULE istate |
|---|
