Changeset 4192


Ignore:
Timestamp:
Dec 20, 2019, 12:56:43 AM (5 years ago)
Author:
dubos
Message:

simple_physics : cleanup astronomy

Location:
dynamico_lmdz/simple_physics/phyparam
Files:
2 deleted
4 edited

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/simple_physics/phyparam/param/iniphyparam.F

    r4191 r4192  
    6464      INTEGER ig,ierr,offset
    6565 
    66       EXTERNAL inifrict,orbite
     66      EXTERNAL inifrict
    6767 
    6868      print*,'INIPHYPARAM'

  • dynamico_lmdz/simple_physics/phyparam/param/phyparam.F

    r4191 r4192  
    164164
    165165      EXTERNAL convadj
    166       EXTERNAL orbite,mucorr
     166      EXTERNAL mucorr
    167167      EXTERNAL ismin,ismax
    168168
    … …  
    482482         ENDDO
    483483
    484          CALL zerophys(ngrid*nlayer,zdum1)
    485          CALL zerophys(ngrid*nlayer,zdum2)
    486 c        CALL zerophys(ngrid*nlayer,zdum3)
     484         zdum1(:,:)=0.
     485         zdum2(:,:)=0.
     486
    487487         do l=1,nlayer
    488488            do ig=1,ngrid
    … …  
    546546            ENDDO
    547547         ENDDO
    548          CALL zerophys(ngrid*nlayer,zdufr)
    549          CALL zerophys(ngrid*nlayer,zdvfr)
    550          CALL zerophys(ngrid*nlayer,zdhfr)
     548
     549         zdufr(:,:)=0.
     550         zdvfr(:,:)=0.
     551         zdhfr(:,:)=0.
     552
    551553         CALL convadj(ngrid,nlayer,ptimestep,
    552554     $                pplay,pplev,zpopsk,

  • dynamico_lmdz/simple_physics/phyparam/param/sw.F

    r4183 r4192  
    204204
    205205      IF (ldiurn) THEN
    206          CALL zerophys(ngrid,fsrfvis)
     206         fsrfvis(:)=0.
    207207         CALL monscatter(ncount,fsrfvis,index,zflux)
    208          CALL zerophys(ngrid*nlayer,dtsw)
     208         dtsw(:,:)=0.
    209209         DO l=1,nlayer
    210210            CALL monscatter(ncount,dtsw(1,l),index,zdtsw(1,l))

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

    r4191 r4192  
    77  REAL, PARAMETER :: unitastr=149.597927, & ! millions of km
    88       pi=2.*ASIN(1.)
     9
    910CONTAINS
    1011 
    … …  
    166167  END SUBROUTINE iniorbit
    167168
     169  PURE SUBROUTINE orbite(pls,pdist_sol,pdecli)
     170    !=======================================================================
     171    !
     172    !   Objet:
     173    !   ------
     174    !
     175    !   Distance from sun and declimation as a function of the solar
     176    !   longitude Ls
     177    !
     178    !   Arguments:
     179    !   ----------
     180    !
     181    !   Input:
     182    !   ------
     183    !   pls          Ls
     184    !
     185    !   Output:
     186    !   -------
     187    !   pdist_sol     Distance Sun-Planet in UA
     188    !   pdecli        declinaison ( en radians )
     189    !
     190    !=======================================================================
     191    !-----------------------------------------------------------------------
     192    !   Declarations:
     193    !   -------------
     194   
     195    ! arguments:
     196    ! ----------
     197
     198    REAL, INTENT(IN) :: pls
     199    REAL, INTENT(OUT) :: pdist_sol,pdecli
     200   
     201    !-----------------------------------------------------------------------
     202   
     203    ! Distance Sun-Planet
     204   
     205    pdist_sol=p_elips/(1.+e_elips*cos(pls+timeperi))
     206   
     207    ! Solar declination
     208   
     209    pdecli= asin (sin(pls)*sin(obliquit*pi/180.))
     210   
     211  END SUBROUTINE orbite
     212
    168213END MODULE astronomy
    169214
Note: See TracChangeset for help on using the changeset viewer.