Changeset 3402 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Oct 16, 2018, 2:35:21 PM (6 years ago)
Author:
lguez
Message:

Set temperature of appearing sea ice to -1.8 Celsius degrees instead
of 0 Celsius degrees, following advice of F. Codron.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

    r3392 r3402  
    32273227                ! F. Codron sensible default values for ocean and sea ice
    32283228                IF (nsrf.EQ.is_oce) THEN
    3229                     tsurf(i,nsrf) = 271.35 ! Freezing sea water
    3230                     DO k=1,nsw
     3229                   tsurf(i,nsrf) = 271.35
     3230                   ! (temperature of sea water under sea ice, so that
     3231                   ! is also the temperature of appearing sea water)
     3232                   DO k=1,nsw
    32313233                      alb_dir(i,k,nsrf) = 0.06 ! typical Ocean albedo
    32323234                      alb_dif(i,k,nsrf) = 0.06
    3233                     ENDDO
     3235                   ENDDO
    32343236                ELSE IF (nsrf.EQ.is_sic) THEN
    3235                     tsurf(i,nsrf) = 273.15 ! Melting ice
    3236                     DO k=1,nsw
     3237                   tsurf(i,nsrf) = 271.35
     3238                   ! (Temperature at base of sea ice. Surface
     3239                   ! temperature could be higher, up to 0 Celsius
     3240                   ! degrees. We set it to -1.8 Celsius degrees for
     3241                   ! consistency with the ocean slab model.)
     3242                   DO k=1,nsw
    32373243                      alb_dir(i,k,nsrf) = 0.3 ! thin ice
    32383244                      alb_dif(i,k,nsrf) = 0.3
    3239                     ENDDO
     3245                   ENDDO
    32403246                ENDIF
    3241                 ! F. Codron
    32423247             ELSE
    32433248                ! The continents have changed. The new fraction receives the mean sum of the existent fractions
Note: See TracChangeset for help on using the changeset viewer.