Context Navigation

dynredem.F @ 1243

Last change on this file since 1243 was 1107, checked in by emillour, 11 years ago

Common dynamics: Updates and modifications to enable running Mars physics with

LMDZ.COMMON dynamics:

For compilation: adapted makelmdz, create_make_gcm and makelmdz_fcm, bld.cfg to compile aeronomy routines in "aerono$physique" if it exists, and added "-P -traditional" preprocessing flags in "arch-linux-ifort*"
Added function "cbrt.F" (cubic root) in 'bibio'
Adapted the reading/writing of dynamics (re)start.nc files for Mars. The main issue is that different information (on time, reference and current) is stored and used differently, hence a few if (planet_type =="mars") here and there. Moreover in the martian case there is the possibility to store fields over multiple times. Some Mars-specific variables (ecritphy,ecritstart,timestart) added in control_mod.F and (hour_ini) in temps.h

File size: 25.3 KB

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

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

Download in other formats:

Original Format