Changeset 1403 for LMDZ4/trunk/libf/bibio/writehist.F
- Timestamp:
- Jul 1, 2010, 11:02:53 AM (14 years ago)
- Location:
- LMDZ4/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk
- Property svn:mergeinfo changed
-
LMDZ4/trunk/libf/bibio/writehist.F
r1279 r1403 2 2 ! $Id$ 3 3 ! 4 subroutine writehist( histid, histvid, time, vcov, 5 , ucov,teta,phi,q,masse,ps,phis) 4 subroutine writehist(time,vcov,ucov,teta,phi,q,masse,ps,phis) 6 5 7 6 #ifdef CPP_IOIPSL … … 9 8 #endif 10 9 USE infotrac, ONLY : nqtot, ttext 10 use com_io_dyn_mod, only : histid,histvid,histuid 11 11 implicit none 12 12 … … 17 17 C 18 18 C Entree: 19 C histid: ID du fichier histoire20 C histvid:ID du fichier histoire pour les vents V (appele a disparaitre)21 19 C time: temps de l'ecriture 22 20 C vcov: vents v covariants … … 29 27 C phis : geopotentiel au sol 30 28 C 31 C32 C Sortie:33 C fileid: ID du fichier netcdf cree34 29 C 35 30 C L. Fairhead, LMD, 03/99 … … 54 49 C 55 50 56 INTEGER histid, histvid57 51 REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) 58 52 REAL teta(ip1jmp1,llm),phi(ip1jmp1,llm) … … 71 65 logical ok_sync 72 66 integer itau_w 67 REAL vnat(ip1jm,llm),unat(ip1jmp1,llm) 68 73 69 C 74 70 C Initialisations … … 79 75 ok_sync =.TRUE. 80 76 itau_w = itau_dyn + time 77 ! Passage aux composantes naturelles du vent 78 call covnat(llm, ucov, vcov, unat, vnat) 81 79 C 82 80 C Appels a histwrite pour l'ecriture des variables a sauvegarder … … 84 82 C Vents U 85 83 C 86 call histwrite(hist id, 'ucov', itau_w, ucov,84 call histwrite(histuid, 'u', itau_w, unat, 87 85 . iip1*jjp1*llm, ndexu) 88 89 86 C 90 87 C Vents V 91 88 C 92 call histwrite(histvid, 'v cov', itau_w, vcov,89 call histwrite(histvid, 'v', itau_w, vnat, 93 90 . iip1*jjm*llm, ndexv) 94 91 … … 106 103 C Traceurs 107 104 C 108 DO iq=1,nqtot109 call histwrite(histid, ttext(iq), itau_w, q(:,:,iq),110 . iip1*jjp1*llm, ndexu)111 enddo112 C105 ! DO iq=1,nqtot 106 ! call histwrite(histid, ttext(iq), itau_w, q(:,:,iq), 107 ! . iip1*jjp1*llm, ndexu) 108 ! enddo 109 !C 113 110 C Masse 114 111 C 115 call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)112 call histwrite(histid,'masse',itau_w, masse,iip1*jjp1*llm,ndexu) 116 113 C 117 114 C Pression au sol … … 121 118 C Geopotentiel au sol 122 119 C 123 call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)120 ! call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d) 124 121 C 125 122 C Fin … … 128 125 call histsync(histid) 129 126 call histsync(histvid) 127 call histsync(histuid) 130 128 endif 131 129 #else
Note: See TracChangeset
for help on using the changeset viewer.