Changeset 1046 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
Sep 21, 2013, 8:33:22 AM (12 years ago)
Author:
emillour
Message:

Mars GCM:

  • Bug fix in yamada4.F: values for near-surface diffusion coefficients were used but not set.
  • NB: FH recommends updating to latest yamada4 from Earth GCM and using their iflag_pbl=11 scheme.

EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/libf/phymars/yamada4.F

    r1036 r1046  
    464464      enddo
    465465      enddo
    466       km(:,nlay+1)=km(:,nlay)
    467       kn(:,nlay+1)=kn(:,nlay)
    468       kq(:,nlay+1)=kq(:,nlay)
     466      ! boundary conditions for km
     467      km(:,nlay+1)=0
     468      km(:,1)=km(:,2) ! km(:,1)=0
     469      ! boundary conditions for kn
     470      kn(:,nlay+1)=0
     471      kn(:,1)=kn(:,2) ! kn(:,1)=0
     472      ! boundary conditions for kq
     473      kq(:,nlay+1)=0  ! zero at top of atmosphere
     474      kq(:,1)=kq(:,2) ! no gradient at surface
    469475
    470476      q2(:,1)=q2(:,2)
Note: See TracChangeset for help on using the changeset viewer.