Ignore:
Timestamp:
2017-09-05T19:53:41+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part2 -

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_da.F90

    r8486 r8498  
    3232CONTAINS 
    3333 
    34 !!gm  even comment line of more than 130 character may cause compilation errors 
    35 !!gm         ===>>> reformat the text 
    3634   SUBROUTINE ice_thd_da 
    3735      !!------------------------------------------------------------------- 
     
    4644      !!                   W = m1 * (Tw -Tf)**m2                    --- originally from Josberger 1979 --- 
    4745      !!                      (Tw - Tf) = elevation of water temp above freezing 
    48       !!                      m1 and m2 = (1.6e-6 , 1.36) best fit from field experiment near the coast of Prince Patrick Island (Perovich 1983) => static ice 
    49       !!                      m1 and m2 = (3.0e-6 , 1.36) best fit from MIZEX 84 experiment (Maykut and Perovich 1987) => moving ice 
     46      !!                      m1 and m2 = (1.6e-6 , 1.36) best fit from field experiment near the coast of Prince Patrick Island 
     47      !!                                                                                           (Perovich 1983) => static ice 
     48      !!                      m1 and m2 = (3.0e-6 , 1.36) best fit from MIZEX 84 experiment 
     49      !!                                                                                (Maykut and Perovich 1987) => moving ice 
    5050      !! 
    5151      !!                   P = N * pi * D                           --- from Rothrock and Thorndike 1984 --- 
     
    5959      !!                      Astar = 1 / ( 1 - (Dmin/Dmax)**(1/beta) ) 
    6060      !!                      Dmin = minimum floe diameter (recommended to be 8m +- 20%) 
    61       !!                      Dmax = maximum floe diameter (recommended to be 300m, but it does not impact melting much except for Dmax<100m) 
     61      !!                      Dmax = maximum floe diameter (recommended to be 300m, 
     62      !!                                                    but it does not impact melting much except for Dmax<100m) 
    6263      !!                      beta = 1.0 +-20% (recommended value) 
    6364      !!                           = 0.3 best fit for western Fram Strait and Antarctica 
     
    113114         ! --- Calculate reduction of total sea ice concentration --- ! 
    114115         zdfloe = rn_dmin * ( zastar / ( zastar - at_i_1d(ji) ) )**rn_beta         ! Mean floe caliper diameter [m] 
    115          zperi  = at_i_1d(ji) * rpi / ( zcs * zdfloe )                             ! Mean perimeter of the floe = N*pi*D = (A/cs*D^2)*pi*D [m.m-2] 
     116         ! 
     117         zperi  = at_i_1d(ji) * rpi / ( zcs * zdfloe )                             ! Mean perimeter of the floe [m.m-2] 
     118         !                                                                         !    = N*pi*D = (A/cs*D^2)*pi*D 
    116119         zwlat  = zm1 * ( MAX( 0._wp, sst_1d(ji) - ( t_bo_1d(ji) - rt0 ) ) )**zm2  ! Melt speed rate [m/s] 
    117           
     120         ! 
    118121         zda_tot(ji) = MIN( zwlat * zperi * rdt_ice, at_i_1d(ji) )                 ! sea ice concentration decrease (>0) 
    119122       
Note: See TracChangeset for help on using the changeset viewer.