Changeset 938 for LMDZ4/trunk/libf/dyn3d
- Timestamp:
- Mar 31, 2008, 6:14:45 PM (17 years ago)
- Location:
- LMDZ4/trunk/libf/dyn3d
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/dyn3d/guide.F
r762 r938 122 122 c calcul de l'humidite saturante 123 123 C----------------------------------------------------------------------- 124 124 c print*,'OK0' 125 125 CALL pression( ip1jmp1, ap, bp, ps, p ) 126 126 call massdair(p,masse) 127 127 c print*,'OK1' 128 128 CALL exner_hyb(ip1jmp1,ps,p,alpha,beta,pks,pk,pkf) 129 129 c print*,'OK2' 130 130 tnat(:,:)=pk(:,:)*teta(:,:)/cpp 131 131 c print*,'OK3' 132 132 unskap = 1./ kappa 133 133 pres(:,:)=preff*(pk(:,:)/cpp)**unskap 134 134 c print*,'OK4' 135 135 call q_sat(iip1*jjp1*llm,tnat,pres,qsat) 136 136 … … 144 144 c----------------------------------------------------------------------- 145 145 146 146 c print*,'ONLINE=',online 147 147 if(online.eq.-1) then 148 148 return -
LMDZ4/trunk/libf/dyn3d/wrgrads.F
r541 r938 29 29 writectl=.false. 30 30 31 31 c print*,if,iid(if),jid(if),ifd(if),jfd(if) 32 32 iii=iid(if) 33 33 iji=jid(if) … … 38 38 lm=lmd(if) 39 39 40 41 40 c print*,'im,jm,lm,name,firsttime(if)' 41 c print*,im,jm,lm,name,firsttime(if) 42 42 43 43 if(firsttime(if)) then … … 57 57 tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar)) 58 58 nld(ivar(if),if)=nl 59 60 59 c print*,'initialisation ecriture de ',var(ivar(if),if) 60 c print*,'if ivar(if) nld ',if,ivar(if),nld(ivar(if),if) 61 61 endif 62 62 writectl=.true. … … 81 81 endif 82 82 83 84 83 c print*,'ivar(if),nvar(if),var(ivar(if),if),writectl' 84 c print*,ivar(if),nvar(if),var(ivar(if),if),writectl 85 85 do l=1,nl 86 86 irec(if)=irec(if)+1
Note: See TracChangeset
for help on using the changeset viewer.