Changeset 50


Ignore:
Timestamp:
Feb 15, 2000, 2:01:27 PM (24 years ago)
Author:
lmdz
Message:

Calcul de valeurs uniques aux poles (suppression du test alphax)
LF

Location:
LMDZ.3.3/trunk/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/dyn3d/etat0_netcdf.F

    r24 r50  
    284284      print *,' ALPHAX ',alphax
    285285
    286       IF( alphax.NE.0. )   THEN
    287        DO  l = 1, llm
    288          DO  i    = 1, iim
     286      DO  l = 1, llm
     287        DO  i    = 1, iim
    289288          xppn(i) = aire( i, 1   ) * masse(  i     ,  1   , l )
    290289          xpps(i) = aire( i,jjp1 ) * masse(  i     , jjp1 , l )
    291          ENDDO
     290        ENDDO
    292291          xpn      = SUM(xppn)/apoln
    293292          xps      = SUM(xpps)/apols
    294          DO i   = 1, iip1
     293        DO i   = 1, iip1
    295294          masse(   i   ,   1     ,  l )   = xpn
    296295          masse(   i   ,   jjp1  ,  l )   = xps
    297          ENDDO
    298        ENDDO
    299       ENDIF
     296        ENDDO
     297      ENDDO
    300298      q3d(iip1,:,:,:) = q3d(1,:,:,:)
    301299      phis(iip1,:) = phis(1,:)

  • LMDZ.3.3/trunk/libf/dyn3d/exner_hyb.F

    r2 r50  
    5353c
    5454      DO   ij  = 1, ngrid
    55        pks(ij) = cpp * ( ps(ij)/preff ) ** kappa
     55        pks(ij) = cpp * ( ps(ij)/preff ) ** kappa
    5656      ENDDO
    5757
    58       IF( alphax.NE.0. )   THEN
    59         DO  ij   = 1, iim
    60          ppn(ij) = aire(   ij   ) * pks(  ij     )
    61          pps(ij) = aire(ij+ip1jm) * pks(ij+ip1jm )
    62         ENDDO
    63          xpn      = SSUM(iim,ppn,1) /apoln
    64          xps      = SSUM(iim,pps,1) /apols
     58      DO  ij   = 1, iim
     59        ppn(ij) = aire(   ij   ) * pks(  ij     )
     60        pps(ij) = aire(ij+ip1jm) * pks(ij+ip1jm )
     61      ENDDO
     62      xpn      = SSUM(iim,ppn,1) /apoln
     63      xps      = SSUM(iim,pps,1) /apols
    6564
    66         DO ij   = 1, iip1
    67          pks(   ij     )  =  xpn
    68          pks( ij+ip1jm )  =  xps
    69         ENDDO
    70       ENDIF
     65      DO ij   = 1, iip1
     66        pks(   ij     )  =  xpn
     67        pks( ij+ip1jm )  =  xps
     68      ENDDO
    7169c
    7270c
Note: See TracChangeset for help on using the changeset viewer.