Changeset 8933 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

Timestamp:

2017-12-07T12:14:23+01:00 (6 years ago)

Author:

clem

Message:

dev_CNRS_2017: raw commit of the Met-Office specifics for coupling with their Earth system model. It is not yet tested and debuging is expected

File:

• branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90 (modified) (2 diffs)

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

@@ -319 +319 @@
          IF(lwp) WRITE(numout,*) '   nn_monocat forced to 0 as jpl>1, i.e. multi-category case is chosen'
       ENDIF
+      !
+     IF( ln_cpl .AND. nn_cats_cpl /= 1 .AND. nn_cats_cpl /= jpl ) THEN
+        CALL ctl_stop( 'STOP', 'par_init: in coupled mode, nn_cats_cpl should be either 1 or jpl' )
+     ENDIF
 !     IF ( jpl == 1 .AND. nn_monocat == 0 ) THEN
 !        CALL ctl_stop( 'STOP', 'par_init : if jpl=1 then nn_monocat should be between 1 and 4' )
@@ -410 +414 @@
 
       ! SIMIP diagnostics
-      diag_fc_bo(:,:)    = 0._wp ; diag_fc_su(:,:)    = 0._wp
-
-      tau_icebfr(:,:) = 0._wp; ! landfast ice param only (clem: important to keep the init here)
+      diag_fc_bo(:,:) = 0._wp ; diag_fc_su(:,:) = 0._wp
+
+      tau_icebfr(:,:)   = 0._wp   ! landfast ice param only (clem: important to keep the init here)
+      cnd_ice   (:,:,:) = 0._wp   ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling)
       
    END SUBROUTINE diag_set0