Context Navigation

ini_archive.F @ 1720

Last change on this file since 1720 was 1422, checked in by milmd, 10 years ago
In GENERIC, MARS and COMMON models replace some include files by modules (usefull for decoupling physics with dynamics).
File size: 17.0 KB

Rev	Line
[38]	1	c=======================================================================
	2	subroutine ini_archive(nid,idayref,phis,ith,tab_cntrl_fi)
	3	c=======================================================================
	4	c
	5	c
	6	c Date: 01/1997
	7	c ----
	8	c
	9	c Objet: ecriture de l'entete du fichier "start_archive"
	10	c -----
	11	c
	12	c Proche de iniwrite.F
	13	c
	14	c On ajoute dans le tableau "tab_cntrl" (dynamique), a partir de 51,
	15	c les valeurs de tab_cntrl_fi (les 38 parametres de controle physiques
	16	c du RUN + ptotal et cotoicetotal)
	17	c
	18	c tab_cntrl(50+l)=tab_cntrl_fi(l)
	19	c
	20	c Arguments:
	21	c ---------
	22	c
	23	c Inputs:
	24	c ------
	25	c
	26	c nid unite logique du fichier "start_archive"
	27	c idayref Valeur du jour initial a mettre dans
	28	c l'entete du fichier "start_archive"
	29	c phis geopotentiel au sol
	30	c ith soil thermal inertia
	31	c tab_cntrl_fi tableau des param physiques
	32	c
	33
	34	c=======================================================================
	35
[1047]	36	use comsoil_h, only: nsoilmx, mlayer
[1422]	37	USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff,presnivs,pseudoalt
	38	USE comconst_mod, ONLY: daysec,dtvr,rad,omeg,g,cpp,kappa,pi
	39	USE logic_mod, ONLY: fxyhypb,ysinus
	40	USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy
	41	USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
[38]	42	implicit none
	43
	44	#include "dimensions.h"
	45	#include "paramet.h"
	46	#include "comgeom.h"
	47	#include "netcdf.inc"
	48
	49	c-----------------------------------------------------------------------
	50	c Declarations
	51	c-----------------------------------------------------------------------
	52
	53	c Local:
	54	c ------
	55	INTEGER length,l
	56	parameter (length = 100)
	57	REAL tab_cntrl(length) ! tableau des parametres du run
	58	INTEGER loop
	59	INTEGER ierr, setvdim, putvdim, putdat, setname,cluvdb
	60	INTEGER setdim
	61	INTEGER ind1,indlast
	62
	63	c Arguments:
	64	c ----------
	65	INTEGER*4 idayref
	66	REAL phis(ip1jmp1)
	67	real ith(ip1jmp1,nsoilmx)
	68	REAL tab_cntrl_fi(length)
	69
	70	!Mars --------Ajouts-----------
	71	c Variables locales pour NetCDF:
	72	c
	73	INTEGER dims2(2), dims3(3) !, dims4(4)
	74	INTEGER idim_index
	75	INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
	76	INTEGER idim_llmp1,idim_llm
	77	INTEGER idim_tim
	78	INTEGER idim_nsoilmx ! "subsurface_layers" dimension ID #
	79	INTEGER nid,nvarid
	80	real sig_s(llm),s(llm)
	81
	82	pi = 2. * ASIN(1.)
	83
	84
	85	c-----------------------------------------------------------------------
	86	c Remplissage du tableau des parametres de controle du RUN (dynamique)
	87	c-----------------------------------------------------------------------
	88
	89	DO l=1,length
	90	tab_cntrl(l)=0.
	91	ENDDO
	92
	93	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
	94	tab_cntrl(1) = REAL(iim) ! nombre de points en longitude
	95	tab_cntrl(2) = REAL(jjm) ! nombre de points en latitude
	96	tab_cntrl(3) = REAL(llm) ! nombre de couches
	97	tab_cntrl(4) = REAL(idayref) ! jour 0
	98	tab_cntrl(5) = rad ! rayon de mars(m) ~3397200
	99	tab_cntrl(6) = omeg ! vitesse de rotation (rad.s-1)
	100	tab_cntrl(7) = g ! gravite (m.s-2) ~3.72
	101	tab_cntrl(8) = cpp
	102	tab_cntrl(8) = 43.49 !mars temporaire Masse molaire de l''atm (g.mol-1) ~43.49
	103	tab_cntrl(9) = kappa ! = r/cp ~0.256793 (=rcp dans physique)
	104	tab_cntrl(10) = daysec ! duree du sol (s) ~88775
	105	tab_cntrl(11) = dtvr ! pas de temps de la dynamique (s)
	106	tab_cntrl(12) = etot0 ! energie totale !
	107	tab_cntrl(13) = ptot0 ! pression totalei ! variables
	108	tab_cntrl(14) = ztot0 ! enstrophie totale ! de controle
	109	tab_cntrl(15) = stot0 ! enthalpie totale ! globales
	110	tab_cntrl(16) = ang0 ! moment cinetique !
	111	tab_cntrl(17) = pa
	112	tab_cntrl(18) = preff
	113
	114	c ..... parametres pour le zoom ......
	115
	116	tab_cntrl(19) = clon ! longitude en degres du centre du zoom
	117	tab_cntrl(20) = clat ! latitude en degres du centre du zoom
	118	tab_cntrl(21) = grossismx ! facteur de grossissement du zoom,selon longitude
	119	tab_cntrl(22) = grossismy ! facteur de grossissement du zoom ,selon latitude
	120
	121	IF ( fxyhypb ) THEN
	122	tab_cntrl(23) = 1.
	123	tab_cntrl(24) = dzoomx ! extension en longitude de la zone du zoom
	124	tab_cntrl(25) = dzoomy ! extension en latitude de la zone du zoom
	125	ELSE
	126	tab_cntrl(23) = 0.
	127	tab_cntrl(24) = dzoomx ! extension en longitude de la zone du zoom
	128	tab_cntrl(25) = dzoomy ! extension en latitude de la zone du zoom
	129	tab_cntrl(26) = 0.
	130	IF ( ysinus) tab_cntrl(26) = 1.
	131	ENDIF
	132
	133	c-----------------------------------------------------------------------
	134	c Copie du tableau des parametres de controle du RUN (physique)
	135	c dans le tableau dynamique
	136	c-----------------------------------------------------------------------
	137
	138	DO l=1,50
	139	tab_cntrl(50+l)=tab_cntrl_fi(l)
	140	ENDDO
	141
	142	c=======================================================================
	143	c Ecriture NetCDF de l''entete du fichier "start_archive"
	144	c=======================================================================
	145
	146	c
	147	c Preciser quelques attributs globaux:
	148	c
	149	ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 21,
	150	& "Fichier start_archive")
	151	c
	152	c Definir les dimensions du fichiers:
	153	c
	154	c CHAMPS AJOUTES POUR LA VISUALISATION T,ps, etc... avec Grads ou ferret:
	155	ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_rlatu)
	156	ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_rlonv)
	157	ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
	158	ierr = NF_DEF_DIM (nid,"subsurface_layers",nsoilmx,idim_nsoilmx)
	159
	160	ierr = NF_DEF_DIM (nid,"index", length, idim_index)
	161	ierr = NF_DEF_DIM (nid,"rlonu", iip1, idim_rlonu)
	162	ierr = NF_DEF_DIM (nid,"rlatv", jjm, idim_rlatv)
	163	ierr = NF_DEF_DIM (nid,"interlayer", llmp1, idim_llmp1)
	164	ierr = NF_DEF_DIM (nid,"Time", NF_UNLIMITED, idim_tim)
	165
	166	c
	167	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
	168
	169	c-----------------------------------------------------------------------
	170	c Ecriture du tableau des parametres du run
	171	c-----------------------------------------------------------------------
	172
	173	call def_var(nid,"Time","Time","days since 00:00:00",1,
	174	. idim_tim,nvarid,ierr)
	175
	176	ierr = NF_REDEF (nid)
	177	#ifdef NC_DOUBLE
	178	ierr = NF_DEF_VAR (nid,"controle",NF_DOUBLE,1,idim_index,nvarid)
	179	#else
	180	ierr = NF_DEF_VAR (nid,"controle",NF_FLOAT,1,idim_index,nvarid)
	181	#endif
	182	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	183	. "Parametres de controle")
	184	ierr = NF_ENDDEF(nid)
	185	#ifdef NC_DOUBLE
	186	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
	187	#else
	188	ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
	189	#endif
	190
	191	c-----------------------------------------------------------------------
	192	c Ecriture des longitudes et latitudes
	193	c-----------------------------------------------------------------------
	194
	195	ierr = NF_REDEF (nid)
	196	#ifdef NC_DOUBLE
	197	ierr = NF_DEF_VAR (nid,"rlonu",NF_DOUBLE,1,idim_rlonu,nvarid)
	198	#else
	199	ierr = NF_DEF_VAR (nid,"rlonu",NF_FLOAT,1,idim_rlonu,nvarid)
	200	#endif
	201	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
	202	. "Longitudes des points U")
	203	ierr = NF_ENDDEF(nid)
	204	#ifdef NC_DOUBLE
	205	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonu)
	206	#else
	207	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu)
	208	#endif
	209	c
	210	ierr = NF_REDEF (nid)
	211	#ifdef NC_DOUBLE
	212	ierr = NF_DEF_VAR (nid,"rlatu",NF_DOUBLE,1,idim_rlatu,nvarid)
	213	#else
	214	ierr = NF_DEF_VAR (nid,"rlatu",NF_FLOAT,1,idim_rlatu,nvarid)
	215	#endif
	216	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	217	. "Latitudes des points U")
	218	ierr = NF_ENDDEF(nid)
	219	#ifdef NC_DOUBLE
	220	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu)
	221	#else
	222	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu)
	223	#endif
	224	c
	225	ierr = NF_REDEF (nid)
	226	#ifdef NC_DOUBLE
	227	ierr = NF_DEF_VAR (nid,"rlonv",NF_DOUBLE,1,idim_rlonv,nvarid)
	228	#else
	229	ierr = NF_DEF_VAR (nid,"rlonv",NF_FLOAT,1,idim_rlonv,nvarid)
	230	#endif
	231	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
	232	. "Longitudes des points V")
	233	ierr = NF_ENDDEF(nid)
	234	#ifdef NC_DOUBLE
	235	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv)
	236	#else
	237	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv)
	238	#endif
	239	c
	240	ierr = NF_REDEF (nid)
	241	#ifdef NC_DOUBLE
	242	ierr = NF_DEF_VAR (nid,"rlatv",NF_DOUBLE,1,idim_rlatv,nvarid)
	243	#else
	244	ierr = NF_DEF_VAR (nid,"rlatv",NF_FLOAT,1,idim_rlatv,nvarid)
	245	#endif
	246	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	247	. "Latitudes des points V")
	248	ierr = NF_ENDDEF(nid)
	249	#ifdef NC_DOUBLE
	250	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatv)
	251	#else
	252	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv)
	253	#endif
	254
	255	c-----------------------------------------------------------------------
	256	c Ecriture des niveaux verticaux
	257	c-----------------------------------------------------------------------
	258
	259	c
	260	ierr = NF_REDEF (nid)
	261	#ifdef NC_DOUBLE
	262	ierr = NF_DEF_VAR (nid,"ap",NF_DOUBLE,1,idim_llmp1,nvarid)
	263	#else
	264	ierr = NF_DEF_VAR (nid,"ap",NF_FLOAT,1,idim_llmp1,nvarid)
	265	#endif
	266	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 32,
	267	. "Coef A: niveaux pression hybride")
	268	ierr = NF_ENDDEF(nid)
	269	#ifdef NC_DOUBLE
	270	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ap)
	271	#else
	272	ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)
	273	#endif
	274	c
	275	ierr = NF_REDEF (nid)
	276	#ifdef NC_DOUBLE
	277	ierr = NF_DEF_VAR (nid,"bp",NF_DOUBLE,1,idim_llmp1,nvarid)
	278	#else
	279	ierr = NF_DEF_VAR (nid,"bp",NF_FLOAT,1,idim_llmp1,nvarid)
	280	#endif
	281	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 35,
	282	. "Coefficient B niveaux sigma hybride")
	283	ierr = NF_ENDDEF(nid)
	284	#ifdef NC_DOUBLE
	285	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bp)
	286	#else
	287	ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)
	288	#endif
	289	c
	290	c ----------------------
	291	ierr = NF_REDEF (nid)
	292	#ifdef NC_DOUBLE
	293	ierr = NF_DEF_VAR (nid,"aps",NF_DOUBLE,1,idim_llm,nvarid)
	294	#else
	295	ierr = NF_DEF_VAR (nid,"aps",NF_FLOAT,1,idim_llm,nvarid)
	296	#endif
	297	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 36,
	298	. "Coef AS: hybrid pressure in midlayers")
	299	ierr = NF_ENDDEF(nid)
	300	#ifdef NC_DOUBLE
	301	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
	302	#else
	303	ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
	304	#endif
	305	c
	306	c ----------------------
	307	ierr = NF_REDEF (nid)
	308	#ifdef NC_DOUBLE
	309	ierr = NF_DEF_VAR (nid,"bps",NF_DOUBLE,1,idim_llm,nvarid)
	310	#else
	311	ierr = NF_DEF_VAR (nid,"bps",NF_FLOAT,1,idim_llm,nvarid)
	312	#endif
	313	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 30,
	314	. "Coef BS: hybrid sigma midlayers")
	315	ierr = NF_ENDDEF(nid)
	316	#ifdef NC_DOUBLE
	317	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
	318	#else
	319	ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
	320	#endif
	321	c
	322	c ----------------------
	323
	324	ierr = NF_REDEF (nid)
	325	#ifdef NC_DOUBLE
	326	ierr = NF_DEF_VAR (nid,"presnivs",NF_DOUBLE,1,idim_llm,nvarid)
	327	#else
	328	ierr = NF_DEF_VAR (nid,"presnivs",NF_FLOAT,1,idim_llm,nvarid)
	329	#endif
	330	ierr = NF_ENDDEF(nid)
	331	#ifdef NC_DOUBLE
	332	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,presnivs)
	333	#else
	334	ierr = NF_PUT_VAR_REAL (nid,nvarid,presnivs)
	335	#endif
	336	c ------------------------------------------------------------------
	337	c Variable uniquement pour visualisation avec Grads ou Ferret
	338	c ------------------------------------------------------------------
	339	ierr = NF_REDEF (nid)
	340	#ifdef NC_DOUBLE
	341	ierr = NF_DEF_VAR (nid,"latitude",NF_DOUBLE,1,idim_rlatu,nvarid)
	342	#else
	343	ierr = NF_DEF_VAR (nid,"latitude",NF_FLOAT,1,idim_rlatu,nvarid)
	344	#endif
	345	ierr =NF_PUT_ATT_TEXT(nid,nvarid,'units',13,"degrees_north")
	346	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
	347	. "North latitude")
	348	ierr = NF_ENDDEF(nid)
	349	#ifdef NC_DOUBLE
	350	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu/pi*180)
	351	#else
	352	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu/pi*180)
	353	#endif
	354	c----------------------
	355	ierr = NF_REDEF (nid)
	356	#ifdef NC_DOUBLE
	357	ierr =NF_DEF_VAR(nid,"longitude", NF_DOUBLE, 1, idim_rlonv,nvarid)
	358	#else
	359	ierr = NF_DEF_VAR(nid,"longitude", NF_FLOAT, 1, idim_rlonv,nvarid)
	360	#endif
	361	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
	362	. "East longitude")
	363	ierr = NF_PUT_ATT_TEXT(nid,nvarid,'units',12,"degrees_east")
	364	ierr = NF_ENDDEF(nid)
	365	#ifdef NC_DOUBLE
	366	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv/pi*180)
	367	#else
	368	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv/pi*180)
	369	#endif
	370	c--------------------------
	371	ierr = NF_REDEF (nid)
	372	#ifdef NC_DOUBLE
	373	ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
	374	. idim_llm,nvarid)
	375	#else
	376	ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
	377	. idim_llm,nvarid)
	378	#endif
	379	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",10,"pseudo-alt")
	380	ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
	381	ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
	382
	383	ierr = NF_ENDDEF(nid)
	384	#ifdef NC_DOUBLE
	385	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
	386	#else
	387	ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
	388	#endif
	389
	390	!-------------------------------
	391	! (soil) depth variable mlayer() (known from comsoil.h)
	392	!-------------------------------
	393	ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
	394	! define variable
	395	#ifdef NC_DOUBLE
	396	ierr=NF_DEF_VAR(nid,"soildepth",NF_DOUBLE,1,idim_nsoilmx,nvarid)
	397	#else
	398	ierr=NF_DEF_VAR(nid,"soildepth",NF_FLOAT,1,idim_nsoilmx,nvarid)
	399	#endif
	400	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 20,
	401	. "Soil mid-layer depth")
	402	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",1,"m")
	403	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"positive",4,"down")
	404	ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
	405	! write variable
	406	#ifdef NC_DOUBLE
	407	ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,mlayer)
	408	#else
	409	ierr=NF_PUT_VAR_REAL (nid,nvarid,mlayer)
	410	#endif
	411
	412	!---------------------
	413	! soil thermal inertia
	414	!---------------------
	415	ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
	416	dims3(1)=idim_rlonv
	417	dims3(2)=idim_rlatu
	418	dims3(3)=idim_nsoilmx
	419	! define variable
	420	#ifdef NC_DOUBLE
	421	ierr=NF_DEF_VAR(nid,"inertiedat",NF_DOUBLE,3,dims3,nvarid)
	422	#else
	423	ierr=NF_DEF_VAR(nid,"inertiedat",NF_FLOAT,3,dims3,nvarid)
	424	#endif
	425	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 20,
	426	& "Soil thermal inertia")
	427	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",15,
	428	& "J.s-1/2.m-2.K-1")
	429	ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
	430	! write variable
	431	#ifdef NC_DOUBLE
	432	ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,ith)
	433	#else
	434	ierr=NF_PUT_VAR_REAL (nid,nvarid,ith)
	435	#endif
	436
	437	c-----------------------------------------------------------------------
	438	c Ecriture aire et coefficients de passage cov. <-> contra. <--> naturel
	439	c-----------------------------------------------------------------------
	440
	441	ierr = NF_REDEF (nid)
	442	dims2(1) = idim_rlonu
	443	dims2(2) = idim_rlatu
	444	#ifdef NC_DOUBLE
	445	ierr = NF_DEF_VAR (nid,"cu",NF_DOUBLE,2,dims2,nvarid)
	446	#else
	447	ierr = NF_DEF_VAR (nid,"cu",NF_FLOAT,2,dims2,nvarid)
	448	#endif
	449	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
	450	. "Coefficient de passage pour U")
	451	ierr = NF_ENDDEF(nid)
	452	#ifdef NC_DOUBLE
	453	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cu)
	454	#else
	455	ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)
	456	#endif
	457	c
	458	ierr = NF_REDEF (nid)
	459	dims2(1) = idim_rlonv
	460	dims2(2) = idim_rlatv
	461	#ifdef NC_DOUBLE
	462	ierr = NF_DEF_VAR (nid,"cv",NF_DOUBLE,2,dims2,nvarid)
	463	#else
	464	ierr = NF_DEF_VAR (nid,"cv",NF_FLOAT,2,dims2,nvarid)
	465	#endif
	466	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
	467	. "Coefficient de passage pour V")
	468	ierr = NF_ENDDEF(nid)
	469	#ifdef NC_DOUBLE
	470	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cv)
	471	#else
	472	ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)
	473	#endif
	474	c
	475	c Aire de chaque maille:
	476	c
	477	ierr = NF_REDEF (nid)
	478	dims2(1) = idim_rlonv
	479	dims2(2) = idim_rlatu
	480	#ifdef NC_DOUBLE
	481	ierr = NF_DEF_VAR (nid,"aire",NF_DOUBLE,2,dims2,nvarid)
	482	#else
	483	ierr = NF_DEF_VAR (nid,"aire",NF_FLOAT,2,dims2,nvarid)
	484	#endif
	485	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	486	. "Aires de chaque maille")
	487	ierr = NF_ENDDEF(nid)
	488	#ifdef NC_DOUBLE
	489	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aire)
	490	#else
	491	ierr = NF_PUT_VAR_REAL (nid,nvarid,aire)
	492	#endif
	493
	494	c-----------------------------------------------------------------------
	495	c Ecriture du geopentiel au sol
	496	c-----------------------------------------------------------------------
	497
	498	ierr = NF_REDEF (nid)
	499	dims2(1) = idim_rlonv
	500	dims2(2) = idim_rlatu
	501	#ifdef NC_DOUBLE
	502	ierr = NF_DEF_VAR (nid,"phisinit",NF_DOUBLE,2,dims2,nvarid)
	503	#else
	504	ierr = NF_DEF_VAR (nid,"phisinit",NF_FLOAT,2,dims2,nvarid)
	505	#endif
	506	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
	507	. "Geopotentiel au sol")
	508	ierr = NF_ENDDEF(nid)
	509	#ifdef NC_DOUBLE
	510	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,phis)
	511	#else
	512	ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)
	513	#endif
	514
	515	PRINT*,'iim,jjm,llm,idayref',iim,jjm,llm,idayref
	516	PRINT*,'rad,omeg,g,mugaz,kappa',
	517	s rad,omeg,g,43.49,kappa !mars temporaire (ecrire mugaz ensuite)
	518	PRINT*,'daysec,dtvr',daysec,dtvr
	519
	520	RETURN
	521	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format