Ignore:
Timestamp:
2017-11-27T14:24:20+01:00 (3 years ago)
Author:
jamrae
Message:

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

File:
1 edited

Legend:

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

    r8670 r8819  
    6868#endif 
    6969   USE ice_forcing, only: frcvdr,frcvdf,frcidr,frcidf 
    70    USE ice_atmo, only: calc_strair, z0atm 
     70   USE ice_atmo, only: calc_strair, z0atm, z0atmskin 
    7171 
    7272   USE CICE_InitMod 
     
    187187      CALL CICE_Initialize 
    188188 
    189       IF(lwp) THEN                        ! control print 
    190         WRITE(numout,*) 
    191         WRITE(numout,*)'=============================================================' 
    192         WRITE(numout,*)'In cice_sbc_init, after call to CICE_Initialize' 
    193         WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm) 
    194         WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm) 
    195         WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm) 
    196         WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm) 
    197         WRITE(numout,*)'=============================================================' 
    198         WRITE(numout,*) 
    199       ENDIF 
    200  
    201189      ! Do some CICE consistency checks 
    202190      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 
     
    272260 
    273261      CALL cice2nemo(z0atm,fmdice,'T', 1. ) 
     262      CALL cice2nemo(z0atmskin,rough_ice_skin,'T', 1. ) 
    274263 
    275264      IF( .NOT. ln_rstart ) THEN 
     
    323312      IF( nn_timing == 1 )  CALL timing_stop('cice_sbc_init') 
    324313      ! 
     314 
     315 
     316      IF(lwp) THEN                        ! control print 
     317        WRITE(numout,*) 
     318        WRITE(numout,*)'=============================================================' 
     319        WRITE(numout,*)'At end of subroutine cice_sbc_init' 
     320        WRITE(numout,*)'MAXVAL(z0atm): ', MAXVAL(z0atm) 
     321        WRITE(numout,*)'MAXLOC(z0atm): ', MAXLOC(z0atm) 
     322        WRITE(numout,*)'MINVAL(z0atm): ', MINVAL(z0atm) 
     323        WRITE(numout,*)'MINLOC(z0atm): ', MINLOC(z0atm) 
     324        WRITE(numout,*)'MAXVAL(z0atmskin): ', MAXVAL(z0atmskin) 
     325        WRITE(numout,*)'MAXLOC(z0atmskin): ', MAXLOC(z0atmskin) 
     326        WRITE(numout,*)'MINVAL(z0atmskin): ', MINVAL(z0atmskin) 
     327        WRITE(numout,*)'MINLOC(z0atmskin): ', MINLOC(z0atmskin) 
     328        WRITE(numout,*)'=============================================================' 
     329        WRITE(numout,*) 
     330      ENDIF 
     331 
     332 
     333 
    325334   END SUBROUTINE cice_sbc_init 
    326335 
     
    803812      CALL cice2nemo(z0atm,fmdice,'T', 1. ) 
    804813 
     814! Surface skin roughness length 
     815      CALL cice2nemo(z0atmskin,rough_ice_skin,'T', 1. ) 
     816 
    805817      ! 
    806818      IF( nn_timing == 1 )  CALL timing_stop('cice_sbc_hadgam') 
Note: See TracChangeset for help on using the changeset viewer.