Changeset 4233 for dynamico_lmdz


Ignore:
Timestamp:
Jan 18, 2020, 5:59:32 PM (5 years ago)
Author:
dubos
Message:

simple_physics : enforce F2003 strictly

Location:
dynamico_lmdz/simple_physics
Files:
2 added
10 edited

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/simple_physics/bash/beautify.sh

    r4229 r4233  
    33    # indent source file using emacs in batch mode
    44    # following standard settings of emacs f90-mode
    5     emacs -batch $x --eval '(f90-mode)' --eval '(indent-region (point-min) (point-max) nil)' --eval '(delete-trailing-whitespace)' -f save-buffer
     5    emacs -batch $x --eval '(f90-mode)' --eval '(indent-region (point-min) (point-max) nil)' --eval '(delete-trailing-whitespace)' -f save-buffer 2>/dev/null
    66    # remove blank lines at end of file
    77    sed -i '' -e :a -e '/^\n*$/{$d;N;};/\n$/ba' $x

  • dynamico_lmdz/simple_physics/phyparam/Makefile

    r4229 r4233  
    2222clean :
    2323        rm -f obj/* include/* lib/* xml/*
     24
     25nice:
    2426        ../bash/beautify.sh physics/*.F90
     27
    2528%.so : $(OBJECTS)
    2629        $(F90) -shared $^ -o $@

  • dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

    r4229 r4233  
    1414
    1515  SUBROUTINE solarlong(pday,psollong)
    16     REAL, INTENT(IN) :: pday               ! jour de l'annee (le jour 0 correspondant a l'equinoxe)
     16    REAL, INTENT(IN) :: pday               ! jour de l annee (le jour 0 correspondant a l equinoxe)
    1717    REAL, INTENT(OUT) :: psollong          ! solar longitude
    1818    LOGICAL, PARAMETER ::  lwrite=.TRUE.
    … …  
    2424
    2525    !--------------------------------------------------------
    26     ! calcul de l'angle polaire et de la distance au soleil :
     26    ! calcul de l angle polaire et de la distance au soleil :
    2727    ! -------------------------------------------------------
    2828
    29     !  calcul de l'zanomalie moyenne
     29    !  calcul de l zanomalie moyenne
    3030
    3131    zz=(pday-peri_day)/year_day
    … …  
    3333    xref=abs(zanom)
    3434
    35     !  resolution de l'equation horaire  zx0 - e * sin (zx0) = xref
     35    !  resolution de l equation horaire  zx0 - e * sin (zx0) = xref
    3636    !  methode de Newton
    3737
    … …  
    7474    !   ------
    7575    !    Initialisation du sous programme orbite qui calcule
    76     !    a une date donnee de l'annee de duree year_day commencant
    77     !    a l'equinoxe de printemps et dont le perihelie se situe
     76    !    a une date donnee de l annee de duree year_day commencant
     77    !    a l equinoxe de printemps et dont le perihelie se situe
    7878    !    a la date peri_day, la distance au soleil et la declinaison.
    7979    !
    … …  
    8181    !   ----------
    8282    !   - initialise certaines variables de ce module
    83     !   - Doit etre appele avant d'utiliser orbite.
     83    !   - Doit etre appele avant d utiliser orbite.
    8484    !
    8585    !   Arguments:
    … …  
    8888    !   Input:
    8989    !   ------
    90     !   aphelie       \   aphelie et perihelie de l'orbite
     90    !   aphelie       \   aphelie et perihelie de l orbite
    9191    !   periheli      /   en millions de kilometres.
    9292    !
    … …  
    113113
    114114    !-----------------------------------------------------------------------
    115     ! calcul de l'angle polaire et de la distance au soleil :
     115    ! calcul de l angle polaire et de la distance au soleil :
    116116    ! -------------------------------------------------------
    117117
    118     !  calcul de l'zanomalie moyenne
     118    !  calcul de l zanomalie moyenne
    119119
    120120    zz=(year_day-peri_day)/year_day
    … …  
    123123    WRITELOG(*,*) 'zanom  ',zanom
    124124
    125     !  resolution de l'equation horaire  zx0 - e * sin (zx0) = zxref
     125    !  resolution de l equation horaire  zx0 - e * sin (zx0) = zxref
    126126    !  methode de Newton
    127127

  • dynamico_lmdz/simple_physics/phyparam/physics/logging.F90

    r4230 r4233  
    9898    INTEGER, INTENT(IN) :: index ! index of gridpoint
    9999    CHARACTER(*), INTENT(OUT) :: line
    100     line=''
     100    line=' '
    101101  END SUBROUTINE default_log_gridpoint
    102102

  • dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90

    r4229 r4233  
    2222    !=======================================================================
    2323    !
    24     !   calcul de l'evolution de la temperature sous l'effet du rayonnement
     24    !   calcul de l evolution de la temperature sous l effet du rayonnement
    2525    !   infra-rouge.
    2626    !   Pour simplifier, les transmissions sont precalculees et ne
    27     !   dependent que de l'altitude.
     27    !   dependent que de l altitude.
    2828    !
    2929    !   arguments:
    … …  
    7676
    7777    !-----------------------------------------------------------------------
    78     !   2. calcul des quantites d'absorbants:
     78    !   2. calcul des quantites d absorbants:
    7979    !   -------------------------------------
    8080

  • dynamico_lmdz/simple_physics/phyparam/physics/radiative_sw.F90

    r4229 r4233  
    4242    !
    4343    !   Rayonnement solaire en atmosphere non diffusante avec un
    44     !   coefficient d'absoprption gris.
     44    !   coefficient d absoprption gris.
    4545    !
    4646    !=======================================================================
    … …  
    113113    tau0=-.5*log(coefvis)
    114114
    115     ! calcul de la partie homogene de l'opacite
     115    ! calcul de la partie homogene de l opacite
    116116    tau0=tau0/ps_rad
    117117    DO l=1,nlayer+1
    … …  
    122122
    123123    !-----------------------------------------------------------------------
    124     !   2. calcul de la transmission depuis le sommet de l'atmosphere:
     124    !   2. calcul de la transmission depuis le sommet de l atmosphere:
    125125    !   -----------------------------------------------------------
    126126
    … …  
    202202
    203203    !-----------------------------------------------------------------------
    204     !   6.ajout a l'echauffement de la contribution du ray. sol. reflechit:
     204    !   6.ajout a l echauffement de la contribution du ray. sol. reflechit:
    205205    !   -------------------------------------------------------------------
    206206

  • dynamico_lmdz/simple_physics/phyparam/physics/read_param_mod.F90

    r4230 r4233  
    9696    LOG_INFO('read_param')
    9797  END SUBROUTINE read_paramb
    98  
     98
    9999END MODULE read_param_mod

  • dynamico_lmdz/simple_physics/phyparam/physics/solar.F90

    r4229 r4233  
    8989          pfract(jl)=1.
    9090       ELSE
    91           pmu0(jl)=0.
     91          pmu0(jl)=0.
    9292          pfract(jl)=0.
    9393       ENDIF
    … …  
    102102    !=============================================================
    103103    ! Auteur : O. Boucher (LMD/CNRS)
    104     !          d'apres les routines zenith et angle de Z.X. Li
    105     ! Objet  : calculer les valeurs moyennes du cos de l'angle zenithal
    106     !          et l'ensoleillement moyen entre gmtime1 et gmtime2
     104    !          d apres les routines zenith et angle de Z.X. Li
     105    ! Objet  : calculer les valeurs moyennes du cos de l angle zenithal
     106    !          et l ensoleillement moyen entre gmtime1 et gmtime2
    107107    !          connaissant la declinaison, la latitude et la longitude.
    108108    ! Rque   : Different de la routine angle en ce sens que zenang
    … …  
    114114    !===============================================================
    115115    ! longi----INPUT : la longitude vraie de la terre dans son plan
    116     !                  solaire a partir de l'equinoxe de printemps (degre)
     116    !                  solaire a partir de l equinoxe de printemps (degre)
    117117    ! gmtime---INPUT : temps universel en fraction de jour
    118118    ! pdtrad---INPUT : pas de temps du rayonnement (secondes)
    … …  
    133133    ! omega1, omega2 : temps 1 et 2 exprime en radian avec 0 a midi.
    134134    ! omega : heure en radian du coucher de soleil
    135     ! -omega est donc l'heure en radian de lever du soleil
     135    ! -omega est donc l heure en radian de lever du soleil
    136136    real omegadeb, omegafin
    137137    real zfrac1, zfrac2, z1_mu, z2_mu
    … …  
    275275    !         pdeclin          solar declinaison
    276276    !         plat(npts)        latitude
    277     !         phaut            hauteur typique de l'atmosphere
     277    !         phaut            hauteur typique de l atmosphere
    278278    !         prad             rayon planetaire
    279279    !

  • dynamico_lmdz/simple_physics/phyparam/physics/surface.F90

    r4229 r4233  
    2828  SUBROUTINE init_soil(ngrid,nsoil)
    2929    INTEGER, INTENT(IN) :: ngrid, nsoil
    30     REAL min_period,dalph_soil
    31     REAL fz,rk,fz1,rk1,rk2
     30    REAL :: min_period,dalph_soil, rk,fz1,rk1,rk2
    3231    INTEGER :: jk
    33 
    34     ! this is a function definition
    35     fz(rk)=fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
    3632
    3733    !-----------------------------------------------------------------------
    … …  
    7167    ENDDO
    7268    LOG_INFO('init_soil')
     69
     70  CONTAINS
     71
     72    FUNCTION fz(rk) RESULT(val)
     73      REAL :: val, rk
     74      val = fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
     75    END FUNCTION fz
    7376
    7477  END SUBROUTINE init_soil

  • dynamico_lmdz/simple_physics/phyparam/physics/turbulence.F90

    r4229 r4233  
    388388
    389389    !-----------------------------------------------------------------------
    390     !   calcul le l'evolution de la temperature du sol:
     390    !   calcul le l evolution de la temperature du sol:
    391391    !   -----------------------------------------------
    392392
Note: See TracChangeset for help on using the changeset viewer.