Changeset 1048


Ignore:
Timestamp:
Sep 23, 2013, 7:02:50 PM (11 years ago)
Author:
slebonnois
Message:

SL: update pour divers details titan + quelques modifs arch et makelmdz

Location:
trunk
Files:
4 added
11 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/makelmdz

    r1019 r1048  
    595595########################################################################
    596596cd $LMDGCM
    597 find libf -name '*.[Fh]' -print > tmp77
    598 find libf -name '*.[Fh]' -exec egrep -i " *use *ioipsl" {} \; -print > tmp90
    599 find libf -name '*.[Fh90]' -print >> tmp90
     597find libf -name '*.[Fh]' -print | sort > tmp77
     598find libf -name '*.[Fh]' -exec egrep -i " *use *ioipsl" {} \; -print | sort > tmp90
     599find libf -name '*.[Fh90]' -print | sort >> tmp90
    600600######### CAS PARTICULIER CHIMIE TITAN
    601601if [[ "$physique" == "titan" ]]
    602602then
    603   find libf -name '*.[ch]' -print > tmpC
     603  find libf -name '*.[ch]' -print | sort > tmpC
    604604fi
    605605#########

  • trunk/LMDZ.TITAN/libf/phytitan/YOMCST.h

    r102 r1048  
    1919      REAL R,RMD,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
    2020      REAL RKAPPA,RETV
    21 ! ADAPTATION GCM POUR CP(T)
    22       real cpdet
    23       external cpdet
    2421! A1.5,6 Thermodynamic liquid,solid phases
    2522      REAL RCW,RCS

  • trunk/LMDZ.TITAN/libf/phytitan/ajsec.F

    r102 r1048  
    77
    88      use dimphy
     9      use cpdet_mod, only: t2tpot, tpot2t
    910      IMPLICIT none
    1011c======================================================================

  • trunk/LMDZ.TITAN/libf/phytitan/clesphys.h

    r495 r1048  
    3838
    3939       COMMON/clesphys_r/                                               &
    40      &     ecritphy, solaire, z0, lmixmin, ksta, inertie,               &
     40     &     ecritphy, solaire, z0, lmixmin, ksta, inertie, emis,         &
    4141     &     tx,tcorrect,p_prodaer,xnuf
    4242

  • trunk/LMDZ.TITAN/libf/phytitan/clmain.F

    r495 r1048  
    3535      USE interface_surf
    3636      use dimphy
     37      use cpdet_mod, only: t2tpot
    3738      IMPLICIT none
    3839c======================================================================
    … …  
    342343         y_cd_m(1:knon) = ycoefh(1:knon,1)
    343344         endif
     345
    344346         call ustarhb(knon,yu,yv,y_cd_m, yustar)
    345347
    … …  
    478480      USE interface_surf
    479481      use dimphy
     482      use cpdet_mod, only: t2tpot,tpot2t,cpdet
    480483
    481484      IMPLICIT none
    … …  
    886889
    887890      use dimphy
     891      use cpdet_mod, only: cpdet,t2tpot
    888892      IMPLICIT none
    889893c======================================================================
    … …  
    11691173
    11701174      use dimphy
     1175      use cpdet_mod, only: cpdet
    11711176      IMPLICIT none
    11721177c======================================================================

  • trunk/LMDZ.TITAN/libf/phytitan/comcstfi.h

    r3 r1048  
    66
    77      REAL pi,rad,g,r,cpp,rcp,dtphys,daysec,mugaz,omeg
    8       real cpdet
    9       external cpdet
    108
    119c-----------------------------------------------------------------------

  • trunk/LMDZ.TITAN/libf/phytitan/diagphy.F

    r102 r1048  
    209209 
    210210      use dimphy
     211      use cpdet_mod, only: cpdet
    211212      IMPLICIT NONE
    212213C

  • trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90

    r495 r1048  
    4545      & tsol_rad, tsurf_new, alb_new)
    4646
     47      use write_field_phy
     48      use cpdet_mod, only: cpdet
     49
     50      IMPLICIT none
    4751
    4852! Cette routine sert d'aiguillage entre l'atmosphere et la surface en general
    … …  
    141145
    142146  if (check) write(*,*) 'Entree ', modname
    143  
     147
    144148! Initialisations diverses
    145149!
    … …  
    149153
    150154! ADAPTATION GCM POUR CP(T)
    151        zcp(1:klon) = cpdet(tsurf(1:klon))
     155       do ii=1,klon
     156         zcp(ii)=cpdet(tsurf(ii))
     157       enddo
    152158
    153159       IF (soil_model) THEN
    154160           CALL soil(dtime, knon, tsurf, tsoil,soilcap, soilflux)
    155161           cal(1:knon) = zcp(1:knon) / soilcap(1:knon)
     162! for tests:
     163!  call writefield_phy('interfsurf_hq_zcp',zcp,1)
     164!  call writefield_phy('interfsurf_hq_cal',cal,1)
     165!  call writefield_phy('interfsurf_hq_soilcap',soilcap,1)
    156166!       print*,"DIAGNOSTIC SOIL"
    157167!       print*,"soilcap=",soilcap
    … …  
    186196     & tsurf_new, fluxsens, dflux_s)
    187197
     198  use write_field_phy
     199  use cpdet_mod, only: t2tpot, tpot2t
     200
     201  IMPLICIT none
     202
    188203! Cette routine calcule les fluxs en h a l'interface et eventuellement
    189204! une temperature de surface (au cas ou ok_veget = false)
    … …  
    317332  ENDDO
    318333
     334! for tests: write output fields...
     335!  call writefield_phy('calcul_fluxs_d_ts',d_ts,1)
     336!  call writefield_phy('calcul_fluxs_fluxsens',fluxsens,1)
     337!  call writefield_phy('calcul_fluxs_dflux_s',dflux_s,1)
     338
    319339  END SUBROUTINE calcul_fluxs
    320340!

  • trunk/LMDZ.TITAN/libf/phytitan/physiq.F

    r953 r1048  
    6161      use dimphy
    6262      USE comgeomphy
     63      use cpdet_mod, only: cpdet, t2tpot
    6364      IMPLICIT none
    6465c======================================================================
    … …  
    279280      REAL zday
    280281      REAL zls,zlsdeg,zlsm1
     282      save zlsm1
    281283c
    282284      INTEGER i, k, iq, ig, j, ll, l

  • trunk/LMDZ.TITAN/libf/phytitan/rcm1d.F

    r904 r1048  
    1 
    21      PROGRAM rcm1d
    32     
    … …  
    54      use control_mod
    65      use comgeomphy
     6      use cpdet_mod, only: ini_cpdet
    77      IMPLICIT NONE
    88
    … …  
    474474#include "../dyn3d/disvert_noterre.F"
    475475#include "../dyn3d/abort_gcm.F"
    476 #include "../dyn3d/dump2d.F"
    477 #include "../dyn3d/cpdet.F"
     476!#include "../dyn3d/dump2d.F"
    478477
    479478c***********************************************************************

  • trunk/LMDZ.VENUS/libf/phyvenus/grid_noro.F

    r808 r1048  
    5252C=======================================================================
    5353
    54       IMPLICIT INTEGER (I,J)
    55       IMPLICIT REAL(X,Z)
     54      IMPLICIT none
    5655     
    5756#include "dimensions.h"
     57#include "YOMCST.h"
    5858
    5959      INTEGER imdep, jmdep
    … …  
    6363      INTEGER imar, jmar
    6464c parametres lies au fichier d entree... A documenter...
    65       parameter(iext=216, epsfra = 1.e-5)
     65      integer iext
     66      parameter(iext=216)
    6667      REAL xusn(imdep+2*iext),yusn(jmdep+2)
    6768      REAL zusn(imdep+2*iext,jmdep+2)
    6869 
     70c local var
     71      real zdeltax,zdeltay,zlenx,zleny,xincr
     72      real zbordnor,zbordsud,zbordest,zbordoue,weighx,weighy
     73      real zllmmea,zllmstd,zllmsig,zllmgam,zllmpic,zllmval,zllmthe
     74      real zminthe,xk,xl,xm,xp,xq,xw
     75      real zmeanor,zmeasud,zstdnor,zstdsud,zsignor,zsigsud
     76      real zweinor,zweisud,zpicnor,zpicsud,zvalnor,zvalsud
     77      integer i,j,ii,jj
     78
    6979C INTERMEDIATE FIELDS  (CORRELATIONS OF OROGRAPHY GRADIENT)
    7080
    … …  
    8797c
    8898      print *,' parametres de l orographie a l echelle sous maille'
    89       xpi=acos(-1.)
    90       rad    = 6 371 229.
    91       zdeltay=2.*xpi/REAL(jmdep)*rad
     99
     100      zdeltay=2.*RPI/REAL(jmdep)*RA
    92101c
    93102c  quelques tests de dimensions:
    … …  
    123132      DO i=1,iext
    124133        zusn(i,j+1)=zdata(imdep-iext+i,j)
    125         xusn(i)=xdata(imdep-iext+i)-2.*xpi
     134        xusn(i)=xdata(imdep-iext+i)-2.*RPI
    126135        zusn(imdep+iext+i,j+1)=zdata(i,j)
    127         xusn(imdep+iext+i)=xdata(i)+2.*xpi
     136        xusn(imdep+iext+i)=xdata(i)+2.*RPI
    128137      ENDDO
    129138      ENDDO
    … …  
    197206c  SUMMATION OVER GRIDPOINT AREA
    198207c
    199       zleny=xpi/REAL(jmdep)*rad
    200       xincr=xpi/2./REAL(jmdep)
     208      zleny=RPI/REAL(jmdep)*RA
     209      xincr=RPI/2./REAL(jmdep)
    201210       DO ii = 1, imar+1
    202211       DO jj = 1, jmar
    … …  
    208217            zlenx=zleny*cos(yusn(j))
    209218            zdeltax=zdeltay*cos(yusn(j))
    210             zbordnor=(c(jj)-yusn(j)+xincr)*rad
    211             zbordsud=(yusn(j)-d(jj)+xincr)*rad
     219            zbordnor=(c(jj)-yusn(j)+xincr)*RA
     220            zbordsud=(yusn(j)-d(jj)+xincr)*RA
    212221            weighy=AMAX1(0.,
    213222     *             amin1(zbordnor,zbordsud,zleny))
    214223         IF(weighy.ne.0)THEN
    215224         DO i = 2, imdep+2*iext-1
    216             zbordest=(xusn(i)-a(ii)+xincr)*rad*cos(yusn(j))
    217             zbordoue=(b(ii)+xincr-xusn(i))*rad*cos(yusn(j))
     225            zbordest=(xusn(i)-a(ii)+xincr)*RA*cos(yusn(j))
     226            zbordoue=(b(ii)+xincr-xusn(i))*RA*cos(yusn(j))
    218227            weighx=AMAX1(0.,
    219228     *             amin1(zbordest,zbordoue,zlenx))
Note: See TracChangeset for help on using the changeset viewer.