Changeset 8561 for branches/UKMO/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/step.F90

Timestamp:

2017-09-22T17:45:41+02:00 (7 years ago)

Author:

jgraham

Message:

Updates for operational diagnostics:
25h mean diagnostics - bottom temperature (and insitu temp)
Operational foam diagnostics - diaopfoam and DIU routines added.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/step.F90 (modified) (2 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/step.F90

@@ -225 +225 @@
       ENDIF
 
+      ! Cool skin
+      IF ( ln_diurnal ) THEN  
+         IF ( ln_blk_core ) THEN
+            CALL diurnal_sst_coolskin_step( &  
+                    qns(:,:)+(rn_abs*qsr(:,:)), taum, rhop(:,:,1), rdt) 
+         ELSE
+            CALL diurnal_sst_coolskin_step( &  
+                    qns, taum, rhop(:,:,1), rdt) 
+         ENDIF
+      ENDIF
+
       !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
       ! diagnostics and outputs             (ua, va, tsa used as workspace)
@@ -237 +248 @@
       !
       IF( ln_crs     )      CALL crs_fld( kstp )         ! ocean model: online field coarsening & output
+
+      !Diurnal warm layer model        
+      IF ( ln_diurnal ) THEN
+         IF ( ln_blk_core ) THEN
+            IF( kstp == nit000 )THEN  
+               CALL diurnal_sst_takaya_step( &  
+               &    qsr(:,:)-(rn_abs*qsr(:,:)), qns(:,:)+(rn_abs*qsr(:,:)), &
+               &    taum, rhop(:,:,1), &
+               &    rdt, ld_calcfrac = .TRUE.)  
+            ELSE  
+               CALL diurnal_sst_takaya_step( &  
+               &    qsr(:,:)-(rn_abs*qsr(:,:)), qns(:,:)+(rn_abs*qsr(:,:)), &
+               &    taum, rhop(:,:,1), rdt )  
+            ENDIF 
+         ELSE
+            IF( kstp == nit000 )THEN  
+               CALL diurnal_sst_takaya_step( &  
+               &    qsr, qns, taum, rhop(:,:,1), &
+               &    rdt, ld_calcfrac = .TRUE.)  
+            ELSE  
+               CALL diurnal_sst_takaya_step( &  
+               &    qsr, qns, taum, rhop(:,:,1), rdt )  
+            ENDIF 
+         ENDIF
+      ENDIF
 
 #if defined key_top