Ignore:
Timestamp:
Jul 24, 2024, 4:23:34 PM (7 months ago)
Author:
abarral
Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/inlandsis.f90

    r5116 r5117  
    183183  IMPLICIT NONE
    184184
    185   logical :: SnoMod
    186   logical :: BloMod
    187   logical :: debut
     185  LOGICAL :: SnoMod
     186  LOGICAL :: BloMod
     187  LOGICAL :: debut
    188188  INTEGER :: jjtime
    189189
    … …  
    409409    IF (BloMod) THEN
    410410
    411     if (klonv==1) THEN
    412       IF(isnoSV(1)>=2                   .and. &
    413             TsisSV(1,max(1,isnoSV(1)))<273.  .and. &
    414             ro__SV(1,max(1,isnoSV(1)))<500.  .and. &
     411    IF (klonv==1) THEN
     412      IF(isnoSV(1)>=2                   .AND. &
     413            TsisSV(1,max(1,isnoSV(1)))<273.  .AND. &
     414            ro__SV(1,max(1,isnoSV(1)))<500.  .AND. &
    415415            eta_SV(1,max(1,isnoSV(1)))<epsi) THEN
    416416  ! +                       **********
    … …  
    518518  ! +     Gallee et al., 2001    eq 5, p5
    519519
    520       if (usth_param == "gal") THEN
     520      IF (usth_param == "gal") THEN
    521521        Salt_us   =   (log(2.868) - log(1 + SaltMo)) * rCd10n/0.085
    522522        Salt_us   = Salt_us * Fac_Mo
    … …  
    525525      endif
    526526
    527       if (usth_param == "lis") then !Liston et al. 2007
     527      IF (usth_param == "lis") then !Liston et al. 2007
    528528        IF(ro__SV(ikl,isn)>300.) THEN
    529529          Salt_us   = 0.005*exp(0.013*ro__SV(ikl,isn))
    … …  
    607607    hSalSV(ikl) = 8.436e-2  * us__SV(ikl)**SblPom
    608608
    609     if (qsalt_param == "pom") THEN
     609    IF (qsalt_param == "pom") THEN
    610610      qSalSV(ikl) = (us__SV(ikl)**2 - usthSV(ikl)**2) *signus &
    611611            / (hSalSV(ikl) * gravit * us__SV(ikl) * 3.25)
    612612    endif
    613613
    614     if (qsalt_param == "bin") THEN
     614    IF (qsalt_param == "bin") THEN
    615615      qSalSV(ikl) = (us__SV(ikl) * us__SV(ikl) &
    616616            -usthSV(ikl) * usthSV(ikl))*signus &
    … …  
    662662  ! #BS       density_kotlyakov = .FALSE.  !C.Amory BS 2018
    663663  ! + ...     Fallen Snow Density, Adapted for Antarctica
    664         if (is_ok_density_kotlyakov) THEN
     664        IF (is_ok_density_kotlyakov) THEN
    665665            tt_tmp = TaT_SV(ikl)-TfSnow
    666666            !vv_tmp = VV10SV(ikl)
    … …  
    668668  ! + ...         [ A compromise between
    669669  ! + ...           Kotlyakov (1961) and Lenaerts (2012, JGR, Part1) ]
    670             if (tt_tmp>=-10) THEN
     670            IF (tt_tmp>=-10) THEN
    671671              BufPro   =  max( rosMin, &
    672672                    104. *sqrt( max( vv_tmp-6.0,0.0))) ! Kotlyakov (1961)
    … …  
    696696  !    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    697697
    698      if (BloMod) THEN
     698     IF (BloMod) THEN
    699699     Bros_N      = frsno
    700700     ro_new      = ro__SV(ikl,max(1,isnoSV(ikl)))
    … …  
    894894
    895895
    896   if (discret_xf.AND.klonv==1) THEN
    897    IF(isnoSV(1)>=1.or.NLaysv(1)>=1) THEN
     896  IF (discret_xf.AND.klonv==1) THEN
     897   IF(isnoSV(1)>=1.OR.NLaysv(1)>=1) THEN
    898898  ! +          **********
    899899     CALL SISVAT_zSn
    … …  
    904904    CALL SISVAT_zSn
    905905  ! +          **********
    906   endif
     906  ENDIF
    907907
    908908  ! +          **********
    … …  
    10801080
    10811081
    1082   if (iflag_temp_inlandsis == 0) THEN
     1082  IF (iflag_temp_inlandsis == 0) THEN
    10831083   CALL SISVAT_TSo
    10841084
    … …  
    12031203  ! Etienne: extrapolation from the two uppermost levels:
    12041204
    1205      if (isnoSV(ikl) >=2) THEN
     1205     IF (isnoSV(ikl) >=2) THEN
    12061206       zm1=-dzsnSV(ikl,isnoSV(ikl))/2.
    12071207       zm2=-(dzsnSV(ikl,isnoSV(ikl)) + dzsnSV(ikl,isnoSV(ikl)-1)/2.)
    1208      else if (isnoSV(ikl) == 1) THEN
     1208     ELSE IF (isnoSV(ikl) == 1) THEN
    12091209       zm1=-dzsnSV(ikl,isnoSV(ikl))/2.
    12101210       zm2=-(dzsnSV(ikl,isnoSV(ikl))+dz_dSV(0)/2.)
    … …  
    12341234  IF (SnoMod)                                                 THEN
    12351235
    1236   if (discret_xf .AND. klonv==1) THEN
     1236  IF (discret_xf .AND. klonv==1) THEN
    12371237  IF(isnoSV(1)>=1) THEN
    12381238  ! +          **********
    12391239  CALL SISVAT_GSn
    12401240  ! +          **********
    1241   endif
     1241  ENDIF
    12421242  else
    12431243  ! +          **********
    12441244    CALL SISVAT_GSn
    12451245  ! +          **********
    1246   endif
     1246  ENDIF
    12471247
    12481248
    … …  
    13341334      coefd = log(z03)-coefc*ta3
    13351335
    1336       if (TaT_SV(ikl) < ta1) THEN
     1336      IF (TaT_SV(ikl) < ta1) THEN
    13371337        Z0_obs = z01
    1338       else if (TaT_SV(ikl)>=ta1 .and. TaT_SV(ikl)<ta2) THEN
     1338      ELSE IF (TaT_SV(ikl)>=ta1 .AND. TaT_SV(ikl)<ta2) THEN
    13391339        Z0_obs = exp(coefa*TaT_SV(ikl) + coefb)
    1340       else if (TaT_SV(ikl)>=ta2 .and. TaT_SV(ikl)<ta3) THEN
     1340      ELSE IF (TaT_SV(ikl)>=ta2 .AND. TaT_SV(ikl)<ta3) THEN
    13411341        ! if st > 0, melting induce smooth surface
    13421342        Z0_obs = exp(coefc*TaT_SV(ikl) + coefd)
    … …  
    14711471
    14721472     IF(ro__SV(ikl,isnoSV(ikl))>50 &
    1473            .and.ro__SV(ikl,isnoSV(ikl))<roSdSV)THEN
     1473           .AND.ro__SV(ikl,isnoSV(ikl))<roSdSV)THEN
    14741474         Z0hnSV(ikl) = max(zero &
    14751475               , sign(unun,zzsnsv(ikl,isnoSV(ikl))-epsi)) &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_bsn.f90

    r5116 r5117  
    2020  ! +  =================
    2121
    22   use VARphy
    23   use VAR_SV
    24   use VARdSV
    25   use VARxSV
    26   use VARySV
     22  USE VARphy
     23  USE VAR_SV
     24  USE VARdSV
     25  USE VARxSV
     26  USE VARySV
    2727
    2828
    … …  
    7171
    7272    ! Densification of the uppermost snow layer if erosion:
    73     if((dzweqo-dzweqn)>0                    .and. &
    74           dzsnSV(ikl,isn)>0                    .and. &
     73    IF((dzweqo-dzweqn)>0                    .AND. &
     74          dzsnSV(ikl,isn)>0                    .AND. &
    7575          ro__SV(ikl,max(1,isnoSV(ikl)))<roBdSV) THEN
    7676    !characteristic time scale for drifting snow compaction set to 24h
    … …  
    8383    Fac         = max(0.,min(1.,Fac))
    8484
    85     if (ro__SV(ikl,max(1,isnoSV(ikl)))>roBdSV) THEN
     85    IF (ro__SV(ikl,max(1,isnoSV(ikl)))>roBdSV) THEN
    8686      densif=densif*Fac
    8787    endif
    … …  
    9393    endif
    9494
    95     IF(dzsnSV(ikl,isn)>0 .and.dzsnSV(ikl,isn)<0.0001)THEN
     95    IF(dzsnSV(ikl,isn)>0 .AND.dzsnSV(ikl,isn)<0.0001)THEN
    9696    dbs_SV(ikl) = dbs_SV(ikl)+ dzsnSV(ikl,isn)*ro__SV(ikl,isn)
    9797    dbs_Er(ikl) = dbs_Er(ikl)+ dzsnSV(ikl,isn)*ro__SV(ikl,isn)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_gsn.f90

    r5116 r5117  
    165165  ! +  ================
    166166
    167   use VARphy
    168   use VAR_SV
    169   use VARdSV
    170   use VAR0SV
    171   use VARxSV
    172   use VARtSV
     167  USE VARphy
     168  USE VAR_SV
     169  USE VARdSV
     170  USE VAR0SV
     171  USE VARxSV
     172  USE VARtSV
    173173
    174174
    … …  
    190190  ! +  ================
    191191
    192   logical :: vector                        !
     192  LOGICAL :: vector                        !
    193193  INTEGER :: ikl                           !
    194194  INTEGER :: isn   ,isnp                   !
    … …  
    806806    DO isn=iiceSV(ikl),isnoSV(ikl)
    807807      IF  (G1snSV(ikl,isn)>=0.)                               THEN
    808         IF(G1snSV(ikl,isn)<vsphe4.and.istoSV(ikl,isn)==0)  THEN
     808        IF(G1snSV(ikl,isn)<vsphe4.AND.istoSV(ikl,isn)==0)  THEN
    809809               istoSV(ikl,isn)=istdSV(1)
    810         ELSEIF(G1_dSV-G1snSV(ikl,isn)         <vsphe4.and. &
     810        ELSEIF(G1_dSV-G1snSV(ikl,isn)         <vsphe4.AND. &
    811811              etaSno(ikl,isn)/dzsnSV(ikl,isn)>vtelv1)       THEN
    812812          IF  (istoSV(ikl,isn)==0) &
    … …  
    858858            *(1-min(1   , abs(istoSV(ikl,isn)-istdSV(1))))
    859859  ! #wp     IF (G1snSV(ikl,isn).gt.0..AND.G1snSV(ikl,isn).lt.vsphe4
    860   ! #wp.                             .AND.istoSV(ikl,isn).eq.     0)
     860  ! #wp.                             .AND.istoSV(ikl,isn).EQ.     0)
    861861  ! #wp.    THEN
    862862  ! #wp       WRITE(6,*) ikl,isn,' G1,G2,hist,OK_Ang  ',
    … …  
    907907
    908908    isn=1
    909     if (dzsnSV(ikl,isn)>0.and.ro__SV(ikl,isn)>0) THEN
     909    IF (dzsnSV(ikl,isn)>0.AND.ro__SV(ikl,isn)>0) THEN
    910910    dzsnSV(ikl,isn) = dzsnSV(ikl,isn) +0.9999*(smb_old-smb_new) &
    911911          / ro__SV(ikl,isn)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_qsn.f90

    r5116 r5117  
    5555  ! +  ================
    5656
    57   use VARphy
    58   use VAR_SV
    59   use VARdSV
    60   use VAR0SV
    61   use VARxSV
    62   use VARySV
    63   use surface_data, ONLY: is_ok_slush,opt_runoff_ac
     57  USE VARphy
     58  USE VAR_SV
     59  USE VARdSV
     60  USE VAR0SV
     61  USE VARxSV
     62  USE VARySV
     63  USE surface_data, ONLY: is_ok_slush,opt_runoff_ac
    6464
    6565
    … …  
    358358
    359359  !XF
    360       IF(ro__SV(ikl,isn) >= roCdSV.and.ro__SV(ikl,1)<900) &
     360      IF(ro__SV(ikl,isn) >= roCdSV.AND.ro__SV(ikl,1)<900) &
    361361            PClose = min(0.50,PClose * &
    362362            (1.-(ro_ice-ro__SV(ikl,isn))/(ro_ice-roCdSV)))
    … …  
    369369      endif
    370370
    371       IF(drr_SV(ikl)    >0  .and.TsisSV(ikl,isn)>273.14) THEN
    372        if((ro__SV(ikl,isn)>900.and.ro__SV(ikl,isn)<920).or. &
     371      IF(drr_SV(ikl)    >0  .AND.TsisSV(ikl,isn)>273.14) THEN
     372       IF((ro__SV(ikl,isn)>900.AND.ro__SV(ikl,isn)<920).OR. &
    373373             ro__SV(ikl,isn)>950) THEN
    374374         dzsnSV(ikl,isn) = dzsnSV(ikl,isn)*ro__SV(ikl,isn)/ro_ice
    … …  
    378378      endif
    379379
    380       ! if (isn>1.and.isn<nsno     .and.
    381   !    .      ro__SV(ikl,isn-1)>900    .and.
    382   !    .      ro__SV(ikl,isn)  >roCdSV .and.
    383   !    .      ro__SV(ikl,isn)  <900    .and.
    384   !    .      TsisSV(ikl,isn)  >273.14 .and.
    385   !    .      TsisSV(ikl,isn+1)<273.15 .and.
     380      ! if (isn>1.AND.isn<nsno     .AND.
     381  !    .      ro__SV(ikl,isn-1)>900    .AND.
     382  !    .      ro__SV(ikl,isn)  >roCdSV .AND.
     383  !    .      ro__SV(ikl,isn)  <900    .AND.
     384  !    .      TsisSV(ikl,isn)  >273.14 .AND.
     385  !    .      TsisSV(ikl,isn+1)<273.15 .AND.
    386386  !    .      drr_SV(ikl)      >0)     THEN
    387387      !  TsisSV(ikl,isn)=273.14
    388388      !  PClose = 1
    389       ! endif
     389      ! END IF
    390390
    391391  !XF
    … …  
    526526
    527527     IF(zt<0.005+(TaT_SV(ikl)-TfSnow)/1000..and. &
    528            isnoSV(ikl)             >0         .and. &
    529            TaT_SV(ikl)             >=TfSnow   .and. &
     528           isnoSV(ikl)             >0         .AND. &
     529           TaT_SV(ikl)             >=TfSnow   .AND. &
    530530           istoSV(ikl,isnoSV(ikl)) >1       ) THEN
    531531      DO isn=1,isnoSV(ikl)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_qso.f90

    r5116 r5117  
    8585  ! +  ================
    8686
    87   use VARphy
    88   use VAR_SV
    89   use VARdSV
    90   use VAR0SV
    91   use VARxSV
    92   use VARySV
     87  USE VARphy
     88  USE VAR_SV
     89  USE VARdSV
     90  USE VAR0SV
     91  USE VARxSV
     92  USE VARySV
    9393
    9494
    … …  
    377377      ikp      = nkhy *  eta_SV(ikl,isl+1) / etadSV(ist)
    378378
    379       IF(ikm<0.or.ik0<0.or.ikp<0)THEN
     379      IF(ikm<0.OR.ik0<0.OR.ikp<0)THEN
    380380       print *,"CRASH1 in sisvat_qso.f on pixel (i,j,n)", &
    381381             ii__SV(ikl),jj__SV(ikl),nn__SV(ikl)
    … …  
    423423      ikp      = nkhy *  eta_SV(ikl,isl+1) / etadSV(ist)
    424424
    425       IF(ik0<0.or.ikp<0)THEN
     425      IF(ik0<0.OR.ikp<0)THEN
    426426       print *,"CRASH2 in sisvat_qso.f on pixel (i,j,n)", &
    427427             ii__SV(ikl),jj__SV(ikl),nn__SV(ikl)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_sno_albedo.f90

    r5116 r5117  
    6161
    6262
    63   use VARphy
    64   use VAR_SV
    65   use VARdSV
    66   use VARxSV
    67   use VARySV
    68   use VARtSV
     63  USE VARphy
     64  USE VAR_SV
     65  USE VARdSV
     66  USE VARxSV
     67  USE VARySV
     68  USE VARtSV
    6969  USE surface_data, ONLY: iflag_albcalc,correc_alb
    7070
    … …  
    464464    ! prescription for each time step with NEMO values
    465465
    466   ! #AO      if (LSmask(ikl) .eq. 0 .and. coupling_ao .eq. .TRUE.) THEN
    467   ! #AO       if (AOmask(ikl) .eq. 0) THEN
     466  ! #AO      if (LSmask(ikl) .EQ. 0 .AND. coupling_ao .EQ. .TRUE.) THEN
     467  ! #AO       if (AOmask(ikl) .EQ. 0) THEN
    468468  ! #AO       albisv(ikl) =  (1.-AOmask(ikl))* albAOsisv(ikl)
    469469  ! #AO.                    +(AOmask(ikl)*albisv(ikl))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.f90

    r5116 r5117  
    7171  ! +  ================
    7272
    73   use VARphy
    74   use VAR_SV
    75   use VARdSV
    76   use VARxSV
    77   use VARySV
    78   use VARtSV
    79   use VAR0SV
     73  USE VARphy
     74  USE VAR_SV
     75  USE VARdSV
     76  USE VARxSV
     77  USE VARySV
     78  USE VARtSV
     79  USE VAR0SV
    8080
    8181
    … …  
    615615      TsisSV(ikl,isl)   = Aux__Q(ikl,isl)  *TsisSV(ikl,isl+1) &
    616616            +TsisSV(ikl,isl)
    617       IF(isl==0.and.isnoSV(ikl)==0) THEN
     617      IF(isl==0.AND.isnoSV(ikl)==0) THEN
    618618       TsisSV(ikl,isl)  = min(TaT_SV(ikl)+30,TsisSV(ikl,isl))
    619619       TsisSV(ikl,isl)  = max(TaT_SV(ikl)-30,TsisSV(ikl,isl))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_weq.f90

    r5116 r5117  
    2222  ! +  ================
    2323
    24   use VARphy
    25   use VAR_SV
    26   use VARxSV
     24  USE VARphy
     25  USE VAR_SV
     26  USE VARxSV
    2727
    2828  IMPLICIT NONE
    … …  
    3333  INTEGER :: istart
    3434
    35   logical :: logWEq
     35  LOGICAL :: logWEq
    3636  common/SISVAT_wEq_L/logWEq
    3737
    … …  
    8484  ! +  ======
    8585
    86   !!      IF (istart.eq.1)                                              THEN
     86  !!      IF (istart.EQ.1)                                              THEN
    8787  !!        WRITE(45,45)dahost,i___SV(lwriSV(1)),j___SV(lwriSV(1)),
    8888  !!     .              n___SV(lwriSV(1))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zag.f90

    r5116 r5117  
    5151  ! +  ================
    5252
    53   use VARphy
    54   use VAR_SV
    55   use VARdSV
    56   use VAR0SV
    57   use VARxSV
     53  USE VARphy
     54  USE VAR_SV
     55  USE VARdSV
     56  USE VAR0SV
     57  USE VARxSV
    5858
    5959  IMPLICIT NONE

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zcr.f90

    r5116 r5117  
    3434  ! +  ================
    3535  !
    36   use VARphy
    37   use VAR_SV
    38   use VARdSV
    39   use VAR0SV
    40   use VARxSV
    41   use VARySV
    42   use VARtSV
     36  USE VARphy
     37  USE VAR_SV
     38  USE VARdSV
     39  USE VAR0SV
     40  USE VARxSV
     41  USE VARySV
     42  USE VARtSV
    4343
    4444  IMPLICIT NONE

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zsn.f90

    r5116 r5117  
    4646
    4747
    48   use VARphy
    49   use VAR_SV
    50   use VARdSV
    51   use VAR0SV
    52   use VARxSV
    53   use VARySV
    54   use surface_data, ONLY: ok_zsn_ii
     48  USE VARphy
     49  USE VAR_SV
     50  USE VARdSV
     51  USE VAR0SV
     52  USE VARxSV
     53  USE VARySV
     54  USE surface_data, ONLY: ok_zsn_ii
    5555
    5656  IMPLICIT NONE
    … …  
    598598  ! + minimum uppermost layer thickness to guarantee a correct reproduction of the snow
    599599  ! + atmosphere coupling
    600     IF(dzsnSV(ikl,max(1,isnoSV(ikl)-0))>0.02 .or. & ! surface layers> 2-5-10
    601           dzsnSV(ikl,max(1,isnoSV(ikl)-1))>0.05 .or. & ! XF 04/07/2019
    602           dzsnSV(ikl,max(1,isnoSV(ikl)-2))>0.10 .or. &
     600    IF(dzsnSV(ikl,max(1,isnoSV(ikl)-0))>0.02 .OR. & ! surface layers> 2-5-10
     601          dzsnSV(ikl,max(1,isnoSV(ikl)-1))>0.05 .OR. & ! XF 04/07/2019
     602          dzsnSV(ikl,max(1,isnoSV(ikl)-2))>0.10 .OR. &
    603603          dzsnSV(ikl,max(1,isnoSV(ikl)-3))>0.30 )THEN
    604604      Agrege(ikl) = min(1, &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90

    r5116 r5117  
    158158        ! CA initialization
    159159        ! dz_profil_15 : 1 m in 15 layers [m]
    160         real, parameter :: dz_profil_15(15) = (/0.005, 0.01, 0.015, 0.02, 0.03, 0.04, 0.05, &
     160        REAL, parameter :: dz_profil_15(15) = (/0.005, 0.01, 0.015, 0.02, 0.03, 0.04, 0.05, &
    161161                                                0.06, 0.07, 0.08, 0.09, 0.1, 0.12, 0.14, 0.17/)
    162162        ! mean_temp : mean annual surface temperature [K]
    163         real, dimension(klon) :: mean_temp
     163        REAL, DIMENSION(klon) :: mean_temp
    164164        ! mean_dens : mean surface density [kg/m3]
    165         real, dimension(klon) :: mean_dens
     165        REAL, DIMENSION(klon) :: mean_dens
    166166        ! lat_scale : temperature lapse rate against latitude [K degree-1]
    167167        REAL :: lat_scale
    … …  
    183183        ! constants for density profile
    184184        ! C0, C1 : constant, 0.07 for z <= 550 kg m-3
    185         real, parameter :: C0 = 0.07
    186         real, parameter :: C1 = 0.03
     185        REAL, parameter :: C0 = 0.07
     186        REAL, parameter :: C1 = 0.03
    187187        ! rho_i : ice density [kg m-3]
    188         real, parameter :: rho_ice = 917.
     188        REAL, parameter :: rho_ice = 917.
    189189        ! E_c : activation energy [J mol-1]
    190         real, parameter :: E_c = 60000.
     190        REAL, parameter :: E_c = 60000.
    191191        ! E_g : activation energy [J mol-1]
    192         real, parameter :: E_g = 42400.
     192        REAL, parameter :: E_g = 42400.
    193193        ! R : gas constant [J mol-1 K-1]
    194         real, parameter :: R = 8.3144621
     194        REAL, parameter :: R = 8.3144621
    195195
    196196     
    … …  
    368368                ! with a moist-adiabatic lapse rate of 5 °C km-1 everywhere except for Antarctica,
    369369                ! for Antarctica, a dry-adiabatic lapse rate of 9.8 °C km-1 is assumed.
    370                 if (lati > 60.) THEN
     370                IF (lati > 60.) THEN
    371371                    ! CA todo : add longitude bounds
    372372                    ! Greenland mean temperature : function of altitude and latitude
    … …  
    379379                    ! surface density: Fausto et al. 2018, https://doi.org/10.3389/feart.2018.00051
    380380                    mean_dens(ikl) = 315.
    381                 else if (lati < -60.) THEN
     381                ELSE IF (lati < -60.) THEN
    382382                    ! Antarctica mean temperature : function of altitude and latitude
    383383                    ! for altitudes 0. to 500. m, lat_scale varies from 1.3 to 0.6 °C °lat-1
    … …  
    463463                    rho0 = exp(E0 * depth) / (rho_ice / mean_dens(ikl) - 1 + exp(E0 * depth)) * rho_ice
    464464                    rho1 = exp(E1 * depth) / (rho_ice / mean_dens(ikl) - 1 + exp(E1 * depth)) * rho_ice
    465                     if (depth <= z550) THEN
     465                    IF (depth <= z550) THEN
    466466                        ro__SV(ikl, isl) = exp(E0 * depth) / (rho_ice / mean_dens(ikl) - 1 + exp(E0 * depth)) * rho_ice
    467467                    else
    … …  
    553553            ! => Upper bound for eroded snow mass
    554554            !        uss_SV(ikl) = SLussl(i,j,n) ! u*qs* (only for Tv in sisvatesbl.f)
    555             ! #BS  IF(dsn_SV(ikl)>eps12.and.erprev(i,j,n).gt.eps9) THEN
     555            ! #BS  IF(dsn_SV(ikl)>eps12.AND.erprev(i,j,n).gt.eps9) THEN
    556556            ! #BS    dsnbSV(ikl) =1.0-min(qsHY(i,j,kB)     !BS neglib. at kb ~100 magl)
    557557            ! #BS.                        /max(qshy(i,j,mz),eps9),unun)
    … …  
    560560            ! #BS  else
    561561            ! #BS    dsnbSV(ikl) = 0.
    562             ! #BS  endif
     562            ! #BS  ENDIF
    563563            !      dsnbSV is the drift fraction of deposited snow updated in sisvat.f
    564564            !      will be used for characterizing the Buffer Layer
    … …  
    586586            ! +--Energy Fluxes                                          (INPUT/OUTPUT)
    587587            ! +  ^^^^^^^^^^^^^                                           ^^^^^^^^^^^^
    588             !IF (.not.firstcall) THEN
     588            !IF (.NOT.firstcall) THEN
    589589            Tsrfsv(ikl)  = tsurf(ikl)                     !hj 12 03 2010
    590590            cld_SV(ikl) = cloudf(ikl)                    ! Cloudiness
    … …  
    932932        !c #kd END IF
    933933
    934         !      IF (nsol.ne.4)                                              THEN
     934        !      IF (nsol.NE.4)                                              THEN
    935935        !        DO isl= 0,-nsol,-1
    936936        !          misl_2 =     -mod(isl,2)
Note: See TracChangeset for help on using the changeset viewer.