Changeset 8670


Ignore:
Timestamp:
2017-11-02T16:02:02+01:00 (3 years ago)
Author:
jamrae
Message:

Added capability to weight sea ice roughness length before coupling.

Location:
branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r8460 r8670  
    878878         CASE ( 'none' ) 
    879879            ssnd(jps_fmdice)%laction = .FALSE. 
    880          CASE ( 'coupled' )  
     880         CASE ( 'ice', 'weighted ice' )  
    881881            ssnd(jps_fmdice)%laction = .TRUE. 
    882882         CASE default   ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_fmd%cldes' ) 
    883883      END SELECT 
    884884      ! 
    885       IF(lwp) THEN                        ! control print 
    886          WRITE(numout,*) 
    887          WRITE(numout,*)'=============================================================' 
    888          WRITE(numout,*)'In sbc_cpl_init, after setting up ssnd(jps_fmdice)' 
    889          WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname 
    890          WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction 
    891          WRITE(numout,*)'=============================================================' 
    892          WRITE(numout,*) 
    893       ENDIF 
    894885      ! 
    895886      !                                                      ! ------------------------------- ! 
     
    970961      ! ================================ ! 
    971962 
    972       IF(lwp) THEN                        ! control print 
    973          WRITE(numout,*) 
    974          WRITE(numout,*)'=============================================================' 
    975          WRITE(numout,*)'In sbc_cpl_init, before call to cpl_define' 
    976          WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname 
    977          WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction 
    978          WRITE(numout,*)'=============================================================' 
    979          WRITE(numout,*) 
    980       ENDIF 
    981  
    982963      CALL cpl_define(jprcv, jpsnd, nn_cplmodel) 
    983        
    984  
    985       IF(lwp) THEN                        ! control print 
    986          WRITE(numout,*) 
    987          WRITE(numout,*)'=============================================================' 
    988          WRITE(numout,*)'In sbc_cpl_init, after call to cpl_define' 
    989          WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname 
    990          WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction 
    991          WRITE(numout,*)'=============================================================' 
    992          WRITE(numout,*) 
    993       ENDIF 
    994964 
    995965      IF (ln_usecplmask) THEN  
     
    26582628      IF( ssnd(jps_rnf   )%laction )  CALL cpl_snd( jps_rnf   , isec, RESHAPE ( rnf , (/jpi,jpj,1/) ), info ) 
    26592629      IF( ssnd(jps_taum  )%laction )  CALL cpl_snd( jps_taum  , isec, RESHAPE ( taum, (/jpi,jpj,1/) ), info ) 
    2660        
    2661       ! 
    2662       IF(lwp) THEN                        ! control print 
    2663          WRITE(numout,*) 
    2664          WRITE(numout,*)'=============================================================' 
    2665          WRITE(numout,*)'In sbc_cpl_snd' 
    2666          WRITE(numout,*)'ssnd(jps_fmdice )%clname: ', ssnd(jps_fmdice )%clname 
    2667          WRITE(numout,*)'ssnd(jps_fmdice)%laction: ',ssnd(jps_fmdice)%laction 
    2668          WRITE(numout,*)'=============================================================' 
    2669          WRITE(numout,*) 
    2670       ENDIF 
    2671       !                                   
    26722630 
    26732631#if defined key_cice 
    26742632      ztmp1(:,:) = sstfrz(:,:) + rt0 
    26752633      IF( ssnd(jps_sstfrz)%laction )  CALL cpl_snd( jps_sstfrz, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info ) 
    2676       ztmp1(:,:) = fmdice(:,:) 
     2634#endif 
     2635      ! 
     2636      !       
     2637      !                                                      ! ------------------------- ! 
     2638      !                                                      !   Sea ice form drag       ! 
     2639      !                                                      ! ------------------------- ! 
     2640#if defined key_cice 
     2641      IF( ssnd(jps_fmdice)%laction ) THEN 
     2642          SELECT CASE( sn_snd_fmd%cldes ) 
     2643             CASE( 'ice' ) 
     2644                ztmp1(:,:) = fmdice(:,:) 
     2645             CASE( 'weighted ice' )  
     2646                ztmp1(:,:) = fmdice(:,:) * fr_i(:,:) 
     2647             CASE default      ;  CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_fmd%cldes' ) 
     2648          END SELECT 
     2649      ENDIF 
    26772650      IF( ssnd(jps_fmdice)%laction )   CALL cpl_snd( jps_fmdice, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info )  
    26782651#endif 
  • branches/UKMO/dev_r5518_GO6_package_r8356_plus_form_drag/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r8484 r8670  
    271271      ENDIF 
    272272 
    273       IF(lwp) THEN                        ! control print 
    274         WRITE(numout,*) 
    275         WRITE(numout,*)'=============================================================' 
    276         WRITE(numout,*)'In cice_sbc_init, before call to cice2nemo for z0atm/fmdice' 
    277         WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm) 
    278         WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm) 
    279         WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm) 
    280         WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm) 
    281         WRITE(numout,*)'=============================================================' 
    282         WRITE(numout,*) 
    283       ENDIF 
    284  
    285273      CALL cice2nemo(z0atm,fmdice,'T', 1. ) 
    286  
    287       IF(lwp) THEN                        ! control print 
    288         WRITE(numout,*) 
    289         WRITE(numout,*)'=============================================================' 
    290         WRITE(numout,*)'In cice_sbc_init, after call to cice2nemo for z0atm/fmdice' 
    291         WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm) 
    292         WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm) 
    293         WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm) 
    294         WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm) 
    295         WRITE(numout,*)'MAXVAL(fmdice): ', MAXVAL(fmdice) 
    296         WRITE(numout,*)'MAXLOC(fmdice): ', MAXLOC(fmdice) 
    297         WRITE(numout,*)'MINVAL(fmdice): ', MINVAL(fmdice) 
    298         WRITE(numout,*)'MINLOC(fmdice): ', MINLOC(fmdice) 
    299         WRITE(numout,*)'============================================================='  
    300         WRITE(numout,*) 
    301       ENDIF 
    302274 
    303275      IF( .NOT. ln_rstart ) THEN 
Note: See TracChangeset for help on using the changeset viewer.