Context Navigation

dynredem.F @ 1033

Last change on this file since 1033 was 960, checked in by lsce, 16 years ago

Ajoute du parametre config_inca dans conf_gcm.F config_inca='none'(sans INCA, par defaut) config_inca='chem'(avec INCA config chemie) config_inca='aero'(avec INCA config aerosol)
Menage parmis les cles CPP INCA
Enleve le calcul d'omega dans calfis.F et active le calcul correspondant dans advtrac.F(avant uniquement pour INCA).

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 20.7 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	c
	5	SUBROUTINE dynredem0(fichnom,iday_end,phis,nq)
	6	USE IOIPSL
	7	IMPLICIT NONE
	8	c=======================================================================
	9	c Ecriture du fichier de redemarrage sous format NetCDF (initialisation)
	10	c=======================================================================
	11	c Declarations:
	12	c -------------
	13	#include "dimensions.h"
	14	#include "paramet.h"
	15	#include "comconst.h"
	16	#include "comvert.h"
	17	#include "comgeom.h"
	18	#include "temps.h"
	19	#include "ener.h"
	20	#include "logic.h"
	21	#include "netcdf.inc"
	22	#include "description.h"
	23	#include "serre.h"
	24	#include "advtrac.h"
	25
	26	c Arguments:
	27	c ----------
	28	INTEGER iday_end
	29	REAL phis(ip1jmp1)
	30	CHARACTER() fichnom
	31	INTEGER nq
	32
	33	c Local:
	34	c ------
	35	INTEGER iq,l
	36	INTEGER length
	37	PARAMETER (length = 100)
	38	REAL tab_cntrl(length) ! tableau des parametres du run
	39	INTEGER ierr
	40	character*20 modname
	41	character*80 abort_message
	42
	43	c Variables locales pour NetCDF:
	44	c
	45	INTEGER dims2(2), dims3(3), dims4(4)
	46	INTEGER idim_index
	47	INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
	48	INTEGER idim_s, idim_sig
	49	INTEGER idim_tim
	50	INTEGER nid,nvarid
	51
	52	REAL zan0,zjulian,hours
	53	INTEGER yyears0,jjour0, mmois0
	54	character*30 unites
	55
	56
	57	c-----------------------------------------------------------------------
	58	modname='dynredem'
	59
	60	call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
	61	call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
	62
	63
	64	DO l=1,length
	65	tab_cntrl(l) = 0.
	66	ENDDO
	67	tab_cntrl(1) = FLOAT(iim)
	68	tab_cntrl(2) = FLOAT(jjm)
	69	tab_cntrl(3) = FLOAT(llm)
	70	tab_cntrl(4) = FLOAT(day_ref)
	71	tab_cntrl(5) = FLOAT(annee_ref)
	72	tab_cntrl(6) = rad
	73	tab_cntrl(7) = omeg
	74	tab_cntrl(8) = g
	75	tab_cntrl(9) = cpp
	76	tab_cntrl(10) = kappa
	77	tab_cntrl(11) = daysec
	78	tab_cntrl(12) = dtvr
	79	tab_cntrl(13) = etot0
	80	tab_cntrl(14) = ptot0
	81	tab_cntrl(15) = ztot0
	82	tab_cntrl(16) = stot0
	83	tab_cntrl(17) = ang0
	84	tab_cntrl(18) = pa
	85	tab_cntrl(19) = preff
	86	c
	87	c ..... parametres pour le zoom ......
	88
	89	tab_cntrl(20) = clon
	90	tab_cntrl(21) = clat
	91	tab_cntrl(22) = grossismx
	92	tab_cntrl(23) = grossismy
	93	c
	94	IF ( fxyhypb ) THEN
	95	tab_cntrl(24) = 1.
	96	tab_cntrl(25) = dzoomx
	97	tab_cntrl(26) = dzoomy
	98	tab_cntrl(27) = 0.
	99	tab_cntrl(28) = taux
	100	tab_cntrl(29) = tauy
	101	ELSE
	102	tab_cntrl(24) = 0.
	103	tab_cntrl(25) = dzoomx
	104	tab_cntrl(26) = dzoomy
	105	tab_cntrl(27) = 0.
	106	tab_cntrl(28) = 0.
	107	tab_cntrl(29) = 0.
	108	IF( ysinus ) tab_cntrl(27) = 1.
	109	ENDIF
	110
	111	tab_cntrl(30) = FLOAT(iday_end)
	112	tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
	113	c
	114	c .........................................................
	115	c
	116	c Creation du fichier:
	117	c
	118	ierr = NF_CREATE(fichnom, NF_CLOBBER, nid)
	119	IF (ierr.NE.NF_NOERR) THEN
	120	WRITE(6,*)" Pb d ouverture du fichier "//fichnom
	121	WRITE(6,*)' ierr = ', ierr
	122	CALL ABORT
	123	ENDIF
	124	c
	125	c Preciser quelques attributs globaux:
	126	c
	127	ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 27,
	128	. "Fichier demmarage dynamique")
	129	c
	130	c Definir les dimensions du fichiers:
	131	c
	132	ierr = NF_DEF_DIM (nid, "index", length, idim_index)
	133	ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
	134	ierr = NF_DEF_DIM (nid, "rlatu", jjp1, idim_rlatu)
	135	ierr = NF_DEF_DIM (nid, "rlonv", iip1, idim_rlonv)
	136	ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
	137	ierr = NF_DEF_DIM (nid, "sigs", llm, idim_s)
	138	ierr = NF_DEF_DIM (nid, "sig", llmp1, idim_sig)
	139	ierr = NF_DEF_DIM (nid, "temps", NF_UNLIMITED, idim_tim)
	140	c
	141	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
	142	c
	143	c Definir et enregistrer certains champs invariants:
	144	c
	145	ierr = NF_REDEF (nid)
[690]	146	cIM 220306 BEG
	147	#ifdef NC_DOUBLE
	148	ierr = NF_DEF_VAR (nid,"controle",NF_DOUBLE,1,idim_index,nvarid)
	149	#else
[524]	150	ierr = NF_DEF_VAR (nid,"controle",NF_FLOAT,1,idim_index,nvarid)
[690]	151	#endif
	152	cIM 220306 END
[524]	153	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	154	. "Parametres de controle")
	155	ierr = NF_ENDDEF(nid)
	156	#ifdef NC_DOUBLE
	157	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
	158	#else
	159	ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
	160	#endif
	161	c
	162	ierr = NF_REDEF (nid)
[690]	163	cIM 220306 BEG
	164	#ifdef NC_DOUBLE
	165	ierr = NF_DEF_VAR (nid,"rlonu",NF_DOUBLE,1,idim_rlonu,nvarid)
	166	#else
[524]	167	ierr = NF_DEF_VAR (nid,"rlonu",NF_FLOAT,1,idim_rlonu,nvarid)
[690]	168	#endif
	169	cIM 220306 END
[524]	170	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
	171	. "Longitudes des points U")
	172	ierr = NF_ENDDEF(nid)
	173	#ifdef NC_DOUBLE
	174	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonu)
	175	#else
	176	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu)
	177	#endif
	178	c
	179	ierr = NF_REDEF (nid)
[690]	180	cIM 220306 BEG
	181	#ifdef NC_DOUBLE
	182	ierr = NF_DEF_VAR (nid,"rlatu",NF_DOUBLE,1,idim_rlatu,nvarid)
	183	#else
[524]	184	ierr = NF_DEF_VAR (nid,"rlatu",NF_FLOAT,1,idim_rlatu,nvarid)
[690]	185	#endif
	186	cIM 220306 END
[524]	187	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	188	. "Latitudes des points U")
	189	ierr = NF_ENDDEF(nid)
	190	#ifdef NC_DOUBLE
	191	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu)
	192	#else
	193	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu)
	194	#endif
	195	c
	196	ierr = NF_REDEF (nid)
[690]	197	cIM 220306 BEG
	198	#ifdef NC_DOUBLE
	199	ierr = NF_DEF_VAR (nid,"rlonv",NF_DOUBLE,1,idim_rlonv,nvarid)
	200	#else
[524]	201	ierr = NF_DEF_VAR (nid,"rlonv",NF_FLOAT,1,idim_rlonv,nvarid)
[690]	202	#endif
	203	cIM 220306 END
[524]	204	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
	205	. "Longitudes des points V")
	206	ierr = NF_ENDDEF(nid)
	207	#ifdef NC_DOUBLE
	208	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv)
	209	#else
	210	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv)
	211	#endif
	212	c
	213	ierr = NF_REDEF (nid)
[690]	214	cIM 220306 BEG
	215	#ifdef NC_DOUBLE
	216	ierr = NF_DEF_VAR (nid,"rlatv",NF_DOUBLE,1,idim_rlatv,nvarid)
	217	#else
[524]	218	ierr = NF_DEF_VAR (nid,"rlatv",NF_FLOAT,1,idim_rlatv,nvarid)
[690]	219	#endif
	220	cIM 220306 END
[524]	221	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	222	. "Latitudes des points V")
	223	ierr = NF_ENDDEF(nid)
	224	#ifdef NC_DOUBLE
	225	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatv)
	226	#else
	227	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv)
	228	#endif
	229	c
	230	ierr = NF_REDEF (nid)
[690]	231	cIM 220306 BEG
	232	#ifdef NC_DOUBLE
	233	ierr = NF_DEF_VAR (nid,"nivsigs",NF_DOUBLE,1,idim_s,nvarid)
	234	#else
[524]	235	ierr = NF_DEF_VAR (nid,"nivsigs",NF_FLOAT,1,idim_s,nvarid)
[690]	236	#endif
	237	cIM 220306 END
[524]	238	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 28,
	239	. "Numero naturel des couches s")
	240	ierr = NF_ENDDEF(nid)
	241	#ifdef NC_DOUBLE
	242	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,nivsigs)
	243	#else
	244	ierr = NF_PUT_VAR_REAL (nid,nvarid,nivsigs)
	245	#endif
	246	c
	247	ierr = NF_REDEF (nid)
[690]	248	cIM 220306 BEG
	249	#ifdef NC_DOUBLE
	250	ierr = NF_DEF_VAR (nid,"nivsig",NF_DOUBLE,1,idim_sig,nvarid)
	251	#else
[524]	252	ierr = NF_DEF_VAR (nid,"nivsig",NF_FLOAT,1,idim_sig,nvarid)
[690]	253	#endif
	254	cIM 220306 END
[524]	255	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 32,
	256	. "Numero naturel des couches sigma")
	257	ierr = NF_ENDDEF(nid)
	258	#ifdef NC_DOUBLE
	259	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,nivsig)
	260	#else
	261	ierr = NF_PUT_VAR_REAL (nid,nvarid,nivsig)
	262	#endif
	263	c
	264	ierr = NF_REDEF (nid)
[690]	265	cIM 220306 BEG
	266	#ifdef NC_DOUBLE
	267	ierr = NF_DEF_VAR (nid,"ap",NF_DOUBLE,1,idim_sig,nvarid)
	268	#else
[524]	269	ierr = NF_DEF_VAR (nid,"ap",NF_FLOAT,1,idim_sig,nvarid)
[690]	270	#endif
	271	cIM 220306 END
[524]	272	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
	273	. "Coefficient A pour hybride")
	274	ierr = NF_ENDDEF(nid)
	275	#ifdef NC_DOUBLE
	276	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ap)
	277	#else
	278	ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)
	279	#endif
	280	c
	281	ierr = NF_REDEF (nid)
[690]	282	cIM 220306 BEG
	283	#ifdef NC_DOUBLE
	284	ierr = NF_DEF_VAR (nid,"bp",NF_DOUBLE,1,idim_sig,nvarid)
	285	#else
[524]	286	ierr = NF_DEF_VAR (nid,"bp",NF_FLOAT,1,idim_sig,nvarid)
[690]	287	#endif
	288	cIM 220306 END
[524]	289	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
	290	. "Coefficient B pour hybride")
	291	ierr = NF_ENDDEF(nid)
	292	#ifdef NC_DOUBLE
	293	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bp)
	294	#else
	295	ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)
	296	#endif
	297	c
	298	ierr = NF_REDEF (nid)
[690]	299	cIM 220306 BEG
	300	#ifdef NC_DOUBLE
	301	ierr = NF_DEF_VAR (nid,"presnivs",NF_DOUBLE,1,idim_s,nvarid)
	302	#else
[524]	303	ierr = NF_DEF_VAR (nid,"presnivs",NF_FLOAT,1,idim_s,nvarid)
[690]	304	#endif
	305	cIM 220306 END
[524]	306	ierr = NF_ENDDEF(nid)
	307	#ifdef NC_DOUBLE
	308	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,presnivs)
	309	#else
	310	ierr = NF_PUT_VAR_REAL (nid,nvarid,presnivs)
	311	#endif
	312	c
	313	c Coefficients de passage cov. <-> contra. <--> naturel
	314	c
	315	ierr = NF_REDEF (nid)
	316	dims2(1) = idim_rlonu
	317	dims2(2) = idim_rlatu
[690]	318	cIM 220306 BEG
	319	#ifdef NC_DOUBLE
	320	ierr = NF_DEF_VAR (nid,"cu",NF_DOUBLE,2,dims2,nvarid)
	321	#else
[524]	322	ierr = NF_DEF_VAR (nid,"cu",NF_FLOAT,2,dims2,nvarid)
[690]	323	#endif
	324	cIM 220306 END
[524]	325	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
	326	. "Coefficient de passage pour U")
	327	ierr = NF_ENDDEF(nid)
	328	#ifdef NC_DOUBLE
	329	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cu)
	330	#else
	331	ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)
	332	#endif
	333	c
	334	ierr = NF_REDEF (nid)
	335	dims2(1) = idim_rlonv
	336	dims2(2) = idim_rlatv
[690]	337	cIM 220306 BEG
	338	#ifdef NC_DOUBLE
	339	ierr = NF_DEF_VAR (nid,"cv",NF_DOUBLE,2,dims2,nvarid)
	340	#else
[524]	341	ierr = NF_DEF_VAR (nid,"cv",NF_FLOAT,2,dims2,nvarid)
[690]	342	#endif
	343	cIM 220306 END
[524]	344	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
	345	. "Coefficient de passage pour V")
	346	ierr = NF_ENDDEF(nid)
	347	#ifdef NC_DOUBLE
	348	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cv)
	349	#else
	350	ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)
	351	#endif
	352	c
	353	c Aire de chaque maille:
	354	c
	355	ierr = NF_REDEF (nid)
	356	dims2(1) = idim_rlonv
	357	dims2(2) = idim_rlatu
[690]	358	cIM 220306 BEG
	359	#ifdef NC_DOUBLE
	360	ierr = NF_DEF_VAR (nid,"aire",NF_DOUBLE,2,dims2,nvarid)
	361	#else
[524]	362	ierr = NF_DEF_VAR (nid,"aire",NF_FLOAT,2,dims2,nvarid)
[690]	363	#endif
	364	cIM 220306 END
[524]	365	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
	366	. "Aires de chaque maille")
	367	ierr = NF_ENDDEF(nid)
	368	#ifdef NC_DOUBLE
	369	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aire)
	370	#else
	371	ierr = NF_PUT_VAR_REAL (nid,nvarid,aire)
	372	#endif
	373	c
	374	c Geopentiel au sol:
	375	c
	376	ierr = NF_REDEF (nid)
	377	dims2(1) = idim_rlonv
	378	dims2(2) = idim_rlatu
[690]	379	cIM 220306 BEG
	380	#ifdef NC_DOUBLE
	381	ierr = NF_DEF_VAR (nid,"phisinit",NF_DOUBLE,2,dims2,nvarid)
	382	#else
[524]	383	ierr = NF_DEF_VAR (nid,"phisinit",NF_FLOAT,2,dims2,nvarid)
[690]	384	#endif
	385	cIM 220306 END
[524]	386	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
	387	. "Geopotentiel au sol")
	388	ierr = NF_ENDDEF(nid)
	389	#ifdef NC_DOUBLE
	390	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,phis)
	391	#else
	392	ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)
	393	#endif
	394	c
	395	c Definir les variables pour pouvoir les enregistrer plus tard:
	396	c
	397	ierr = NF_REDEF (nid) ! entrer dans le mode de definition
	398	c
[690]	399	cIM 220306 BEG
	400	#ifdef NC_DOUBLE
	401	ierr = NF_DEF_VAR (nid,"temps",NF_DOUBLE,1,idim_tim,nvarid)
	402	#else
[524]	403	ierr = NF_DEF_VAR (nid,"temps",NF_FLOAT,1,idim_tim,nvarid)
[690]	404	#endif
	405	cIM 220306 END
[524]	406	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
	407	. "Temps de simulation")
	408	write(unites,200)yyears0,mmois0,jjour0
	409	200 format('days since ',i4,'-',i2.2,'-',i2.2,' 00:00:00')
	410	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "units", 30,
	411	. unites)
	412
	413	c
	414	dims4(1) = idim_rlonu
	415	dims4(2) = idim_rlatu
	416	dims4(3) = idim_s
	417	dims4(4) = idim_tim
[690]	418	cIM 220306 BEG
	419	#ifdef NC_DOUBLE
	420	ierr = NF_DEF_VAR (nid,"ucov",NF_DOUBLE,4,dims4,nvarid)
	421	#else
[524]	422	ierr = NF_DEF_VAR (nid,"ucov",NF_FLOAT,4,dims4,nvarid)
[690]	423	#endif
	424	cIM 220306 END
[524]	425	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
	426	. "Vitesse U")
	427	c
	428	dims4(1) = idim_rlonv
	429	dims4(2) = idim_rlatv
	430	dims4(3) = idim_s
	431	dims4(4) = idim_tim
[690]	432	cIM 220306 BEG
	433	#ifdef NC_DOUBLE
	434	ierr = NF_DEF_VAR (nid,"vcov",NF_DOUBLE,4,dims4,nvarid)
	435	#else
[524]	436	ierr = NF_DEF_VAR (nid,"vcov",NF_FLOAT,4,dims4,nvarid)
[690]	437	#endif
	438	cIM 220306 END
[524]	439	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
	440	. "Vitesse V")
	441	c
	442	dims4(1) = idim_rlonv
	443	dims4(2) = idim_rlatu
	444	dims4(3) = idim_s
	445	dims4(4) = idim_tim
[690]	446	cIM 220306 BEG
	447	#ifdef NC_DOUBLE
	448	ierr = NF_DEF_VAR (nid,"teta",NF_DOUBLE,4,dims4,nvarid)
	449	#else
[524]	450	ierr = NF_DEF_VAR (nid,"teta",NF_FLOAT,4,dims4,nvarid)
[690]	451	#endif
	452	cIM 220306 END
[524]	453	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 11,
	454	. "Temperature")
	455	c
	456	dims4(1) = idim_rlonv
	457	dims4(2) = idim_rlatu
	458	dims4(3) = idim_s
	459	dims4(4) = idim_tim
	460	IF(nq.GE.1) THEN
	461	DO iq=1,nq
[690]	462	cIM 220306 BEG
	463	#ifdef NC_DOUBLE
	464	ierr = NF_DEF_VAR (nid,tname(iq),NF_DOUBLE,4,dims4,nvarid)
	465	#else
[524]	466	ierr = NF_DEF_VAR (nid,tname(iq),NF_FLOAT,4,dims4,nvarid)
[690]	467	#endif
	468	cIM 220306 END
[524]	469	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))
	470	ENDDO
	471	ENDIF
	472	c
	473	dims4(1) = idim_rlonv
	474	dims4(2) = idim_rlatu
	475	dims4(3) = idim_s
	476	dims4(4) = idim_tim
[690]	477	cIM 220306 BEG
	478	#ifdef NC_DOUBLE
	479	ierr = NF_DEF_VAR (nid,"masse",NF_DOUBLE,4,dims4,nvarid)
	480	#else
[524]	481	ierr = NF_DEF_VAR (nid,"masse",NF_FLOAT,4,dims4,nvarid)
[690]	482	#endif
	483	cIM 220306 END
[524]	484	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,
	485	. "C est quoi ?")
	486	c
	487	dims3(1) = idim_rlonv
	488	dims3(2) = idim_rlatu
	489	dims3(3) = idim_tim
[690]	490	cIM 220306 BEG
	491	#ifdef NC_DOUBLE
	492	ierr = NF_DEF_VAR (nid,"ps",NF_DOUBLE,3,dims3,nvarid)
	493	#else
[524]	494	ierr = NF_DEF_VAR (nid,"ps",NF_FLOAT,3,dims3,nvarid)
[690]	495	#endif
	496	cIM 220306 END
[524]	497	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 15,
	498	. "Pression au sol")
	499	c
	500	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
	501	ierr = NF_CLOSE(nid) ! fermer le fichier
	502
	503	PRINT*,'iim,jjm,llm,iday_end',iim,jjm,llm,iday_end
	504	PRINT*,'rad,omeg,g,cpp,kappa',
	505	, rad,omeg,g,cpp,kappa
	506
	507	RETURN
	508	END
	509	SUBROUTINE dynredem1(fichnom,time,
	510	. vcov,ucov,teta,q,nq,masse,ps)
	511	IMPLICIT NONE
	512	c=================================================================
	513	c Ecriture du fichier de redemarrage sous format NetCDF
	514	c=================================================================
	515	#include "dimensions.h"
	516	#include "paramet.h"
	517	#include "description.h"
	518	#include "netcdf.inc"
	519	#include "comvert.h"
	520	#include "comgeom.h"
	521	#include "advtrac.h"
[541]	522	#include "temps.h"
[960]	523	#include "control.h"
[524]	524
	525	INTEGER nq, l
	526	REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm)
	527	REAL teta(ip1jmp1,llm)
	528	REAL ps(ip1jmp1),masse(ip1jmp1,llm)
	529	REAL q(ip1jmp1,llm,nq)
	530	CHARACTER() fichnom
	531
	532	REAL time
[822]	533	INTEGER nid, nvarid, nid_trac, nvarid_trac
	534	REAL trac_tmp(ip1jmp1,llm)
	535	INTEGER ierr, ierr_file
[524]	536	INTEGER iq
[541]	537	INTEGER length
	538	PARAMETER (length = 100)
	539	REAL tab_cntrl(length) ! tableau des parametres du run
[524]	540	character*20 modname
	541	character*80 abort_message
	542	c
	543	INTEGER nb
	544	SAVE nb
	545	DATA nb / 0 /
	546
	547	modname = 'dynredem1'
	548	ierr = NF_OPEN(fichnom, NF_WRITE, nid)
	549	IF (ierr .NE. NF_NOERR) THEN
	550	PRINT*, "Pb. d ouverture "//fichnom
	551	CALL abort
	552	ENDIF
	553
	554	c Ecriture/extension de la coordonnee temps
	555
	556	nb = nb + 1
	557	ierr = NF_INQ_VARID(nid, "temps", nvarid)
	558	IF (ierr .NE. NF_NOERR) THEN
	559	print *, NF_STRERROR(ierr)
	560	abort_message='Variable temps n est pas definie'
	561	CALL abort_gcm(modname,abort_message,ierr)
	562	ENDIF
	563	#ifdef NC_DOUBLE
	564	ierr = NF_PUT_VAR1_DOUBLE (nid,nvarid,nb,time)
	565	#else
	566	ierr = NF_PUT_VAR1_REAL (nid,nvarid,nb,time)
	567	#endif
	568	PRINT*, "Enregistrement pour ", nb, time
	569
[541]	570	c
	571	c Re-ecriture du tableau de controle, itaufin n'est plus defini quand
	572	c on passe dans dynredem0
	573	ierr = NF_INQ_VARID (nid, "controle", nvarid)
	574	IF (ierr .NE. NF_NOERR) THEN
	575	abort_message="dynredem1: Le champ <controle> est absent"
	576	ierr = 1
	577	CALL abort_gcm(modname,abort_message,ierr)
	578	ENDIF
	579	#ifdef NC_DOUBLE
	580	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
	581	#else
	582	ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
	583	#endif
	584	tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
	585	#ifdef NC_DOUBLE
	586	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
	587	#else
	588	ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
	589	#endif
	590
[524]	591	c Ecriture des champs
	592	c
	593	ierr = NF_INQ_VARID(nid, "ucov", nvarid)
	594	IF (ierr .NE. NF_NOERR) THEN
	595	PRINT*, "Variable ucov n est pas definie"
	596	CALL abort
	597	ENDIF
	598	#ifdef NC_DOUBLE
	599	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ucov)
	600	#else
	601	ierr = NF_PUT_VAR_REAL (nid,nvarid,ucov)
	602	#endif
	603
	604	ierr = NF_INQ_VARID(nid, "vcov", nvarid)
	605	IF (ierr .NE. NF_NOERR) THEN
	606	PRINT*, "Variable vcov n est pas definie"
	607	CALL abort
	608	ENDIF
	609	#ifdef NC_DOUBLE
	610	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,vcov)
	611	#else
	612	ierr = NF_PUT_VAR_REAL (nid,nvarid,vcov)
	613	#endif
	614
	615	ierr = NF_INQ_VARID(nid, "teta", nvarid)
	616	IF (ierr .NE. NF_NOERR) THEN
	617	PRINT*, "Variable teta n est pas definie"
	618	CALL abort
	619	ENDIF
	620	#ifdef NC_DOUBLE
	621	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,teta)
	622	#else
	623	ierr = NF_PUT_VAR_REAL (nid,nvarid,teta)
	624	#endif
	625
[960]	626	IF (config_inca /= 'none') THEN
[822]	627	! Ajout Anne pour lecture valeurs traceurs dans un fichier start_trac.nc
[960]	628	ierr_file = NF_OPEN ("start_trac.nc", NF_NOWRITE,nid_trac)
	629	IF (ierr_file .NE.NF_NOERR) THEN
	630	write(6,*)' Pb d''ouverture du fichier start_trac.nc'
	631	write(6,*)' ierr = ', ierr_file
	632	ENDIF
	633	END IF
[822]	634
[524]	635	IF(nq.GE.1) THEN
[960]	636	do iq=1,nq
	637
	638	IF (config_inca == 'none') THEN
	639	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
[822]	640	IF (ierr .NE. NF_NOERR) THEN
[960]	641	PRINT*, "Variable tname(iq) n est pas definie"
	642	CALL abort
	643	ENDIF
[822]	644	#ifdef NC_DOUBLE
[960]	645	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
[822]	646	#else
[960]	647	ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
[822]	648	#endif
[960]	649	ELSE ! config_inca = 'chem' ou 'aero'
	650	! lecture de la valeur du traceur dans start_trac.nc
	651	IF (ierr_file .ne. 2) THEN
	652	ierr = NF_INQ_VARID (nid_trac, tname(iq), nvarid_trac)
	653	IF (ierr .NE. NF_NOERR) THEN
	654	PRINT*, tname(iq),"est absent de start_trac.nc"
	655	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
	656	IF (ierr .NE. NF_NOERR) THEN
	657	PRINT*, "Variable ", tname(iq)," n est pas definie"
	658	CALL abort
	659	ENDIF
[822]	660	#ifdef NC_DOUBLE
[960]	661	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
	662	#else
	663	ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
	664	#endif
	665
	666	ELSE
	667	PRINT*, tname(iq), "est present dans start_trac.nc"
	668	#ifdef NC_DOUBLE
[822]	669	ierr = NF_GET_VAR_DOUBLE(nid_trac, nvarid_trac, trac_tmp)
	670	#else
	671	ierr = NF_GET_VAR_REAL(nid_trac, nvarid_trac, trac_tmp)
	672	#endif
[960]	673	IF (ierr .NE. NF_NOERR) THEN
	674	PRINT*, "Lecture echouee pour", tname(iq)
	675	CALL abort
	676	ENDIF
	677	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
	678	IF (ierr .NE. NF_NOERR) THEN
	679	PRINT*, "Variable ", tname(iq)," n est pas definie"
	680	CALL abort
	681	ENDIF
[822]	682	#ifdef NC_DOUBLE
[960]	683	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,trac_tmp)
[822]	684	#else
[960]	685	ierr = NF_PUT_VAR_REAL (nid,nvarid,trac_tmp)
[822]	686	#endif
[960]	687
	688	ENDIF ! IF (ierr .NE. NF_NOERR)
[822]	689	! fin lecture du traceur
[960]	690	ELSE ! si il n'y a pas de fichier start_trac.nc
[822]	691	! print *, 'il n y a pas de fichier start_trac'
	692	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
	693	IF (ierr .NE. NF_NOERR) THEN
	694	PRINT*, "Variable tname(iq) n est pas definie"
	695	CALL abort
	696	ENDIF
	697	#ifdef NC_DOUBLE
	698	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
	699	#else
	700	ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
	701	#endif
[960]	702	ENDIF ! (ierr_file .ne. 2)
	703	END IF ! config_inca
	704
[524]	705	ENDDO
	706	ENDIF
	707	c
	708	ierr = NF_INQ_VARID(nid, "masse", nvarid)
	709	IF (ierr .NE. NF_NOERR) THEN
	710	PRINT*, "Variable masse n est pas definie"
	711	CALL abort
	712	ENDIF
	713	#ifdef NC_DOUBLE
	714	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,masse)
	715	#else
	716	ierr = NF_PUT_VAR_REAL (nid,nvarid,masse)
	717	#endif
	718	c
	719	ierr = NF_INQ_VARID(nid, "ps", nvarid)
	720	IF (ierr .NE. NF_NOERR) THEN
	721	PRINT*, "Variable ps n est pas definie"
	722	CALL abort
	723	ENDIF
	724	#ifdef NC_DOUBLE
	725	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ps)
	726	#else
	727	ierr = NF_PUT_VAR_REAL (nid,nvarid,ps)
	728	#endif
	729
	730	ierr = NF_CLOSE(nid)
	731	c
	732	RETURN
	733	END
	734

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

Download in other formats:

Original Format