Changeset 1608 for LMDZ5


Ignore:
Timestamp:
Jan 10, 2012, 2:33:09 PM (13 years ago)
Author:
lguez
Message:

In "phy1d/*", replaced obsolete calls to "lnblnk" by calls to "trim".

In "lmdz1d.F", use collector module "ioipsl" rather than specific
module "calendar", which is an internal module of IOIPSL. Removed line
containing only "#" (causes compilation error). Bug fix in call to
"init_phys_lmdz".

Location:
LMDZ5/trunk/libf/phy1d
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phy1d/1DUTILS.h_no_writelim

    r1607 r1608  
    735735!   local
    736736 
    737       integer im,jm,lm,i,j,l,lnblnk,iv,iii,iji,iif,ijf
     737      integer im,jm,lm,i,j,l,iv,iii,iji,iif,ijf
    738738 
    739739      logical writectl
    … …  
    769769            nvar(if)=ivar(if)
    770770            var(ivar(if),if)=name
    771             tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar))
     771            tvar(ivar(if),if)=trim(titlevar)
    772772            nld(ivar(if),if)=nl
    773773            print*,'initialisation ecriture de ',var(ivar(if),if)
    … …  
    810810      file=fichier(if)
    811811!   WARNING! on reecrase le fichier .ctl a chaque ecriture
    812       open(unit(if),file=file(1:lnblnk(file))//'.ctl',
     812      open(unit(if),file=trim(file)//'.ctl',
    813813     &         form='formatted',status='unknown')
    814814      write(unit(if),'(a5,1x,a40)')
    815      &       'DSET ','^'//file(1:lnblnk(file))//'.dat'
     815     &       'DSET ','^'//trim(file)//'.dat'
    816816 
    817817      write(unit(if),'(a12)') 'UNDEF 1.0E30'
    … …  
    848848      implicit none
    849849 
    850       integer if,im,jm,lm,i,j,l,lnblnk
     850      integer if,im,jm,lm,i,j,l
    851851      real x(im),y(jm),z(lm),fx,fy,fz,dt
    852852      real xmin,xmax,ymin,ymax
    … …  
    868868      ivar(if)=0
    869869 
    870       fichier(if)=file(1:lnblnk(file))
     870      fichier(if)=trim(file)
    871871 
    872872      firsttime(if)=.true.
    … …  
    898898 
    899899      print*,4*(ifd(if)-iid(if))*(jfd(if)-jid(if))
    900       print*,file(1:lnblnk(file))//'.dat'
    901  
    902       OPEN (unit(if)+1,FILE=file(1:lnblnk(file))//'.dat',
     900      print*,trim(file)//'.dat'
     901 
     902      OPEN (unit(if)+1,FILE=trim(file)//'.dat',
    903903     s   FORM='UNFORMATTED',
    904904     s   ACCESS='DIRECT'
    … …  
    13601360      integer iy,im,id,ih
    13611361     
    1362       INTEGER, EXTERNAL :: lnblnk
    1363 
    13641362       real plev_min
    13651363
    13661364       plev_min = 55.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
    13671365
    1368       open(21,file=fich_toga(:lnblnk(fich_toga)),form='formatted')
     1366      open(21,file=trim(fich_toga),form='formatted')
    13691367      read(21,'(a)')
    13701368      do ip = 1, nt_toga
    … …  
    19181916
    19191917     
    1920       INTEGER, EXTERNAL :: lnblnk
    1921 
    19221918      print*,fich_rico
    1923       open(21,file=fich_rico(:lnblnk(fich_rico)),form='formatted')
     1919      open(21,file=trim(fich_rico),form='formatted')
    19241920        do k=1,llm
    19251921      zlay(k)=0.
    … …  
    20722068      integer iy,im,id,ih,in
    20732069
    2074       INTEGER, EXTERNAL :: lnblnk
    2075 
    2076       open(21,file=fich_armcu(:lnblnk(fich_armcu)),form='formatted')
     2070      open(21,file=trim(fich_armcu),form='formatted')
    20772071      read(21,'(a)')
    20782072      do ip = 1, nt_armcu

  • LMDZ5/trunk/libf/phy1d/1DUTILS.h_with_writelim

    r1607 r1608  
    735735!   local
    736736 
    737       integer im,jm,lm,i,j,l,lnblnk,iv,iii,iji,iif,ijf
     737      integer im,jm,lm,i,j,l,iv,iii,iji,iif,ijf
    738738 
    739739      logical writectl
    … …  
    769769            nvar(if)=ivar(if)
    770770            var(ivar(if),if)=name
    771             tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar))
     771            tvar(ivar(if),if)=trim(titlevar)
    772772            nld(ivar(if),if)=nl
    773773            print*,'initialisation ecriture de ',var(ivar(if),if)
    … …  
    810810      file=fichier(if)
    811811!   WARNING! on reecrase le fichier .ctl a chaque ecriture
    812       open(unit(if),file=file(1:lnblnk(file))//'.ctl',
     812      open(unit(if),file=trim(file)//'.ctl',
    813813     &         form='formatted',status='unknown')
    814814      write(unit(if),'(a5,1x,a40)')
    815      &       'DSET ','^'//file(1:lnblnk(file))//'.dat'
     815     &       'DSET ','^'//trim(file)//'.dat'
    816816 
    817817      write(unit(if),'(a12)') 'UNDEF 1.0E30'
    … …  
    848848      implicit none
    849849 
    850       integer if,im,jm,lm,i,j,l,lnblnk
     850      integer if,im,jm,lm,i,j,l
    851851      real x(im),y(jm),z(lm),fx,fy,fz,dt
    852852      real xmin,xmax,ymin,ymax
    … …  
    868868      ivar(if)=0
    869869 
    870       fichier(if)=file(1:lnblnk(file))
     870      fichier(if)=trim(file)
    871871 
    872872      firsttime(if)=.true.
    … …  
    898898 
    899899      print*,4*(ifd(if)-iid(if))*(jfd(if)-jid(if))
    900       print*,file(1:lnblnk(file))//'.dat'
    901  
    902       OPEN (unit(if)+1,FILE=file(1:lnblnk(file))//'.dat',
     900      print*,trim(file)//'.dat'
     901 
     902      OPEN (unit(if)+1,FILE=trim(file)//'.dat',
    903903     s   FORM='UNFORMATTED',
    904904     s   ACCESS='DIRECT'
    … …  
    14831483      integer iy,im,id,ih
    14841484     
    1485       INTEGER, EXTERNAL :: lnblnk
    1486 
    14871485       real plev_min
    14881486
    14891487       plev_min = 55.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
    14901488
    1491       open(21,file=fich_toga(:lnblnk(fich_toga)),form='formatted')
     1489      open(21,file=trim(fich_toga),form='formatted')
    14921490      read(21,'(a)')
    14931491      do ip = 1, nt_toga
    … …  
    20412039
    20422040     
    2043       INTEGER, EXTERNAL :: lnblnk
    2044 
    20452041      print*,fich_rico
    2046       open(21,file=fich_rico(:lnblnk(fich_rico)),form='formatted')
     2042      open(21,file=trim(fich_rico),form='formatted')
    20472043        do k=1,llm
    20482044      zlay(k)=0.
    … …  
    21952191      integer iy,im,id,ih,in
    21962192
    2197       INTEGER, EXTERNAL :: lnblnk
    2198 
    2199       open(21,file=fich_armcu(:lnblnk(fich_armcu)),form='formatted')
     2193      open(21,file=trim(fich_armcu),form='formatted')
    22002194      read(21,'(a)')
    22012195      do ip = 1, nt_armcu

  • LMDZ5/trunk/libf/phy1d/1DUTILS.h_with_writelim_old

    r1607 r1608  
    735735!   local
    736736 
    737       integer im,jm,lm,i,j,l,lnblnk,iv,iii,iji,iif,ijf
     737      integer im,jm,lm,i,j,l,iv,iii,iji,iif,ijf
    738738 
    739739      logical writectl
    … …  
    769769            nvar(if)=ivar(if)
    770770            var(ivar(if),if)=name
    771             tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar))
     771            tvar(ivar(if),if)=trim(titlevar)
    772772            nld(ivar(if),if)=nl
    773773            print*,'initialisation ecriture de ',var(ivar(if),if)
    … …  
    810810      file=fichier(if)
    811811!   WARNING! on reecrase le fichier .ctl a chaque ecriture
    812       open(unit(if),file=file(1:lnblnk(file))//'.ctl',
     812      open(unit(if),file=trim(file)//'.ctl',
    813813     &         form='formatted',status='unknown')
    814814      write(unit(if),'(a5,1x,a40)')
    815      &       'DSET ','^'//file(1:lnblnk(file))//'.dat'
     815     &       'DSET ','^'//trim(file)//'.dat'
    816816 
    817817      write(unit(if),'(a12)') 'UNDEF 1.0E30'
    … …  
    848848      implicit none
    849849 
    850       integer if,im,jm,lm,i,j,l,lnblnk
     850      integer if,im,jm,lm,i,j,l
    851851      real x(im),y(jm),z(lm),fx,fy,fz,dt
    852852      real xmin,xmax,ymin,ymax
    … …  
    868868      ivar(if)=0
    869869 
    870       fichier(if)=file(1:lnblnk(file))
     870      fichier(if)=trim(file)
    871871 
    872872      firsttime(if)=.true.
    … …  
    898898 
    899899      print*,4*(ifd(if)-iid(if))*(jfd(if)-jid(if))
    900       print*,file(1:lnblnk(file))//'.dat'
    901  
    902       OPEN (unit(if)+1,FILE=file(1:lnblnk(file))//'.dat',
     900      print*,trim(file)//'.dat'
     901 
     902      OPEN (unit(if)+1,FILE=trim(file)//'.dat',
    903903     s   FORM='UNFORMATTED',
    904904     s   ACCESS='DIRECT'
    … …  
    14831483      integer iy,im,id,ih
    14841484     
    1485       INTEGER, EXTERNAL :: lnblnk
    14861485
    14871486       real plev_min
    … …  
    14891488       plev_min = 55.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
    14901489
    1491       open(21,file=fich_toga(:lnblnk(fich_toga)),form='formatted')
     1490      open(21,file=trim(fich_toga),form='formatted')
    14921491      read(21,'(a)')
    14931492      do ip = 1, nt_toga
    … …  
    20412040
    20422041     
    2043       INTEGER, EXTERNAL :: lnblnk
    2044 
    20452042      print*,fich_rico
    2046       open(21,file=fich_rico(:lnblnk(fich_rico)),form='formatted')
     2043      open(21,file=trim(fich_rico),form='formatted')
    20472044        do k=1,llm
    20482045      zlay(k)=0.
    … …  
    21952192      integer iy,im,id,ih,in
    21962193
    2197       INTEGER, EXTERNAL :: lnblnk
    2198 
    2199       open(21,file=fich_armcu(:lnblnk(fich_armcu)),form='formatted')
     2194      open(21,file=trim(fich_armcu),form='formatted')
    22002195      read(21,'(a)')
    22012196      do ip = 1, nt_armcu

  • LMDZ5/trunk/libf/phy1d/lmdz1d.F

    r1607 r1608  
    11      PROGRAM lmdz1d
    22
    3       USE ioipsl, only: ju2ymds, ymds2ju
    4       USE calendar
     3      USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar
    54      use phys_state_var_mod
    65      use comgeomphy
    … …  
    298297
    299298c      calend = 'earth_365d'
    300 #
    301299      if (calend == 'earth_360d') then
    302300        call ioconf_calendar('360d')
    … …  
    358356! Initialization of dimensions, geometry and initial state
    359357!---------------------------------------------------------------------
    360       call init_phys_lmdz(1,1,llm,1,1)
     358      call init_phys_lmdz(1,1,llm,1,(/1/))
    361359      call suphel
    362360      call initcomgeomphy
Note: See TracChangeset for help on using the changeset viewer.