Changeset 8670 for branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2017-11-02T16:02:02+01:00 (6 years ago)

Author:

jamrae

Message:

Added capability to weight sea ice roughness length before coupling.

File:

• branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (3 diffs)

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

@@ -878 +878 @@
          CASE ( 'none' )
             ssnd(jps_fmdice)%laction = .FALSE.
-         CASE ( 'coupled' ) 
+         CASE ( 'ice', 'weighted ice' ) 
             ssnd(jps_fmdice)%laction = .TRUE.
          CASE default   ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_fmd%cldes' )
       END SELECT
       !
-      IF(lwp) THEN                        ! control print
-         WRITE(numout,*)
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)'In sbc_cpl_init, after setting up ssnd(jps_fmdice)'
-         WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname
-         WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)
-      ENDIF
       !
       !                                                      ! ------------------------------- !
@@ -970 +961 @@
       ! ================================ !
 
-      IF(lwp) THEN                        ! control print
-         WRITE(numout,*)
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)'In sbc_cpl_init, before call to cpl_define'
-         WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname
-         WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)
-      ENDIF
-
       CALL cpl_define(jprcv, jpsnd, nn_cplmodel)
-      
-
-      IF(lwp) THEN                        ! control print
-         WRITE(numout,*)
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)'In sbc_cpl_init, after call to cpl_define'
-         WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname
-         WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)
-      ENDIF
 
       IF (ln_usecplmask) THEN 
@@ -2658 +2628 @@
       IF( ssnd(jps_rnf   )%laction )  CALL cpl_snd( jps_rnf   , isec, RESHAPE ( rnf , (/jpi,jpj,1/) ), info )
       IF( ssnd(jps_taum  )%laction )  CALL cpl_snd( jps_taum  , isec, RESHAPE ( taum, (/jpi,jpj,1/) ), info )
-      
-      !
-      IF(lwp) THEN                        ! control print
-         WRITE(numout,*)
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)'In sbc_cpl_snd'
-         WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname
-         WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction
-         WRITE(numout,*)'============================================================='
-         WRITE(numout,*)
-      ENDIF
-      !                                  
 
 #if defined key_cice
       ztmp1(:,:) = sstfrz(:,:) + rt0
       IF( ssnd(jps_sstfrz)%laction )  CALL cpl_snd( jps_sstfrz, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info )
-      ztmp1(:,:) = fmdice(:,:)
+#endif
+      !
+      !      
+      !                                                      ! ------------------------- !
+      !                                                      !   Sea ice form drag       !
+      !                                                      ! ------------------------- !
+#if defined key_cice
+      IF( ssnd(jps_fmdice)%laction ) THEN
+          SELECT CASE( sn_snd_fmd%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' )
+          END SELECT
+      ENDIF
       IF( ssnd(jps_fmdice)%laction )   CALL cpl_snd( jps_fmdice, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info ) 
 #endif