Changeset 4270


Ignore:
Timestamp:
2013-11-20T10:01:12+01:00 (8 years ago)
Author:
clem
Message:

few rewritings limcat_1D

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90

    r4161 r4270  
    5151      !-------------------------------------------------------------------- 
    5252      ijpij = SIZE(zhti,1) 
    53       zht_i(1:ijpij,1:jpl) = 0.d0 
    54       zht_s(1:ijpij,1:jpl) = 0.d0 
    55       za_i (1:ijpij,1:jpl) = 0.d0 
     53      zht_i(1:ijpij,1:jpl) = 0._wp 
     54      zht_s(1:ijpij,1:jpl) = 0._wp 
     55      za_i (1:ijpij,1:jpl) = 0._wp 
    5656 
    5757      !------------------------------------------------------------------------------------ 
     
    8787             
    8888            ! initialisation of ice variables for each try 
    89             zht_i(ji,1:jpl) = 0.d0 
    90             za_i (ji,1:jpl) = 0.d0 
     89            zht_i(ji,1:jpl) = 0._wp 
     90            za_i (ji,1:jpl) = 0._wp 
    9191             
    9292            ! *** case very thin ice: fill only category 1 
     
    100100               ! Fill ice thicknesses except the last one (i_fill) by (hmax-hmin)/2  
    101101               DO jl = 1, i_fill - 1 
    102                   zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) / 2. 
     102                  zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) * 0.5_wp 
    103103               END DO 
    104104                
     
    116116               DO jl = 1, i_fill - 1 
    117117                  IF ( jl == jl0 ) CYCLE 
    118                   zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) / 2. ) 
     118                  zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) * 0.5_wp ) 
    119119                  za_i(ji,jl) =   za_i (ji,jl0) * EXP(-zarg**2) 
    120120               END DO 
     
    146146            ztest_4 = 1 
    147147            DO jl = 1, i_fill 
    148                IF ( za_i(ji,jl) < 0.0d0 ) ztest_4 = 0 
     148               IF ( za_i(ji,jl) < 0._wp ) ztest_4 = 0 
    149149            END DO 
    150150             
Note: See TracChangeset for help on using the changeset viewer.