Changeset 8819 for branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2017-11-27T14:24:20+01:00 (6 years ago)

Author:

jamrae

Message:

Implemented code to pass sea ice skin roughness length to atmosphere.

File:

• branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (5 diffs)

branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -68 +68 @@
 #endif
    USE ice_forcing, only: frcvdr,frcvdf,frcidr,frcidf
-   USE ice_atmo, only: calc_strair, z0atm
+   USE ice_atmo, only: calc_strair, z0atm, z0atmskin
 
    USE CICE_InitMod
@@ -187 +187 @@
       CALL CICE_Initialize
 
-      IF(lwp) THEN                        ! control print
-        WRITE(numout,*)
-        WRITE(numout,*)'============================================================='
-        WRITE(numout,*)'In cice_sbc_init, after call to CICE_Initialize'
-        WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm)
-        WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm)
-        WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm)
-        WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm)
-        WRITE(numout,*)'============================================================='
-        WRITE(numout,*)
-      ENDIF
-
       ! Do some CICE consistency checks
       IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN
@@ -272 +260 @@
 
       CALL cice2nemo(z0atm,fmdice,'T', 1. )
+      CALL cice2nemo(z0atmskin,rough_ice_skin,'T', 1. )
 
       IF( .NOT. ln_rstart ) THEN
@@ -323 +312 @@
       IF( nn_timing == 1 )  CALL timing_stop('cice_sbc_init')
       !
+
+
+      IF(lwp) THEN                        ! control print
+        WRITE(numout,*)
+        WRITE(numout,*)'============================================================='
+        WRITE(numout,*)'At end of subroutine cice_sbc_init'
+        WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm)
+        WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm)
+        WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm)
+        WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm)
+        WRITE(numout,*)'MAXVAL(z0atmskin): ', MAXVAL(z0atmskin)
+        WRITE(numout,*)'MAXLOC(z0atmskin): ', MAXLOC(z0atmskin)
+        WRITE(numout,*)'MINVAL(z0atmskin): ', MINVAL(z0atmskin)
+        WRITE(numout,*)'MINLOC(z0atmskin): ', MINLOC(z0atmskin)
+        WRITE(numout,*)'============================================================='
+        WRITE(numout,*)
+      ENDIF
+
+
+
    END SUBROUTINE cice_sbc_init
 
@@ -803 +812 @@
       CALL cice2nemo(z0atm,fmdice,'T', 1. )
 
+! Surface skin roughness length
+      CALL cice2nemo(z0atmskin,rough_ice_skin,'T', 1. )
+
       !
       IF( nn_timing == 1 )  CALL timing_stop('cice_sbc_hadgam')