Changeset 8819 for branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.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/sbccpl.F90 (modified) (6 diffs)

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

@@ -150 +150 @@
    INTEGER, PARAMETER ::   jps_bio_chloro = 36        ! MEDUSA chlorophyll surface concentration
    INTEGER, PARAMETER ::   jps_fmdice = 37            ! ice form drag
-   INTEGER, PARAMETER ::   jpsnd      = 37            ! total number of fields sent
+   INTEGER, PARAMETER ::   jps_rough_ice_skin = 38    ! ice skin roughness length
+   INTEGER, PARAMETER ::   jpsnd      = 38            ! total number of fields sent
 
    REAL(wp), PARAMETER :: dms_unit_conv = 1.0e+6      ! Coversion factor to get outgong DMS in standard units for coupling
@@ -164 +165 @@
    END TYPE FLD_C
    ! Send to the atmosphere                           !
-   TYPE(FLD_C) ::   sn_snd_temp, sn_snd_alb, sn_snd_thick, sn_snd_crt, sn_snd_co2, sn_snd_cond, sn_snd_mpnd, sn_snd_sstfrz, sn_snd_thick1, sn_snd_fmd
+   TYPE(FLD_C) ::   sn_snd_temp, sn_snd_alb, sn_snd_thick, sn_snd_crt, sn_snd_co2, sn_snd_cond, sn_snd_mpnd, sn_snd_sstfrz, sn_snd_thick1, sn_snd_icerough
    TYPE(FLD_C) ::   sn_snd_bio_co2, sn_snd_bio_dms, sn_snd_bio_chloro                   
 
@@ -242 +243 @@
       !!
       NAMELIST/namsbc_cpl/  sn_snd_temp, sn_snd_alb   , sn_snd_thick , sn_snd_crt   , sn_snd_co2,     &
-         &                  sn_snd_cond, sn_snd_mpnd  , sn_snd_sstfrz, sn_snd_thick1, sn_snd_fmd,     &
+         &                  sn_snd_cond, sn_snd_mpnd  , sn_snd_sstfrz, sn_snd_thick1, sn_snd_icerough,&
          &                  sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau   , sn_rcv_dqnsdt, sn_rcv_qsr,     &
          &                  sn_rcv_qns , sn_rcv_emp   , sn_rcv_rnf   , sn_rcv_cal   , sn_rcv_iceflx,  &
@@ -315 +316 @@
          WRITE(numout,*)'      meltponds fraction & depth      = ', TRIM(sn_snd_mpnd%cldes  ), ' (', TRIM(sn_snd_mpnd%clcat   ), ')'
          WRITE(numout,*)'      sea surface freezing temp       = ', TRIM(sn_snd_sstfrz%cldes   ), ' (', TRIM(sn_snd_sstfrz%clcat   ), ')'
-         WRITE(numout,*)'      ice formdrag                    = ', TRIM(sn_snd_fmd%cldes   ), ' (', TRIM(sn_snd_fmd%clcat   ), ')'
+         WRITE(numout,*)'      ice formdrag                    = ', TRIM(sn_snd_icerough%cldes   ), ' (', TRIM(sn_snd_icerough%clcat   ), ')'
 
          WRITE(numout,*)'  nn_cplmodel                         = ', nn_cplmodel
@@ -875 +876 @@
       !                                                      ! ------------------------- !
       ssnd(jps_fmdice )%clname = 'OIceFmd'
-      SELECT CASE ( TRIM( sn_snd_fmd%cldes ) )
+      ssnd(jps_rough_ice_skin )%clname = 'ORghIcSk'
+      SELECT CASE ( TRIM( sn_snd_icerough%cldes ) )
          CASE ( 'none' )
             ssnd(jps_fmdice)%laction = .FALSE.
+            ssnd(jps_rough_ice_skin)%laction = .FALSE.
          CASE ( 'ice', 'weighted ice' ) 
             ssnd(jps_fmdice)%laction = .TRUE.
-         CASE default   ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_fmd%cldes' )
+            ssnd(jps_rough_ice_skin)%laction = .TRUE.
+         CASE default   ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_icerough%cldes' )
       END SELECT
       !
@@ -2640 +2644 @@
 #if defined key_cice
       IF( ssnd(jps_fmdice)%laction ) THEN
-          SELECT CASE( sn_snd_fmd%cldes )
+          SELECT CASE( sn_snd_icerough%cldes )
              CASE( 'ice' )
                 ztmp1(:,:) = fmdice(:,:)
              CASE( 'weighted ice' ) 
                 ztmp1(:,:) = fmdice(:,:) * fr_i(:,:)
-             CASE default      ;  CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_fmd%cldes' )
+             CASE default      ;  CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_icerough%cldes' )
           END SELECT
       ENDIF
       IF( ssnd(jps_fmdice)%laction )   CALL cpl_snd( jps_fmdice, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info ) 
+
+      IF( ssnd(jps_rough_ice_skin)%laction ) THEN
+          SELECT CASE( sn_snd_icerough%cldes )
+             CASE( 'ice' )
+                ztmp1(:,:) = rough_ice_skin(:,:)
+             CASE( 'weighted ice' ) 
+                ztmp1(:,:) = rough_ice_skin(:,:) * fr_i(:,:)
+             CASE default      ;  CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_icerough%cldes' )
+          END SELECT
+      ENDIF
+      IF( ssnd(jps_rough_ice_skin)%laction )   CALL cpl_snd( jps_rough_ice_skin, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info ) 
 #endif
       !