Changeset 143 for LMDZ.3.3


Ignore:
Timestamp:
Oct 19, 2000, 5:23:00 PM (24 years ago)
Author:
lmdzadmin
Message:

Corrections bugs couplage
LF

Location:
LMDZ.3.3/branches/rel-LF/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

    r112 r143  
    313313     .            yts, yrugos, yu, yv, yt, yq,
    314314     .            ycoefm, ycoefh)
    315       CALL coefkz2(nsrf, knon, paprs, pplay,t,
     315      CALL coefkz2(nsrf, knon, ypaprs, ypplay,yt,
    316316     .                  ycoefm0, ycoefh0)
    317317      DO k = 1, klev
     
    705705     s evap, fluxsens, fluxlat, dflux_l, dflux_s,             
    706706     s tsol_rad, tsurf_new, alb_new, emis_new, z0_new,
    707      s pctsrf_NEW, agesno)
     707     s pctsrf_new, agesno)
    708708
    709709      flux_t(:,1) = fluxsens
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

    r140 r143  
    792792      call abort_gcm(modname,abort_message,1)
    793793    endif
    794     allocate(pctsrf_sav(klon,2), stat = error)
     794    allocate(pctsrf_sav(klon,nbsrf), stat = error)
    795795    if ( error  /=0) then
    796796      abort_message='Pb allocation variable pctsrf_sav'
    797797      call abort_gcm(modname,abort_message,1)
    798798    endif
     799    pctsrf_sav = 0.
    799800
    800801    do ig = 1, klon
     
    909910      IF (pctsrf(ig,is_oce) > epsfra .OR.            &
    910911     &             pctsrf(ig,is_sic) > epsfra) THEN
    911             pctsrf_sav(ig,is_sic) = tamp_sic(ig) * pctsrf(ig,is_sic)
    912             pctsrf_sav(ig,is_oce) = pctsrf(ig,is_oce)    &
    913      &                        - (pctsrf_sav(ig,is_sic)-pctsrf(ig,is_sic))
     912          pctsrf_sav(ig,is_sic) = (pctsrf(ig,is_oce) + pctsrf(ig,is_sic)) &
     913     &                               * tamp_sic(ig)
     914          pctsrf_sav(ig,is_oce) = (pctsrf(ig,is_oce) + pctsrf(ig,is_sic)) &
     915     &                        - pctsrf_sav(ig,is_sic)
     916        endif
     917      enddo
     918      if (check) then
     919        write(46,*)'pct_srf_sav_ice = '
     920        write(46,'(72f8.3)')pctsrf_sav(:,is_sic)
    914921      endif
    915     enddo
    916     if (minval(pctsrf_new(:,is_oce)) < 0.) then
    917       write(*,*)'Pb fraction ocean inferieure a 0'
    918       write(*,*)'au point ',minloc(pctsrf_new(:,is_oce))
    919       write(*,*)'valeur = ',minval(pctsrf_new(:,is_oce))
    920       abort_message = 'voir ci-dessus'
    921       call abort_gcm(modname,abort_message,1)
    922     endif
    923     if (minval(pctsrf_new(:,is_sic)) < 0.) then
    924       write(*,*)'Pb fraction glace inferieure a 0'
    925       write(*,*)'au point ',minloc(pctsrf_new(:,is_sic))
    926       write(*,*)'valeur = ',minval(pctsrf_new(:,is_sic))
    927       abort_message = 'voir ci-dessus'
    928       call abort_gcm(modname,abort_message,1)
    929     endif
    930        if (check) then
    931          write(47,*)'Sortie fromcpl apres bidouille'
    932          write(47,*)' read_sst = '
    933          write(47,'(72f8.3)')read_sst
    934          call flush(47)
    935 !        allocate(isst(iim, jjm+1), stat = error)
    936 !        isst = 0
    937 !        where (read_sst >0.) isst = 1
    938 !        write(46,'(72i1)')isst
    939        endif
     922      if (minval(pctsrf_new(:,is_oce)) < 0.) then
     923        write(*,*)'Pb fraction ocean inferieure a 0'
     924        write(*,*)'au point ',minloc(pctsrf_new(:,is_oce))
     925        write(*,*)'valeur = ',minval(pctsrf_new(:,is_oce))
     926        abort_message = 'voir ci-dessus'
     927        call abort_gcm(modname,abort_message,1)
     928      endif
     929      if (minval(pctsrf_new(:,is_sic)) < 0.) then
     930        write(*,*)'Pb fraction glace inferieure a 0'
     931        write(*,*)'au point ',minloc(pctsrf_new(:,is_sic))
     932        write(*,*)'valeur = ',minval(pctsrf_new(:,is_sic))
     933        abort_message = 'voir ci-dessus'
     934        call abort_gcm(modname,abort_message,1)
     935      endif
    940936    endif
    941937  endif                         ! fin mod(itime, nexca) == 1
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r132 r143  
    121121c======================================================================
    122122      LOGICAL ok_journe ! sortir le fichier journalier
    123       PARAMETER (ok_journe=.FALSE.)
     123      PARAMETER (ok_journe=.true.)
    124124c
    125125      LOGICAL ok_mensuel ! sortir le fichier mensuel
     
    127127c
    128128      LOGICAL ok_instan ! sortir le fichier instantane
    129       PARAMETER (ok_instan=.false.)
     129      PARAMETER (ok_instan=.true.)
    130130c
    131131      LOGICAL ok_region ! sortir le fichier regional
Note: See TracChangeset for help on using the changeset viewer.