Changeset 11180 for NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DOM/istate.F90

Timestamp:

2019-06-25T18:50:27+02:00 (5 years ago)

Author:

clne

Message:

Initial commit of code for 2d (surge) work in NEMO4.
This is aiming to replicate the 3.6 version in branches/UKMO/dev_r5518_Surge_Modelling

File:

• NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DOM/istate.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DOM/istate.F90

@@ -96 +96 @@
          !                                    ! Initialization of ocean to zero
          !
-         IF( ln_tsd_init ) THEN               
+         IF (ln_2d) THEN 
+            IF(lwp) WRITE(numout,*) 'istate_init : 2D case, setting tracers to contants and ocean at rest' 
+            tsb(:,:,:,:)= 15._wp               ! No tracers, but can't set salinity to 0 otherwise it triggers crash message
+            sshb(:,:)   = 0._wp                ! set the ocean at rest
+            ub  (:,:,:) = 0._wp                
+            vb  (:,:,:) = 0._wp           
+         ELSE IF( ln_tsd_init ) THEN   
             CALL dta_tsd( nit000, tsb )       ! read 3D T and S data at nit000
             !