Ignore:
Timestamp:
2020-10-14T17:45:46+02:00 (5 months ago)
Author:
techene
Message:

#2385 add some comments and cleaning

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/tests/AM98/MY_SRC/usrdef_istate.F90

    r13503 r13599  
    77   !! User defined : set the initial state of a user configuration 
    88   !!====================================================================== 
    9    !! History :  4.0 ! 2016-03  (S. Flavoni) Original code 
     9   !! History :  4.0  ! 2016-03  (S. Flavoni) Original code 
     10   !!             -   ! 2020-03  (A. Nasser) Shallow Water Eq. configuration 
    1011   !!---------------------------------------------------------------------- 
    1112 
     
    1314   !!  usr_def_istate : initial state in Temperature and salinity 
    1415   !!---------------------------------------------------------------------- 
    15    USE par_oce        ! ocean space and time domain 
     16   USE par_oce        ! ocean space and time domain master parameters 
    1617   USE phycst         ! physical constants 
    1718   ! 
    1819   USE in_out_manager ! I/O manager 
    1920   USE lib_mpp        ! MPP library 
    20     
     21 
    2122   IMPLICIT NONE 
    2223   PRIVATE 
     
    4344      REAL(wp), DIMENSION(jpi,jpj,jpk)     , INTENT(in   ) ::   pdept   ! depth of t-point               [m] 
    4445      REAL(wp), DIMENSION(jpi,jpj,jpk)     , INTENT(in   ) ::   ptmask  ! t-point ocean mask             [m] 
     46      ! 
    4547      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(  out) ::   pts     ! T & S fields      [Celsius ; g/kg] 
    4648      REAL(wp), DIMENSION(jpi,jpj,jpk)     , INTENT(  out) ::   pu      ! i-component of the velocity  [m/s]  
     
    5557      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   Ocean at rest, with an horizontally uniform T and S profiles' 
    5658      ! 
    57       pu  (:,:,:) = 0._wp        ! ocean at rest 
    58       pv  (:,:,:) = 0._wp 
    59       pssh(:,:)   = 0._wp 
    60 !!an    
    61        IF(lwp) WRITE(numout,*) "end istate" 
    62         call flush(numout) 
     59      pu  (:,:,:)   = 0._wp        ! ocean at rest 
     60      pv  (:,:,:)   = 0._wp 
     61      pssh(:,:)     = 0._wp 
     62      ts  (:,:,:,:) = 0._wp            ! not used in SWE 
    6363 
     64      IF(lwp) WRITE(numout,*) "end istate" 
     65      CALL FLUSH(numout) 
    6466      !    
    6567   END SUBROUTINE usr_def_istate 
Note: See TracChangeset for help on using the changeset viewer.